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ABSTRACT 


Acoustic  instability  in  solid  propellant  rockets  is  observed  as  a 
high-amplitude  oscillation  in  the  acoustic  modes  of  the  gas  cavity,  the 
combustion  processes  providing  the  driving  energy.  In  general,  the 
solid  phase,  l.e.,  the  propellant,  participates  in  the  wave  behavior. 

In  this  work,  the  behavior  of  both  the  gas  phase  and  the  solid  phase 
were  experimentally  studied  in  a  small  rocket  motor. 

The  motor  used  in  this  study  is  known  as  a  T-burner  and  is, 
basically,  a  centrally- vented  cylindrical  cavity  with  propellant  being 
burned  at  one  or  both  ends.  In  such  a  motor  the  longitudinal  modes  of 
oscillation  are  excited.  The  T-burner  has  been  fairly  widely  used  to 
estimate  the  amplifying  characteristics,  in  the  form  of  acoustic  admit¬ 
tances,  of  the  burning  surface.  In  the  gas-phase  portion  of  the  author's 
woi*k,  pressure  measurements  were  made  at  both  ends  of  a  T-burner  while 
a  thin  sample  of  propellant  was  burned  in  one  end.  In  the  solid-phase 
work,  similar  measurements  were  made  but  a  long  cylinder  of  propellant 
was  burned,  the  propellant  being  advanced  as  it  burned  to  maintain  the 
burning  interface  at  a  fixed  position.  The  propellant  used  was  com¬ 
posed  principally  of  ammonium  perchlorate  and  polybutadiene-acrylic 
acid  copolymer. 

The  data  from  the  gas-phase  study  showed  evidence  of  unusual 
acoustic  behavior,  unusual  in  the  sense  that  it  was  considerably  dif¬ 
ferent  from  previous  descriptions.  In  the  past  the  acoustic  behavior 
in  the  gas  cavity  has  been  thought  to  closely  approximate  the  behavior 
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that  would  be  present  in  the  absence  of  gains,  losses,  or  mean  flow. 

Such  behavior  requires  the  pressure  amplitudes  to  be  nearly  equal,  at 
both  ends  of  the  motor.  The  data  obtained  showed  the  amplitude  at 
the  burning  end  of  the  motor  to  be  lower  by  a  factor  of  .  nominally. 

A  careful  evaluation  of  the  experimental  system  verified  the  observation. 
Subsequent  consideration  led  to  the  conclusion  that  an  asymmetric 
reflection  was  occurring  in  the  vicinity  of  the  vent,  the  reflection 
being  such  that  a  portion  of  the  wave  traveling  toward  the  propellant 
was  reflected  so  as  to  reinforce  the  wave  traveling  in  the  opposite 
direction.  There  appears  to  be  no  other  way  of  explaining  the  data. 

The  usual  method  of  using  the  T-burner  is  to  measure  rates  of 
growth  and  decay  of  oscillations  and  subsequently  infer  values  for 
the  acoustic  admittance  from  these  rates .  Such  measurements  were  also 
made  in  the  author's  work.  In  addition,  a  technique  was  devised  to 
extend  the  method  in  the  low  frequency  region  where  the  oscillations 
would  not  grow,  i.e.,  the  motor  was  stable.  Firings  were  made  with  a 
nominal  frequency  of  200  cps  using  a  T-burner.  It  was  found  that  the 
pressure  pulse  from  the  igniter  system  excited  oscillations,  but  the 
combustion  provided  insufficient  energy  to  maintain  the  oscillations. 
Decay  rates  were  measured  with  propellant  burning  in  one,  both,  or 
neither  end  of  the  motor.  These  decay  rates  wore  used  to  estimate  the 
acoustic  admittance  of  the  burning  surface. 

In  the  study  of  solid  phase  participation  in  the  oscillations, 
the  pressure  measurement  at  the  non-burning  end  of  the  motor  was  used 
to  infer  the  pressure  at  the  burning  surface.  The  inferred  value  and 


x 


the  measurement  under  the  propellant  were  subsequently  used  to  infer 
the  acoustic  behavior  of  the  propellant.  Initially,  the  data  appeared 
to  indicate  a  degradation  of  the  propellant  with  burning  time.  However, 
further  examination  showed  that  the  apparent  changes  could  equally  well 
be  due  to  wall  effects  or  changes  in  the  amount  of  asymmetric  reflection 
during  the  burning.  The  results  were  rather  inconclusive  because  of 
the  inability  to  separate  the  various  effects. 


CHAPTER  I 


INTRODUCTION 


It  seems  appropriate  to  begin  with  a  description  of  what  is  meant 
by  acoustic  instability.  If  a  rocket  motor  is  exhibiting  acoustic 
instability,  it  Is  observed  that  the  acoustic  modes  of  oscillation 
for  the  gas  cavity  are  excited  to  a  high  amplitude.  More  specifically, 
if  one  were  to  make  multiple  pressure  measurements  in  the  gas  cavity 
of  a  motor  during  a  period  of  unstable  but  steady  operation,  one 
would  measure  a  pressure  of  the  type 


p  *  P0  +  pi 

where  p1  is  of  the  type 

P1  -  f(x,y,  z)  cos  foot  +  0{x,y,  z)] 


(1.1) 


(1.2) 


The  nomenclature  may  be  found  in  Appendix  0.  The  observed  variation 
of  the  oscillating  pressure,  p^,  with  position  and  the  observed  variation 
of  frequency  with  the  dimensions  of  the  cavity  would  be  well  described 
by  the  corresponding  solutions  of  the  acoustic  wave  equation  for  the 
particular  geometry  [e.g.,  17,21] .  Apparently  Smith  and  Sprenger  were 
among  the  first  to  observe  this  relationship  [52).  Because  of  the 
ever-present  damping  in  all  real  systems,  the  pressure  oscillations 
in  an  unstable  motor  would  decay  out  if  no  acoustic  energy  source  were 
present.  The  combustion  process  provides  a  source.  The  combustion  of 
solid  propellants  at  normal  operating  pressures  takes  place  in  a  very 
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thin  zone,  the  solid  being  heated  sufficiently  for  gasification  and 
burned  within  a  region  about  100  microns  in  thickness  [e.g.,  21,56]. 
Leaving  this  region  are  gaseous  products  which  are  very  nearly  in 
chemical  equilibrium.  It  is  this  region,  the  flame  or  combustion  zone, 
that  appears  to  provide  the  most  likely  source  of  acoustic  energy; 
and  in  fact,  it  is  believed  to  be  the  primary  source  [e.g.,  5,21,22, 
23,36].  The  only  other  possible  sources  are  expected  to  be  unimportant, 
being  due  vortex  generation  or  reactions  in  the  main  body  of  the  gas. 

In  a  practical  motor  acoustic  instability  may  lead  to  many 
undesirable  effects  [e.g.,  1,21,43].  At  times  the  pressure  amplitude 
may  become  very  large.  Amplitudes  of  the  same  order  of  magnitude  as 
the  mean  pressure,  with  mean  pressure  levels  of  70  or  100  atm.,  have 
been  reported.  Occasionally  the  propellant  breaks  up  due  to  very  large 
stresses,  and  excessive  mean  pressures  result.  Burning  rate  changes 
frequently  occur.  Large  roll  torques  have  resulted.  Heat  transfer 
to  the  hardware  is  enhanced.  Even  at  low  amplitudes,  vibrations  are 
transmitted  throughout  the  motor.  Obviously  these  phemomena  can  lead 
to  inefficient  operation  and  even  disastrous  results. 

Acoustic  instability  has  been  an  ever-increasing  problem  since 
about  1942.  As  the  propellants  have  become  more  energetic,  so  too 
have  the  instances  of  instability.  In  earlier  years,  and  in  some  engines 
even  today,  such  things  as  inert  baffles  or  rods,  and  holes  in  the 
propellant  were  used  as  remedies  for  the  problem  [e.g.,  1,30].  In  every 
case  the  remedy  was  obtained  by  a  trial-end-error  process  because  so 
little  was  known  about  the  mechanism  of  the  instability.  More  recently, 
adding  powdered  metals  to  the  propellant  has  proven  very  effective  in 
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preventing  or  reducing  the  level  of  instability.  There  is  no  assurance, 
however,  that  powdered  metals  will  be  adequate  with  new  and  exotic 
propellants  of  the  future.  The  instability  problem  has  proven  so 
troublesome  that  there  are  many  laboratories  devoting  considerable 
effort  to  various  aspects  of  the  problem. 

Other  Research:  Current  and  Past 

A  number  of  theoretical  analyses  of  the  combustion  zone  during 
acoustic  instability  have  been  presented.  Because  of  the  extreme 
complexity  of  the  behavior  in  the  combustion  zone,  each  theory  has 
treated  a  simplified  model.  The  first  theory  put  forth  was  by  Grad 
in  19^  [19  ]•  His  theory  included  a  kind  of  relaxation  time  or  time 
lag  that  was  to  be  associated  with  finite  reaction  rates.  Also  the 
later  theories  of  Cheng  [8]  and  Green  [20]  included  time  lags.  Each 
of  these  models  predicted  a  very  peaked  frequency  response  and  thus 
indicated  instability  in  selected  frequency  ranges .  Experimental 
data  which  are  now  available  show  clearly  that  this  is  not  the  case, 
oscillations  being  observed  over  wide  ranges  of  frequency. 

Other  models  have  been  analyzed  and  presented  more  recently 
which  do  not  incorporate  a  time  lag,  for  example  the  theories  of  Smith 
[ 51 of  Shinnar  and  Dishon  [50],  and  of  McClure,  Hart,  and  Bird 

[2,2k,  36]. 

The  most  prominent  and  complete  theory  is  that  of  McClure,  Hart, 
and  Bird.  They  have  treated  a  simplified  model  also;  however,  they 
have  gone  to  considerable  lengths  to  retain  as  much  realism  as  possible 
and  still  have  a  tractable  analysis.  Their  model  for  the  burning  zone 
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was  one  dimensional  and  all  reactions  were  assumed  to  be  localized  at 
two  planes,  the  surface  and  the  flame.  Diffusion  effects  were  neglected, 
as  were  momentum  effects.  Even  with  a  model  this  restricted,  the  solutions 
to  the  pertinent  equations  became  very  complex.  The  final  solutions  con¬ 
tain  a  good  many  (seventeen)  physical  parameter.-  which  are,  for  the  most 
part,  unknown  for  most  propellants.  By  means  of  what  is  hoped  to  be 
judicious  guessing  at  the  values  of  these  parameters,  they  predict 
instability  over  a  broad  range  of  frequencies.  This  is  held  to  be  a 
reasonable  degree  of  success  for  the  model. 

Also  McClure  and  associates  demonstrated  that  one  should  expect 
the  solid  propellant  to  participate  in  the  acoustic  behavior  [2 ] . 

Angelus  was  able  to  experimentally  verify  this  kind  of  behavior  [l]. 

The  participation  of  the  solid  phase  was  also  dramatically  demonstrated 
by  Ryan  et  al .  [47]. 

Experimentally,  there  has  been  considerable  effort  expended  in 
recent  years  in  attempting  to  measure  various  contributions  to  the 
overall  instability  picture.  This  is  demonstrated  by  the  large  number 
of  people  engaged  in  instability  research  [see  for  instance,  references 
31,57*58,  also  the  June  and  July,  1964,  issues  of  the  AIAA  Journal]. 

There  are  many  review  articles  available  which  discuss  the  current 
status  of  both  theoretical  and  experimental  research.  The  most  current 
are  those  by  Hart  and  McClure  [21]  on  theory  and  Price  [44]  on  experi¬ 
mented  research  which  were  presented  at  the  Tenth  Symposium  on  Combustion. 
Also  references  [2,43]  are  good  review  articles. 

The  theory  of  the  unstable  behavior  in  the  flame  zone  has  advanced 
to  such  a  stage  that  many  of  the  general  features  of  experimental  findings 
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may  be  explained  [37l»  On  the  other  hand,  it  is  at  a  stage  where 
little  more  can  be  learned  without  detailed  and  reasonably  accurate 
experimental  information.  The  experimental  information  has  been  slow 
in  becoming  available  and  is  currently  rather  limited.  The  reason 
for  the  lack  of  adequate  experimental  information  is  the  great  diffi¬ 
culty  in  making  the  necessary  measurements. 

Both  theory  and  experiment  have  indicated  that  the  acoustic 
behavior  of  the  flame  zone  may  be  treated  as  a  boundary  condition  for 
the  acoustic  problem  in  the  gas  cavity  [e.g.,  2,21].  Theoretical 
models  may  be  related  to  this  boundary  condition  in  order  to  evaluate 
their  success.  From  the  practical  standpoint,  the  boundary  condition 
must  be  used  to  assess  the  stability  of  a  practical  motor  [21].  In 
acoustic  terms  the  boundary  condition  may  be  characterized  by  the 
acoustic  admittance,  y,  which  is  defined  as 


Consequently  it  would  be  desirable  to  directly  measure  the  acoustic 
velocity  and  acoustic  pressure  in  the  flame  zone  so  as  to  directly 
obtain  the  admittance.  Unfortunately  the  flame  zone  is  hostile  and 
resists  such  measurements  because  of  the  small  thickness  (—  100  microns), 
and  the  high  mean  temperature  (5»2500°K)  and  pressure(:»  15  atm.).  The 
acoustic  behavior  of  the  flame  is  not  easily  measured  either,  as  the 
amplification  that  an  impinging  wave  receives  upon  reflection  is  of  the 
order  of  one  per  cent  [2,21].  Consequently,  the  appropriate  measure¬ 
ments  are  very  difficult  and  the  ingenuity  of  the  experimentalist  is 


greatly  taxed. 
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Several  kinds  of  experiments  have  been  tried  in  order  to  measure 
the  admittance.  For  instance,  several  attempts  have  been  made  to 
measure  the  effect  of  the  combustion  on  a  single  reflected  pulse  [32, 

40, 4l].  The  idea  is  to  compare  the  shape  of  a  pulse  or  weak  shock 
wave  as  it  impinges  on  the  burning  surface  with  the  shape  of  the  same 
pulse  after  reflection.  To  date,  no  positive  results  have  been  produced 
although  great  care  and  diligence  have  been  exhibited  by  the  experi¬ 
mentalists.  Similarly,  no  results  have  been  published  from  attempts 
which  involved  passing  a  wave  through  the  flame  from  the  propellant 
side  [4].  Sirens  have  been  used  to  create  an  oscillating  environment, 
but  here  again,  no  admittance  values  were  produced.  Similarly,  a 
modification  of  a  conventional  acoustic  approach,  the  Mawardt  Method, 
has  failed  [ 33,  34 ] .  Only  a  device  known  as  a  T-burner  has  had  a 
reasonable  degree  of  success.  A  T-burner  was  used  in  this  investigation. 

T- Burners 

A  typical  T-burne -  if  .•icted  in  Figure  1.  Basically  it  is  a 
cylindrical  rocket  motor  in  *  icn  p  -opellant  samples  may  be  burned  in 
one  or  both  ends.  The  combustion  products  are  allowed  to  escape 
through  a  centrally  located  vent  For  convenience  we  shall  describe 
the  burner  as  a  "single-"  or  "double-ended"  burner  when  propellant  is 
burned  in  one  or  both  ends,  respectively.  Horton,  while  at  the  Uni¬ 
versity  of  Utah,  showed  that  such  a  motor  would  exhibit  instability  for 
sufficiently  unstable  propellants  [26].  In  this  configuration  the 
longitudinal  modes  of  oscillation  are  excited,  and  for  the  most  part, 
only  the  fundamental  mode  is  observed.  Following  ignition  of  the 
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Figure  1.  Sketch  of  T-burner  configuration 
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propellant  the  oscillatory  pressure  grows  in  an  exponential  fashion, 

1  •  e  • , 

p1  =  A  cos  cot  (l.4) 

The  quantity  a  is  known  as  the  growth  constant.  It  is  experimentally 
determined  by  plotting  the  logarithm  of  the  oscillatory  pressure  ampli¬ 
tude  against  time  and  measuring  the  slope.  After  the  propellant  is 
consumed  the  oscillatory  pressure  decays  in  a  similar  exponential  fashion. 
The  corresponding  a  is  known  as  the  decay  constant.  These  growth  and 
decay  constants  have  been  used  to  infer  values  of  the  acoustic  admittance. 
The  usual  methods  for  using  these  data  are  described  by  Coates  et  al  - 
[10,  also  Appendix  E]. 

In  Figure  2  are  shown  oscilloscope  traces  which  typify  the  kinds 
of  pressure  records  obtained  from  a  T-burner.  The  upper  trace  is  the 
oscillating  component  of  pressure  and  the  lower  trace  is  the  total 
pressure.  This  particular  record  was  obtained  from  a  T-burner  whose 
exhaust  vent  was  a  sonic  nozzle;  the  equilibrium  mean  pressure  was, 
therefore,  determined  by  the  nozzle  size.  Frequently  T-burners  are 
vented  into  a  pre-pressurized  surge  tank  so  that  there  is  little  change 
in  mean  pressure  during  ignition  and  burning.  In  Figure  3  is  shown  a 
typical  oscillatory  growth  transient. 

Introduction  to  this  Investigation 

This  investigation  is  divided  into  two  rather  distinct  parts, 
the  first  concerned  with  the  solid  phase  and  the  second  with  the  gas 
phase.  The  solid  phase  work  developed  from  experimental  work  by  Coates  [11 ] 
and  an  analysis  of  his  data  by  Ryan  and  Coates  [49,  also  see  Appendix  L]. 
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Figure  2.  Typical  trace  from  a  single-ended  T-burner  firing. 
Run  No.  O-255.  Top  trace  Is  oscillatory  pressure  followed  by 
reference  at  .27  atm/cm.  Bottom  trace  is  total  pressure  at 
6.8  atm/cm.  Sweep  rate  is  .5  sec/cm.  from  left. 


Figure  3«  Typical  oscillatory  growth  transient.  Run  No.  0-251. 
Sweep  rate  100  msec/cm.  from  left.  No  vertical  scale. 
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Coates '  experiment  made  use  of  a  variation  of  the  T-burner  in 
which  one  fairly  long  grain  of  propellant  was  burned.  This  grain  was 
advanced  during  burning  at  a  rate  commensurate  with  its  burning  rate 
so  that  the  burning  interface  remained  at  a  nearly  fixed  position. 

Coates'  data  showed  rather  dramatic  evidence  of  the  participation  of 
the  propellant  in  the  oscillatory  behavior  [^7]*  In  addition  the 
analysis  of  his  data  suggested  that  the  viscoelastic  properties  of 
propellant  were  drastically  deteriorating  while  the  propellant  oscillated. 
Coates'  experiment  was  repeated  in  this  investigation  in  order  to  study 
the  apparent  degradation  in  viscoelastic  properties.  Further  discussion 
of  this  part  of  the  work  will  be  deferred  until  later. 

The  second  part  of  this  work  developed  quite  naturally  from  the 
first.  In  order  to  understand  the  behavior  of  the  solid  propellant 
it  was  necessary  to  know  the  oscillating  pressure  at  both  ends  of  the 
propellant.  The  pressure  at  the  burning  surface  had  been  inferred 
from  pressure  measurements  at  the  non-burning  end  of  the  motor  through 
acoustic  theory.  However,,  an  experimental  check  of  the  approach  showed 
the  pressure  inferred  was  greatly  in  error.  The  second  part  of  this 
work  became  one  of  trying  to  ascertain  why  the  oscillating  pressure 
distribution  was  much  different  than  expected  on  the  basis  of  previously 
reported  information.  The  experimental  data  suggested  that  the  acoustic 
behavior  of  the  burning  surface  was  possibly  quite  different  from  any¬ 
thing  previously  report  The  implications  were  sufficiently  important 
to  warrant  a  thorough  investigation.  In  addition  a  portion  of  the  gas- 
phase  work  was  concerned  with  obtaining  the  acoustic  admittance  values 
for  several  propellants. 


CHAPTER  II 


GAS  PHASE  THEORY 

The  acoustic  behavior  in  the  gas  cavity  of  a  T-burner  is  usually 
thought  to  approximate  the  acoustic  behavior  of  a  rigid,  closed  tube 
which  contains  an  inviscid,  non-reacting  gas  with  the  same  velocity  of 
sound  as  the  combusion  products.  For  convenience  we  shall  describe 
this  no-loss -no-gain  acoustic  field  as  the  "zero-order"  acoustic  field. 
The  fundamental,  longitudinal  mode  of  oscillation  of  the  gas  cavity 
in  a  T-burner  is  thus  expected  to  be  approximated  by 

p  -  A  <>*  =08  2S  (2.1) 

c 

where  the  origin  of  coordinates  has  been  taken  at  one  end  of  the  motor 
and  the  radial  dependence  has  been  neglected.  In  particular  if  one 
were  to  measure  the  acoustic  pressure  at  both  ends  of  a  T-burner,  he 
would  expect  to  find  equal  amplitudes  and  a  phase  difference  of  «  radians. 
In  this  work  preliminary  experiments  had  shown  the  pressure  amplitudes 
at  the  two  ends  of  a  single-ended  T-burner  to  be  substantially  different 
rather  them  nearly  equal.  The  primary  concern  of  this  study  was  to 
ascertain  why  the  amplitudes  were  unequal.  If  It  was  due  to  the 
combustion  processes,  then  it  would  be  of  great  importance  since  previous 
estimates  of  the  acoustic  admittance  for  the  burning  surface  could  be  in 
considerable  error.  If  it  was  due  to  other  acoustic  phenomena,  it  would 
at  least  be  valuable  information.  Rather  fortunately,  the  latter  will 
be  shown  to  be  the  case. 
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The  approach  used  in  this  work  was  to  develop  an  analytical  des¬ 
cription  a’  the  acoustic  behavior  in  terms  of  several  parameters.  These 
parameters  were  then  evaluated  from  experimental  data.  The  analysis 
was  sufficient  to  describe  the  effects  of  combustion,  if  these  were  the 
source  of  the  unequal  amplitudes,  but  was  insufficient  to  describe 
strong  reflection  phenomena.  The  acoustic  field  defined  by  the  experi¬ 
mentally-determined  parameters  was  then  tested  for  its  ability  to  describe 
the  real  behavior.  In  so  doing,  sufficient  information  was  obtained  to 
determine  that  the  deviations  from  an  approximately  zero-order  field 
were,  in  all  likelihood,  caused  by  an  asymmetric  reflection  in  the  vicinity 
of  the  nozzle. 

We  shall  proceed  to  develop  a  theoretical  analysis  of  the  acoustic 
behavior  in  the  gas  cavity  as  well  as  other  needed  relationships. 

Although  the  analytical  description  will  actually  prove  inadequate  to 
describe  the  real  behavior,  the  analysis  will  prove  useful  in  elucidating 
the  situation. 

Qualitative  Description  of  Motor  Behavior 

Before  embarking  on  an  analytical  description  of  the  acoustic 
behavior  in  a  T-burner,  the  complexity  of  the  system  will  be  noted  in 
a  qualitative  manner.  Although  the  T-burner  was  originally  conceived 
as  a  system  which  would  be  relatively  simple  to  describe  analytically, 
the  system  is  nevertheless  still  very  complex.  In  Figure  k  is  a  sketch 
of  a  T-burner;  the  dashed  lines  are  meant  to  depict  the  mean  flow.  Gen 
is  generated  by  the  burning  propellant  and  subsequently  flows  out  through 
a  sonic  nozzle.  In  the  vicinity  of  the  nozzle  the  mean  velocity  must  be 


BURNING 

PROPELLANT 


Figure  4.  Sketch  of  T-burner  configuration  shoving  (expected) 
lines  of  mean  flow. 
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large.  At  the  non-burning  or  "cold"  end  the  gas  is  expected  to  be 
relatively  stagnant  except  for  some  circulation  due  to  the  mean  flow 
near  the  nozzle.  This  circulation  should  be  adequate  to  maintain  the 
mean  gas  temperature  nearly  rniform  throughout  the  motor.  Following 
ignition  there  may  be  some  cold  gas  trapped  in  the  non- burning  end, 
but  it  is  eventually  swept  out  by  the  circulation.  At  the  other  end 
there  will  be  a  boundary  layer  at  the  curved  walls  which  will  build  up 
from  the  burning  surface  in  a  manner  similar  to  the  entrance  region  in 
a  tube.  The  wall  temperature  near  the  propellant  will  be  greater  than 
farther  down  the  cylinder.  There  is  possibly  a  considerable  degree 
of  small-scale  turbulence  in  the  vicinity  of  the  burning  surface 
because  of  the  heterogeneous  nature  of  the  propellant. 

It  is  readily  seen  that  a  detailed  mechanistic  description  is 
impossible;  simplifying  assumptions  will  be  required.  In  the  past  it 
has  frequently  been  assumed  that  the  acoustic  behavior  in  the  gas 
cavity  of  a  rocket  motor  deviates  only  slightly  from  the  behavior  that 
would  occur  in  the  absence  of  combustion  and  transpiration  [e.g.,  3^]. 
Similarly,  it  has  been  assumed  that  the  acoustic  behavior  in  a  T- 
burner  is  very  nearly  the  same  as  that  of  a  closed  cylinder  [e.g.,  10, 
55].  The  effects  of  the  mean  flow  have  either  been  neglected  entirely 
or  have  been  assumed  to  modify  only  the  boundai’y  conditions  [e.g.,  39]  • 
At  times,  the  one-dimens ional  damped-wave  equation  has  been  applied, 
with  the  effects  of  mean  flow  being  neglected  [e.g.,  10,55]*  Thus,  it 
may  be  observed  that  previous  theoretical  descriptions  assume  that  the 
behavior  in  the  main  body  of  the  gas  is  reasonably  well  described  by 
simple  acoustics. 
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We  will  now  consider  some  of  the  regions  of  particular  analytical 
difficulty. 

1.  The  nozzle  is  particularly  difficult  to  handle.  There  has 
been  some  theoretical  and  experimental  work  done  on  acoustic 
losses  associated  with  vents.  Hart  [25]  and  McClure  [39] 
have  concluded  from  their  theoretical  studies  that  the  losses 
are  negligible  for  a  subsonic  orifice.  Cold-flow  experi¬ 
ments  carried  out  at  Bolt,  Beranek,  and  Newman,  Inc.  [l8] 
indicate  that  the  losses  are  approximately  proportional  to 
the  vent  cross-sectional  area  and  inversely  proportional 

to  the  Mach  Number  in  the  vent,  at  least  for  low  velocities 
(subsonic)  in  the  vent.  Also  it  is  strongly  dependent  upon 
displacement  from  the  nodal  position,  i.e.,  the  position  of 
the  pressure  node.  These  studies  are  inadequate  for  the 
problem  at  hand  but  they  do  permit  some  hope  that  the  nozzle 
may  be  neglected  without  serious  error.  In  the  literature 
there  appears  to  be  no  consideration  of  major  distortion  of 
the  acoustic  field  attributable  to  the  nozzle.  Experimentally 
the  overall  losses  in  the  T-burner  used  in  this  study  have  not 
been  found  to  be  highly  sensitive  to  nozzle  location  so  that 
these  effects  would  appear  to  be  reasonably  small.  In  any 
case  the  specific  effects  of  the  nozzle  will  be  neglected  in 
this  analysis. 

2.  The  mean  flow  can  usually  be  neglected  in  acoustic  analyses 
as  long  as  the  Mach  Number  is  small.  The  average  Mach  Number 
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for  the  system  used  was  about  .01  and  should  be  negligible. 

The  effects  of  mean  flow  will  be  estimated. 

3«  The  damping  due  to  the  walls  is  dependent  upon  the  wall 

temperature  so  that  the  damping  should  be  position  dependent. 
However,  the  wall  damping  with  a  steep  thermal  gradient  and 
uniform  wall  temperature  is  very  difficult  to  treat 
analytically  [6j  and,  therefore,  one  can  only  hope  that  an 
average  value  for  the  damping  is  adequate. 

h.  The  effect  of  the  turbulence  is  unknown;  however,  it  is  pre¬ 
sumed  that  it  may  be  neglected  because  the  scale  of  turbulence 
is  much  smaller  than  the  acoustic  wave  length. 

5.  The  effect  of  a  mean  temperature  variation  along  the  axiB 

of  the  motor  could  be  significant,  but  the  gradient  is  thought 
to  be  negligible  for  a  motor  less  than  30  or  UO  cm.  in  length. 
This  effect  is  expected  to  increase  with  the  length  of  the 
motor.  In  Appendix  F  analyses  are  made  of  several  simple 
systems  with  mean  temperature  variations. 

It  is  in  this  setting  then  that  an  approximate  theoretical  des¬ 
cription  will  be  developed. 

Acoustic  Relationships 

We  wish  to  use  acoustic  theory  to  infer  the  nature  of  the  acoustic 
field  from  measurements  in  the  motor.  There  are  two  classes  of  measure¬ 
ments  which  can  be  made,  those  made  before  or  after  the  firing  and  those 
made  during  the  firing.  Those  made  during  the  firing  will  be  restricted 
to  dynamic  pressure  measurements  since  it  is  very  difficult  to  measure 
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anything  else.  However,  the  dynamic  pressure  data  will  provide  frequency 
and  amplitude  information  throughout  the  firing  for  a  particular  location 
and  phase  differences  between  locations.  Before  or  after  a  firing  it 
is  possible  to  measure  various  dimensions  in  the  motor.  Also  the 
properties  of  the  propellant  such  as  density  and  burning  rate  may  be 
obtained  in  supplemental  tests.  Now  we  shall  proceed  with  the  theo¬ 
retical  description,  recognizing  that  it  must  be  expressed  in  terms 
which  may  be  obtained  from  these  kinds  of  information. 

In  order  to  develop  an  analysis,  the  following  simplifications 
will  be  made: 

1.  The  main  body  of  the  gas  will  be  assumed  to  be  homogeneous, 
i.e.,  there  are  no  variations  in  properties  or  mean 
temperature . 

2.  The  gas  will  be  assumed  to  behave  as  a  non-reactive  ideal 
gas. 

3.  Initially  the  effects  of  mean  flow  will  be  neglected  and 
later  estimated. 

4.  The  effects  of  the  nozzle  will  be  neglected. 

5.  The  oscillations  will  be  assumed  to  be  sufficiently  small 

to  allow  application  of  the  conservation  equations  in  linear 
perturbation  form.  (See  Appendix  C  for  a  discussion  of  the 
limitations  of  this  assumption.) 

6.  The  dissipation  function  will  be  neglected,  in  accord  with 
Rayleigh  [45]. 

7*  Finally,  cylindrical  symmetry  will  be  assumed. 
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Under  these  restrictions  the  conservation  equations  may  be  written 

as 

motion- -x  direction 

°o  it +  if  ■  +  3 4 1?  v-  “  (2-s) 

continuity 

N  — * 

+  PQ  V  '  "u  =  0  (2.3) 

energy 

P0CV  If  =  kV2  T  -  p0  V  •  5  (2.4) 

Rayleigh  gives  the  solutions  to  these  equations  for  which  he  credits 
Kirchhoff  [45  ].  These  solutions,  as  modified  for  the  needs  of  the 
problem  at  hand,  are  obtained  in  Appendix  A.  The  acoustic  pressure 
and  admittance  of  the  burning  boundary  are  found  to  be 

v  =  p0  ax  Jo^kir^  eJu)t  cosh  «37x  (2-5) 

y  =  -  tanh  j7x  (2,6) 

M0 

The  origin  of  coordinates  in  these  solutions  has  been  taken  at  the 
non-burning  end  of  the  motor.  Also,  the  acoustic  velocity  has  been 
set  equal  to  zero  at  the  origin  (the  rigid,  steel  fl?nge).  In 
addition,  small  terms  have  been  neglected.  These  solutions  are  expected 
to  be  valid  everywhere  except  near  the  curved  boundaries. 
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Equation  (2.5)  will  form  the  basis  for  inferences  as  to  the 
nature  of  the  acoustic  field.  It  will  be  possible  to  determine  the 
parameters  in  this  expression  from  data.  The  determined  parameters 
may  then  be  used  to  calculate  the  acoustic  admittance  from  Equation 
(2.6).  The  radial  dependence  in  Equation  (2.5)  will  never  actually 
be  evaluated  although  it  could  be  calculated,  at  least  in  principle. 
The  argument  of  the  Bessel  function  should  be  small  enough  to  permit 
approximation  of  the  function  by  unity.  In  any  case  we  shall  give  only 
token  verbal  recognition  to  the  existence  of  a  radial  dependence. 

An  expression  for  the  propagation  constant,  7,  is  also  obtained 
in  Appendix  A.  The  expression  is  necessarily  approximate,  however, 
because  it  is  necessary  to  apply  boundary  conditions  at  the  curved 
boundaries  in  order  to  obtain  the  expression.  This  expression  is: 
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(2.7) 


Cantrell  [6]  has  obtained  a  very  similar  expression  (see  Appendix  A) 
for  a  system  in  which  the  steep  temperature  gradient  was  retained. 

The  hyperbolic  cosine  in  Equation  (2.5)  may  be  rewritten  in 
exponential  form.  In  that  case,  each  exponential  then  represents  a 
wave,  one  traveling  in  each  direction.  When  a  uniform  mean  flow  term 

V 

is  introduced,  the  effect  is  to  augment  or  suppress  the  velocity  of 
wave  propagation  depending  on  whether  the  wave  moves  with  or  contrary 

t 

to  the  mean  flow.*  The  propagation  constant  in  this  case  becomes 


approximately. 
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7  -  +  (1  +  M)  7  (2.8) 

T  — 

where  the  plus  subscript  corresponds  to  the  positive  wave,  i.e.,  the 
wave  traveling  in  the  positive  direction.  The  correction  is  small 
and  will  be  neglected. 

The  acoustic  admittance  given  in  Equation  (2.6)  corresponds  to 
the  virtual,  admittance  described  by  McClure  [ 39  J -  He  has  pointed  out 
that  there  is  a  convection  of  acoustic  energy  into  the  system  due  to 
the  mean  flow.  Equation  (2.6)  gives  a  virtual  admittance  rather  than 
an  actual  admittance  because  the  mean  flow  was  neglected.  In  Appendix 
D  appear  the  details  of  why  this  arises. 

It  is  noteworthy  that  the  expressions  for  the  acoustic  pressure 
and  admittance  (Equations  (2.5)  and  (2.6))  are  in  a  sense  more  general 
than  the  analysis  used,  if  the  parameters  are  experimentally  determined. 
The  solution  is  of  the  form 


AeJa)t  ej7X 


(2.9) 


This  is  a  valid  solution  (along  a  linear  coordinate)  for  a  great  many 
systems  with  constant  properties.  In  the  section  on  the  analysis  of 
the  solid  phase,  the  same  solution  will  be  found  to  satisfy  the  visco¬ 
elastic  equations  of  motion  for  the  solid.  Also,  as  long  as  the  viscous 
and  convective  effects  are  small,  the  equation  of  motion  (x  direction) 
may  be  approximated  as 

’0 1 ♦ !■ -  0 


or 
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(2.11) 
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where  a  harmonic  time  dependence  has  been  assumed  to  obtain  the  latter 
equation.  With  the  aid  of  Equations  (2.9)  and  (2.10)  it  is  necessary- 
only  to  apply  the  appropriate  boundary  conditions  to  obtain  Equations 
(2.5)  (except  for  the  small  radial  dependence)  and  (2.6).  In  Appendix 
F  similar  solutions  are  obtained  for  systems  with  a  mean  temperature 
variation  as  long  as  the  variations  are  not  too  large. 

Let  us  now  consider  how  the  parameters  in  Equation  (2.5)  may  be 
related  to  experimental  measurements.  The  product  of  terms  is 

simply  the  absolute  value  of  the  acoustic  pressure  at  the  origin.  The 
length  L  may  be  calculated  from  the  length  of  the  motor,  the  sample 
length,  and  the  linear  burning  rate.  We  are  left  with  the  evaluation 
of  7*  The  imaginary  part  of  7  (7  5  b-ja)  is  related  to  the  damping 
suffered  by  the  wave  and  will  always  be  very  small.  If  advantage  is 
taken  of  this  fact,  then  an  excellent  approximation  is 

—  I cos  bx |  (2.12) 

where  both  pressures  are  measured  along  the  axis  of  the  motor.  The 
value  of  aL  may  be  obtained  from  a  phase-angle  measurement  of  the 
oscillatory  pressure  at  some  point  along  the  axis  referred  tc  «he 
origin.  Thus 

tan  ♦  =  tanh  ax  tan  bx  —  ax  tan  bx  (2 .13) 

The  value  of  aL  may  also  be  approximated  from  the  oscillatory- 
decay  rate  after  the  propellant  is  consumed.  Ihis  decay  rate  may  be 
obtained  from  the  amplitude  data.  The  amplitude  of  the  wave  traveling 
in  the  positive  direction  is  attenuated  as  e 


and,  similarly,  the 


6LX 

negative  wave  is  attenuated  as  e  .  If  the  wave  is  reflected  perfectly 
from  each  end  of  the  motor  during  the  decay,  then  the  amplitude  at  a 
particular  position  will  decay  by  an  amount  e  in  one  cycle  of 
oscillation.  If  the  rate  of  decay  may  be  represented  as  e  ,  where 
is  negative, 

a 

aL  =  -  tS  (2.14) 


Thus  the  rate  of  decay  is  a  measure  of  a  also.  However,  is  some¬ 
what  dependent  upon  the  amplitude  and  in  order  to  be  appropriate,  the 
initial  value  must  be  chosen.  We  shall  use  to  denote  the  initial, 
decay.  Equation  (2.13)  is  obtained  analytically  in  Appendix  A. 

In  addition  it  is  possible  to  relate  the  growth  constant,  a, 
and  the  linear  decay  constant,  Q^,  to  a.  The  linear  decay  constant 
will  be  defined  in  the  next  section.  In  Appendix  A  the  relationship 
is  shown  to  be 


a 


(2.15) 


With  a  knowledge  of  a  and  b  the  calculation  of  the  acoustic 
admittance  may  be  essayed.  It  is  convenient  to  express  the  admittance 


in  reduced  form  by  multiplying  by  — *  ,  the  ratio  of  mean  pressure  to 


mean  velocity.  The  steady-state  continuity  equation  .equlres  that 


P0U0  “  r0Ps  (2.16) 

If  advantage  is  taken  of  the  smallness  of  a,  then  the  real  and 
imaginary  parts  of  the  reduced  admittance  may  be  written  as 


-  23  - 


Pft  Pn 

~  y'  -  -  (“)  [bL  +  tan  bL  (l  +  bL  tan  bL)  J  (2.17) 

0  rOps  ® 
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(2.18) 


Growth  Constant  Considerations 

Now  consideration  will  be  shifted  to  the  significance  of  growth 
and  decay  constants .  It  will  be  presumed  that  the  growth  and  decay 
of  the  acoustic  field  at  a  particular  position  may  be  described  by  an 
exponential  in  time,  e.g., 


at  jojt 

e  e 


(2.19) 


For  a  closed  system  the  time  average  acoustic  energy  may  be  written  as 
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wnere  <  >  indicates  a  time  average  and  the  subscripts  indicate  the 
order  of  perturbation  for  real  quantities.  If  a  «  ui,  upon  insertion 
of  equations  of  the  same  type  as  Equation  (2.19), 
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Thus,  for  a  closed  system,  the  growth  constant  is  directly  related  to 
the  rate  of  energy  input.  For  an  open  system  the  result  must  be  modified 
as  shown  by  Cantrell  and  Hart  [7J.  Their  treatment  is  outlined  In 
Appendix  B.  Their  result  is 
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This  equation  is  for  isentropic  and  irrotationai  flow  and  is  valid  to 
second  order  in  the  perturbed  quantities  and  first  order  in  Mach  Number. 

Note  the  implication  in  the  Cantrell  result  that  at  is  really  the 
sum  of  several  contributions  due  to  the  various  surfaces.  Thus 


i 


where  the  sum  is  over  all  surfaces  and  the  Ci^'s  are  positive  if  the 
contributing  surface  is  supplying  acoustic  energy.  Obviously  the  growth 
constants  cannot  be  true  constants  or  the  amplitude  would  grow  without 
bound.  The  a^'s  must  be  dependent  upon  amplitude  in  such  a  way  that 
the  losses  balance  the  gains  at  steady  state.  The  fact  that  experi¬ 
mental  data  may  be  described  by  a  constant  at  in  the  low  amplitude  range 
has  been  taken  as  an  indication  that  the  contribution  of  the  burning 
surface  is  Independent  of  amplitude  and,  similarly,  the  losses  are 
independent  of  amplitude  [27].  If  the  value  of  for  the  burning  sur¬ 
face  may  be  obtained  experimentally  and  the  Integrals  of  the  Cantrell- 
Hart  expression  evaluated,  the  admittance  of  the  buv  z  surface  may 
be  obtained.  In  order  to  do  this,  the  surface  inte£,...  *-•<  j  must  be  re¬ 
lated  to  the  admittance  of  the  burning  surface.  Coates  has  given  the 
most  frequently  used  approaches  to  estimating  from  experimental 
data  (see  Appendix  E).  Cantrell  and  Hart  have  evaluated  the  integrals 
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in  their  expression,  for  the  particular  case  of  a  T-burner,  by  neglec¬ 
ting  small  terms  and  assuming  small  deviations  from  the  zero-order 
field,  i.e.,  the  acoustic  field  for  a  closed  cylinder.  They  obtained 


a  =  - 
b 


7tuo 


(2.2k) 


where  the  relation  o>  =  which  applies  to  the  zero-order  field,  has 
been  used. 

While  considering  growth  and  decay  constants,  it  is  interesting 
to  note  that  at  steady  state,  i.e.,  when  stable  oscillations  exist. 


a  =  ab  +  aD(A)  »  0  (2.25) 

where  CCD(A)  represents  the  losses  and  is  such  a  function  of  amplitude 
that 

lim  a  (A)  =  Q_T  (2.26) 

a  U  h 

A-*0 

It  is  thus  observed  that  if  the  oscillations  decay  from  a  steady  value 
following  burn-out  (i.e.,  completion  of  burning),  the  initial  decay 
constant  will  be  numerically  equal  to  the  growth  constant  due  the 
burning  surface  before  burn-out.  That  is 


ab  -  -  aD  (2.27) 

If  the  zero-order  approximation  to  the  acoustic  field  is  assumed  and 
the  initial  decay  constant  measured,  another  estimate  of  the  acoustic 
admittance  is  thus  obtained. 

It  is  worth  noting  at  this  point  that  the  criterion  for  stability 
of  a  practical  motor  can  be  obtained  from  the  Cantrell-Hart  result. 
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It  Is  simply  required  that 

a  =  )  a.  <0  (2.28) 

i 

for  stability.  With  a  knowledge  of  the  ajyropriate  boundary  conditions 
and  sufficient  knowledge  of  the  acoustic  field  to  evaluate  the  integrals, 
the  degree  of  stability  or  instability  could  be  assessed. 


The  Response  Function 

Most  of  the  currently  available  data  on  instability  has  been 

reported  in  terms  of  the  pressure  "response  function, "  ic,  which  was 

€ 

introduced  by  Hart  and  McClure  [24].  The  response  function  may  be 
simply  related  to  the  acoustic  admittance  by  perturbing  the  mass -flux 
definition,  i.e.,  m  =  pu.  When  this  is  done  the  response  function  is 
found  to  be  numerically  equivalent  to  the  negative  reduced -virtual 
acoustic  admittance,  that  is. 
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From  this  result  Equation  (2.24)  may  be  rewritten  as 
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Thus  as  long  as  a  portion  of  the  perturbation  in  mass  flux  is  in  phase 
with  the  pressure  perturbation  and  in  the  same  direction,  the  surface 
will  have  amplifying  characteristics. 


IrA 


-  28  - 


Summary  of  Equations  to  be  Used 
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CHAPTER  III 


GAS -PHASE  EXPERIMENTAL  APPROACH 

In  Chapter  II  an  analytical  description  of  the  acoustic  behavior 
was  developed.,  In  this  chapter  we  shall  describe  how  experimental 
measurements  were  used  to  determine  the  parameters  in  that  analysis. 
The  acoustic  field  which  corresponds  to  these  pa?  ^ers  will  subse¬ 
quently  be  evaluated. 

There  are  two  classes  of  measurements  which  can  be  made  in  a 
T-burner:  those  made  during  the  firing  and  those  made  before  or  after 
the  firing.  During  a  firing  only  the  dynamic  pressure  can  be  measured 
without  great  difficulty.  Before  a  firing  such  information  as  the 
pertinent  dimensions  in  the  motor  and  physical  properties  of  the  pro¬ 
pellant  may  be  obtained.  Therefore,  for  the  problem  at  hand  one  must 
use  such  measurements  to  infer  the  behavior  of  the  acoustic  field. 

Consider  now  the  analytical  expression  for  the  acoustic  pressure, 
Equation  (2.5), 

P  *  7TP0A1J0  ^kir)  eJCOt  C0Sh  i7x  (2-5) 

It  will  never  actually  be  necessary  to  calculate  the  radial  distri¬ 
bution  but  we  must  recognize  that  It  exists.  The  task  Is  to  examine 
Equation  (2.5)  to  ascertain  the  measurements  necessary  to  determine 
the  parameters.  The  complication  of  radial  dependence,  represented  by 
the  Bessel  function  in  the  above  expression,  may  be  evaded  by  making 
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all  measurements  at  equivalent  radial  positions.  With  this  restriction, 
measurements  as  functions  of  time  of  the  pressure  amplitude  of  two 
positions,  the  phase  angle  between  them,  and  the  frequency  will  provide 
sufficient  information  to  make  the  equation  useful.  One  other  restriction 
must  be  imposed,  this  being  that  both  measurements  cannot  be  made  at 
pressure  antinodes.  The  last  restriction  is  due  to  a  kind  of  degeneracy 
in  the  equations,  i.e.,  the  phase  angle  is  precisely  *  radians  for 
all  values  of  y .  It  may  be  observed,  however,  that  the  phase  angle  need 
not  be  measured  if  the  initial  decay  rate  is  measured.  The  phase  angle 
measurement  is  more  difficult  to  make,  but  it  has  the  virtue  of  pro¬ 
viding  the  necessary  information  without  resorting  to  measurements 
made  during  transient  operations. 

The  pressure  transducers  in  use  when  this  work  was  begun  were  of 
the  uncooled  variety  so  that  it  was  necessary  to  locate  them  in  positions 
where  they  would  not  be  overheated.  This  necessity  prohibited  any 
measurements  in  the  curved  wall  of  the  motor.  However,  if  the  neces¬ 
sary  measurements  could  be  made  with  one  transducer  at  the  non-burning 
end  and  a  second  one  under  the  propellant,  then  the  transducers  would 
not  be  overheated.  Consequently,  initial  firings  were  made  to  ascer¬ 
tain  whether  or  not  a  meaningful  measurement  could  be  made  through  the 
propellant.  It  was  thought  that  when  the  propellant  was  sufficiently 
thin  the  measured  pressure  would  at  least  asymptotically  approach  the 
pressure  existing  on  the  gas  side  of  the  flame.  These  firings  indicated 
that  the  relative  amplitude  and  phase  angle  between  the  two  measure¬ 
ments  did  not  vary  significantly  with  the  thickness  of  the  propellant 
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vhen  the  propellant  thickness  did  not  exceed  1  cm.  Thereby  it  was  con¬ 
cluded  that  the  burning  surface  pressure  could  be  measured  in  this 
manner. 

We  shall  briefly  consider  the  experimented  system  at  this  point. 

A  complete  description  will  be  relegated  to  the  Appendix  (Appendix  G). 

The  T-burner  used  for  this  work  made  use  of  a  centrally  located  sonic 
nozzle.  The  burner  was  constructed  so  that  the  length  could  be  varied 
and.  thereby,  the  oscillatory  frequency.  The  acoxistic  cavity  was  3*86  cm. 
in  diameter  and  ranged  in  length  from  10  to  2 6  cm.,  nominally.  The 
mean  pressure  within  the  cavity  was  nominally  14.5  atm.,  being  con¬ 
trolled  by  the  nozzle  diameter.  The  propellant  chosen  for  the  study 
was  Utah  F  because  of  previous  experience  with  this  propellant  and  data 
available  on  it.  The  characteristics  of  this  propellant  may  be  found 
in  Appendix  H.  Propellant  samples  were  securely  bonded  into  one  end 
of  the  motor,  the  initial  sample  length  being  1.27  cm.,  nominally.  Two 
piezoelectric  pressure  transducers  were  located  axially,  one  at  each 
end.  During  a  firing  the  oscillatory  outputs  from  these  two  pressure 
transducers  were  amplified  and  recorded  on  magnetic  tape.  The  two 
data  channels  were  nominally  identical;  however,  it  proved  best  to 
record  a  reference  signal  through  each  channel  immediately  following 
the  data  from  a  firing.  These  reference  signals  permitted  compensation 
for  differences  in  circuitry  of  the  channels.  The  recorded  data  were 
later  reproduced,  filtered,  and  the  needed  information  extracted. 

The  relative  phase  angle  was  measured  with  an  electronic  phase  meter. 

The  output  of  the  phase  meter,  the  frequency  information,  tnd  all 
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amplitude  information  were  obtained  by  displaying  the  appropriate 
electrical  signal  on  an  oscilloscope  and  photographing  the  trace. 

In  Figure  5  is  shown  a  typical  example  of  the  data  reproduced 
from  the  magnetic  tape.  The  amplitude  data  were  easily  measured  from 
such  photographs.  On  the  other  hand,  the  two  phase-angle  traces 
exhibit  considerable  fluctuations.  Such  fluctuations  limited  the  ability 
to  measure  the  phase  angles  with  confidence.  The  fluctuations  occurred 
even  though  the  data  were  closely  filtered.  It  is  most  likely  that  they 
were  due  to  distortion  of  the  waveform.  The  growth  and  decay  rates  were 
also  much  easier  to  measure  than  the  phase  angles . 

The  experimental  reproducibility  was  qui'^  d.  In  Figure  6  are 
shown  pressure  traces  from  four  consecutive  firings.  These  traces 
demonstrate  the  reproducibility  attained.  The  bonding  of  the  samples 
in  place  proved  to  be  particularly  important. 
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Figure  5.  Ityplcal  phase-angle  and  amplitude  data.  Top  and 
bottom  oscillatory  envelopes  are  from  cold  and  burning  ends, 
respectively.  Top  and  bottom  irregular  traces  are  phase- 
angle  traces  from  data  and  reference,  respectively.  Distance 
between  these  traces  represents  l80°-  <J>  at  4.U°/cm. 
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Run  No.  0-223 

Oscillating  pressure  .3^  atra/cm 


Run  No.  0-225 

Oscillating  pressure  .27  atm/cra 


Run  No.  0-221* 

Oscillating  pressure  .23  atm/cm. 


Run  No.  0-226 

Oscillating  pressure  .27  atm/cm 


Figure  6.  Pressure  traces  from  four  consecutive  firings  vhich  are 
indicative  of  the  reprocucibility  achieved.  For  each  ease,  sweep 
rate  is  .5  sec/cm.  from  right  and  bottom  trace  is  total  pressure 
at  6.8  atm/cm. 


CHAPTER  IV 


GAS -PEASE  RESULTS  AND  DISCUSSION 

In  this  chapter  ve  shall  describe  the  general  character  of  the 
experimental  data.  It  will  be  pointed  out  that  the  data  are  incompatible 
with  an  acoustic  field  waich  approximates  the  zero -order  field.  The 
data  will  be  shown  to  very  strongly  indicate  that  acoustic  energy  is 
being  trapped  at  the  non-burning  end  of  the  motor,  apparently  due  to 
asymmetric  flow  conditions  in  the  vicinity  of  the  nozzle.  Finally, 
a  few  related  experimental  observations  will  be  discussed. 

In  Figure  7  are  shown  pressure  traces  from  two  typical  single- 
ended  firings;  the  first  typifies  the  pressure  behavior  in  the  2.5-  to 
4.5-kilocycle  frequency  range  and  the  second  typifies  the  behavior  in 
the  1-  to  2.5-kilocycle  range.  In  each  case  the  mean  pressure  is  seen 
to  rise  rapidly  following  ignition  and  remain  at  a  fairly  uniform 
level  until  burn-out  (termination  of  burning).  The  level  mean  pres¬ 
sure  trace  indicates  that  the  burning  area  was  uniform  throughout  the 
firing.  The  low-frequency  mean-pressure  trace  shown  in  Figure  7  is 
somewhat  irregular;  however,  this  will  be  attributed  to  effects  other 
than  mean  burning  at  the  end  of  this  chapter.  This  figure  shows  the 
kinds  of  oscillatory  envelopes  observed.  In  all  cases  the  oscillations 
began  within  100  milliseconds  after  the  ignition  transient.  The  waveform 
during  the  growth  period  was  sinusoidal  while  the  amplitude  was  low 
and  became  distorted  as  the  amplitude  became  large.  In  the  high  frequency 
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Run  No.  O-253.  Nominal  frequency  4.2  kc. 
Oscillating  pressure  .27  atm/cm. 

Total  pressure  6.8  atm/cm. 

Sweep  rate  .5  sec/cm.  from  left. 


Run  No.  234  Nominal  frequency  1.9  kc. 

Oscillating  pressure  .27  atm/cm. 

Total  pressure  6.8  atm/cm. 

Sweep  rate  .5  sec/cm.  from  left 

Figure  7.  Typical  pressure  traces  from  single-ended  firings. 
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range  the  waveform  typically  was  distorted  during  the  first  one- third 
of  the  firing  after  the  growth  transient.  The  amplitude  would  then 
drop,  frequently  to  about  half  of  the  initial  amplitude,  and  would  again 
become  sinusoidal.  In  Figure  8  ar  ■>  shown  waveforms  typical  of  these 
two  regions.  On  the  other  hand,  in  the  low  frequency  range  the  ampli¬ 
tude  would  rise  during  the  growth  transient  to  a  maximum.  The  ampli¬ 
tude  would  then  decrease  steadily  to  a  minimum  and  subsequently  grow 
again.  The  minimum  was  observed  in  all  firings  in  the  low-frequency 
range  and  always  occurred  when  about  one-fourth  of  a  1.27-cm.  sample 
remained.  This  peculiar  amplitude  behavior  will  shortly  be  associated 
with  the  energy  trap  concept  alluded  to  in  the  opening  paragraph. 

Ortce  the  data  from  a  firing  had  been  recorded,  it  was  reproduced 
in  order  to  acquire  the  desired  information.  In  Figure  9  are  shown 
the  data  from  two  typical  firings .  The  relative  amplitude  measure¬ 
ment  was  taken  from  such  photographs.  We  shall  define  "relative 
amplitude"  as  JjT|  /  jp^|  for  convenience.  Note  in  Figure  9  that  the 
'illatory  amplitudes  were  considerably  different  at  the  two  ends  of 
the  motor.  This  difference  in  amplitude  was  always  observed. 

The  relative  amplitudes  from  a  number  of  firings  are  plotted 
in  Figure  10.  The  data  shown  are  from  runs  whose  mean  pressure  traces 
indicated  that  the  burning  was  well  behaved.  Furthermore,  the  data 
were  measured  near  burn-out  and  during  stable  oscillations .  The 
effect  of  poor  mean  burning  in  the  latter  part  of  a  firing  was  to  in¬ 
crease  the  relative  amplitude,  which  remained,  however,  significantly 


less  than  unit. 
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Early  portion  of  run  (No.  0~260) 
Frequency  5-3  kc . 


Latter  portion  of  run. 
Frequency  4.6  kc. 


Figure  8.  The  change  in  waveform  frequently  observed  during  a 
firing. 
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Run  No .  0-260 . 

Nominal  frequency  4.1  kc . 

Vertical  scales  are  nominally  equal  sveep 
from  left. 


Run  No.  0-227- 

Nominal  frequency  1.4  kc. 

Vertical  scales  are  nominally  equal, sweep  from  left. 


Figure  9.  Typical  air.pl itude  data  from  single-ended  firings. 
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The  experimental  observations  that  (l)  the  relative  amplitude 
was  much  less  than  unity  and  (2)  +he  measured  phase-angle  was  slightly 
less  than  *  radians  are  of  substantial  importance.  In  particular, 
these  observations  rather  conclusively  indicate  that  the  zero-order 
acoustic -field  approximation  is  -invalid  in  these  circumstances .  In 
order  for  the  zero-order  approximation  to  hold,  the  relative  amplitude 
would  have  to  be  nearly  unity  and  the  phase  single  very  nearly  «  radians. 
These  observations  were  the  real  incentive  for  this  research  work.  If 
the  deviation  from  zero -order  was  due  to  the  combustion  process,  then 
it  would  be  of  great  importance  to  assess  the  implications  in  terms  of 
the  estimated  acoustic  admittance  for  the  burning  surface.  Consequently 
a  number  of  firings  were  made  in  order  to  characterize  the  behavior. 

In  the  past  the  zero-order  approximation  has  been  frequently  used. 
However,  there  appear  to  be  few  experimental  investigations  of  the 
validity  of  the  assumption.  Strittmater  et  al.  [55]  report  measuring 
the  frequency  before  and  after  burn-out.  On  the  basis  of  their  measure¬ 
ments  they  conclude  that  the  acoustic  field  is  nearly  zero-order. 

Horton  [28,29]  reports  measuring  the  pressure  amplitude  at  the  burning 
Interface,  the  cold  end,  and  under  a  thin  sample  of  propellant.  He 
concludes  that  the  amplitudes  are  equal  within  experimental  accuracy 
at  all  three  locations.  Thus  his  results  are  completely  contrary  to 
those  obtained  in  this  work.  However,  if  his  propellant  was  considerably 
different  from  Utah  F,  the  effects  which  are  quite  apparent  in  this 
work  may  not  have  been  evident  in  his  data.  Coates  [10,11],  in  his 
T-burner  work,  apparently  found  the  zero-order  approximation  to  be 
valid  since  he  uses  this  approximation  in  his  analysis  of  the  acoustic 
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behavior.  Therefore,  it  appears  that  this  work  is  the  first  in  which 
substantial  deviations  from  the  zero-order  field  have  been  observed, 
or,  at  any  rate,  recognized. 

In  an  attempt  to  ascertain  why  the  low  relative  amplitude  might 
occur  due  to  other  than  combustion  phenomena,  an  extensive  examination 
of  various  aspects  of  the  problem  was  undertaken.  The  details  of  this 
examination  are  given  in  Appendix  I.  The  conclusions  reached  are  as 
follows : 

1.  The  measured  amplitude  ratio  exists  within  the  motor.  In 
particular,  it  is  not  introduced  by  the  instrumentation, 
it  is  not  introduced  by  either  the  thermal  boundary  layer 
at  the  cold  end  or  the  boundary  layer  associated  with  the 
flame,  and  it  is  not  due  to  the  propellant  behind  which  the 
low  amplitude  is  measured. 

2.  Although  the  gas  may  be  cooler  at  the  non-burning  end  of  the 
motor,  the  predicted  effect  is  in  opposition  to  that  measured. 

3*  Pressure  measurements  made  in  the  curved  wadi  of  the  motor 
(propellant  side  of  the  nozzle)  are  compatible  with  the 
pressure  measured  under  the  propellant. 

4.  Measurements  made  with  the  nozzle  shifted  toward  the  pro¬ 
pellant  from  the  midpoint  of  the  motor  indicate  that  the 
pressure  node  is  located  in  the  vicinity  of  the  midpoint. 

5.  The  low  relative  amplitude  persists  through  most  of  the 
decay  transient  after  burnout.  One  can  readily  see  this  in 
Figure  11,  in  which  a  typical  decay  transient  is  shown. 


-  43  - 


Figure  11.  Oscillatory  decay  transient  from  each  end  of  the 
motor.  The  lower  trace  should  be  shifted  about  .6  cm.  to 
right  in  order  to  correspond  in  time  to  the  top  trace. 

Data  from  Hun  No.  0-260,  time  scale  is  approximately 
6.25  msec/cm. 
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The  latter  two  of  the  above  conclusions  indicate  that  perhaps 
the  analytical  description  developed  in  Chapter  II  is  Inadequate. 

First  consider  the  fourth  conclusion  given  above.  If  the 
analytical  description  were  adequate,  the  pressure  node  would  be  shifted 
toward  the  propellant  from  the  midpoint  in  order  to  satisfy  the  measured 
relative  amplitude.  More  specifically,  the  analytical  description 
indicates 


Accordingly,  a  low  relative  amplitude  requires  that  bL  <  it  and,  therefore, 

x 

the  cosine  term  will  be  zero  when  7  >  .5*  Since  the  fourth  conclusion 

h 

indicates  that  the  pressure  node  is  near  the  midpoint,  it  is  in  conflict 
with  the  description  given  by  Equation  (U.l). 

The  last  conclusion  given  above  seems  to  be  a  particularly  severe 
indictment  of  the  analytical  description  given  in  Chapter  II.  In  order 
for  the  analytical  description  to  be  adequate,  it  must  describe  the 
decay  transient  as  well  as  the  steady  state  behavior.  In  Appendix  A 
it  is  shown  that,  if  the  reflection  is  perfect  at  the  propellant  end 
after  the  propellant  is  consumed,  then  necessarily 

tan  bjjL  =0  (4.2) 

bjjL  =  n*  n  =  1,2,  3,  . . .  (4.3) 

This  also  necessitates  that 

cos  bpL  =1.0  (4.4) 

for  all  component  frequencies .  An  analytical  solution  to  the  decay 
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transient  problem  must  satisfy  the  steady-state  burning  pressure  distri¬ 
bution  as  an  initial  condition;  therefore,  the  complete  solution  would 
be  a  series  of  terms  each  of  which  would  decay  at  a  different  rate. 
However,  the  filtering  would  remove  all  except  the  fundamental  frequency. 
Thus,  in  order  to  meet  the  expected  boundary  condition  at  the  end  of  the 
motor  after  burn-out,  the  measured  relative  amplitude  should  go  to  unity 
as  the  burn-out  occurs.  Since  the  relative  amplitude  does  not  go  to 
unity,  the  analytical  description  does  not  describe  the  real  behavior. 

The  fact  that  the  unusual  acoustic  field  persists  throughout 
the  decay  transient  should  make  it  possible  to  eliminate  all  chemical 
reactions  as  possible  sources  of  the  unusual  behavior.  The  combustion 
processes  must  cease  very  soon  after  burn-out.  The  cooling  that 
accompanies  the  mean  pressure  decay  should  freeze  all  chemical  equilibria 
very  early  in  the  transient. 

On  the  basis  of  the  information  at  hand,  it  appears  tnat  a  con¬ 
siderable  portion  of  a  wave  traveling  in  the  positive  direction  (i.e., 
away  from  the  cold  end)  is  suffering  a  distributed  reflection  such  that 
the  wave  traveling  in  the  opposite  direction  is  reinforced,  only  a  part 
of  the  positive  wave  reaching  the  propellant.  The  reflection  must  be 
distributed  because  essentially  one  frequency  is  present.  The  only 
likely  source  of  such  a  reflection  is  the  flow  field  in  the  neighbor¬ 
hood  of  the  nozzle.  There  is  some  evidence  of  a  similar  reflection 
being  observed  in  a  similar  configuration  [4  S],  The  mean  flow  in 
this  region  is  sufficiently  complicated  that  no  theoretical  explanation 
can  be  presented  at  this  time;  however,  ve  shall  present  substantial 
experimental  support  for  the  postulated  behavior. 


The  analytical  description  that  has  "been  developed  should  be 
sufficiently  general  to  describe  most  phenomena  except  a  reflection 
of  the  type  we  have  described.  Consequently,  it  is  felt  that  the 
failure  of  the  analytical  pressure  distribution  actually  supports  the 
concept  of  an  asymmetric  reflection. 

In  the  light  of  this  postulated  reflection,  some  early  work 
Horton  [26]  did  may  be  significant.  He  measured  a  low  relative  ampli¬ 
tude  in  single-ended  firings .  However,  he  felt  that  his  measured 
amplitude  at  the  propellant  end  was  lew  by  a  factor  of  .45.  This  factor 
was  obtained  from  non-burning  viscoelastic  tests  in  which  a  weak  shock 
wave  was  reflected  from  one  end  of  a  propellant  sample.  In  these  tests 
the  mean  pressure  under  the  propellant  was  low  by  the  factor  of  .45. 

He  then  assumed  this  factor  was  also  appropriate  for  the  oscillating 
pressure  as  well  and  corrected  his  measured  amplitudes  accordingly. 

He  attributed  the  low  measurement  to  poor  transducer  contact.  Unfortu¬ 
nately  the  transducer  was  not  used  to  measure  a  mean  as  well  as  an 
oscillating  pressure  so  that  it  is  not  known  if  the  poor  contact  prob¬ 
lem  existed  in  his  firings.  It  seems  possible  that  Horton's  correction 
may  have  been  needless  and  that  he  was  observing  effects  similar  to 
those  we  are  describing. 

Let  us  now  return  to  consideration  of  oscillatory  pressure 
envelopes  observed.  In  Figure  7,  both  oscillatory  envelopes  were  seen 
to  substantially  drop  in  amplitude  after  the  first  one-third  of  the 
firing.  The  low-frequency  firing  was  seen  to  exhibit  a  minimum  in  the 
oscillatory  envelope.  These  observed  amplitude  reductions  are  presumed 
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to  "be  manifestations  of  the  supposed  asymmetric  reflection.  If  this 
reflection  is  due  to  the  asymmetric  flow  as  supposed,  then  a  double- 
ended  firing  should  not  exhibit  a  reduction  in  amplitude.  In  Figure  12 
is  shown  a  reproduced  oscillatory  envelope  from  a  double-ended  firing 
whose  frequency  was  somewhat  less  than  the  high-frequency  trace  shown 
in  Figure  7-  Clearly  the  envelope  does  not  show  an  amplitude  reduction. 

In  connection  with  these  amplitude  reductions,  it  was  noted  that, 
at  least  at  times,  there  was  a  shift  in  frequency  associated  with  the 
reduction.  An  example  of  this  is  shown  in  Figure  8  since  these  wave¬ 
forms  were  sampled  just  before  and  just  after  the  drop  in  amplitude. 

The  postulated  reflection,  assuming  its  development  to  be  coincident 
with  the  reduction  in  amplitude,  would  be  expected  to  affect  the  frequency. 

The  phase-angle  measurement  may  be  in  opposition  to  the  supposed 
asymmetric  reflection.  In  order  to  explain  this  we  shall  assume  that 
the  effect  of  the  behavior  in  the  vicinity  of  the  nozzle  is  to  decrease 
the  pressure  amplitude  by  an  amount  e  ,  i.e.,  for  x  >  x^ 

P*  ~  |p*g|  e’5  e<^Dt  J0  (k,r)  cosh  J7x  (U.5) 

where  x^  is  the  location  of  the  boundary  of  the  nozzle  effects .  If  the 
pressure  is  given  by  an  expression  of  this  type,  the  expression  for  the 
phase  angle  remains  unchanged,  i.e., 

tan  ♦  =  tanh  aL  tan  bL  (2.13) 

Here  it  is  noted  that  this  equation  requires  the  phase  angle  to  approach 
n  radians  as  bL  approaches  «  radians  in  order  that  aL  be  finite.  According 
to  the  reflection  concept  tne  value  of  bL  would  be  expected  to  approximate 
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Run  No.  0-264 

Nominal  frequency  3*0  kc. 

Sweep  from  left. 


Figure  12.  Oscillatory  envelope  from  a  double 
ended  firing. 
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if  radians.  In  which  case  the  value  of  aL  from  Equation  (2.1 3)  would 
be  unreasonably  large.  On  the  other  hand,  we  must  observe  that  e  5  in 
Equation  (4.5)  may  not  sufficiently  describe  the  effects  of  the  reflection. 
In  addition,  the  phase-angle  measurement  may  be  in  error.  The  measure¬ 
ment  was  made  with  great  care  and  attempts  were  made  to  compensate  for 
electronic  effects;  however,  in  spite  of  this  care,  distortion  of  the 
waveform  by  the  electronic  systems  may  have  caused  a  systematic  error 
in  the  reading.  The  fact  that  the  phase-angle  trace  from  the  reference 
signal  was  unstable  suggests  that  distortion  may  have  been  occurring. 

Real  Part  of  Propagation  Constant 

We  shall  now  examine  the  value  of  b  as  determined  from  the 
relative  amplitude  measurement  and  attempt  to  assess  its  implications 
in  terms  of  the  supposed  asymmetric  reflection. 

It  will  be  convenient  to  estimate  the  value  in  three  ways.  These 

are: 

1.  The  '  '  * ’'•e-amplitude  measurement  in  conjunction  with 

Equation  (2.12),  advantage  being  taken  of  the  sc-  lnes3  of 
a.  This  method  will  give  the  value  of  b  necessary  for  the 
analysis  of  Chapter  II  to  describe  the  end  conditions.  If 
the  low-relative  amplitude  was  due  to  the  combustion, 

then  this  value  should,  in  part,  describe  the  actual  acoustic 
field. 

2.  From  the  theoretical  propagation  constant.  Equation  (2.7), 

b  is  expected  to  approximate  .  The  postulated  effects  of 
the  nozzle  sure  not  expected  to  extend  throughout  the  motor 
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so  that  this  approximation  should  be  valid  near  the  pro¬ 
pellant  surface.  The  frequency  is  available  and  c  will  be 
estimated  by  assuming  the  mean  flow  through  the  nozzle  to 
be  isentropic  and  steady  (see  Appendix  K). 

3«  Finally  the  value  of  b  will  be  estimated  as  jr  .  This  is 
the  value  of  b  that  would  exist  if  the  acoustic  velocity 
was  zero  at  the  burning  surface. 

For  brevity  these  methods  of  estimation  will  be  called  Method  (l), 

(2),  or  (3),  respectively.  In  Figure  13>  the  three  estimates  of  b 
are  shown.  The  plotted  points  were  obtained  by  Method  (l).  Clearly 
this  method  predicts  a  much  smaller  value  of  b  than  the  other  two 
methods . 

The  low  value  of  b  predicted  from  Method  (l)  has  two  interesting 
implications,  these  being: 

1.  The  wave  propagates  faster  than  the  estimated  velocity  of 
sound.  A  wave  propagates  at  a  velocity  of  co/b  and  with 
the  indicated  b  values  this  /elocitv  **  about  thirty  per 
cent  greater  than  the  estimate  for  the  isentropic  velocity 
of  sound.  This  situation  seems  highly  unlikely  even  for 

a  chemically  reactive  system. 

2.  Since  b  is  not  nearly  equal  to  7  ,  it  is  implied  that  the 
acoustic  velocity  is  large  at  the  transpiring  boundary. 

The  veloc'J  litude  is,  in  fact,  estimated  to  be  greater 
than  the  ,  of  mean  flow.  Again,  the  situation  seems 


unlikely. 


Figure  13.  Frequency  dependence  of  estimated  b  values. 
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The  Estimated  Response  Function 

In  this  section  we  shall  show  that  the  response  function  esti¬ 
mated  with  the  developed  analytical  description  of  Chapter  II  and 
the  experimental  data,  i.e.,  relative  amplitudes  and  phase  angles  or 
decay  rates,  does  not  agree  with  an  estimate  obtained  from  growth- 
constant  data.  The  value  of  b  will  then  be  re-estimated  by  Methods 
(2)  and  (3)  and  the  corresponding  response  function  estimate  will  then 
be  shown  to  agree  reasonably  well  with  that  obtained  from  growth- 
constant  data.  This  situation  provides  additional  support  for  the 
asymmetric -reflection  concept.  The  response  function  as  estimated  from 
growth  constant  data  will  be  developed  initially. 

Very  early  in  this  research  work  growth-constant  data  were 
obtained  from  double-  and  single-ended  firings.  After  these  early 
measurements  were  made  it  was  necessary  to  purchase  new  ammonium  per¬ 
chlorate  for  making  the  propellant.  Even  though  the  new  material  was 
nominally  the  same  and  a  size  analysis  showed  nearly  equivalent  particle 
size  distributions,  the  new  propellant  v-3  fourr*  to  burn  at  a  rate 
about  forty  per  cent  higher  than  the  old  propellant  (see  Appendix  H). 

The  new  propellant  did  have  the  same  burning  rate  exponent,  however. 

No  explanation  for  the  change  in  burning  rate  has  been  suggested. 

Because  the  burning  rates  were  different  for  the  old  and  the 
new  propellants,  the  acoustic  admittance  should  be  expected  to  be 
different.  Therefore,  in  addition  to  other  data,  the  growth  constants 
were  measured  from  later  firings.  The  data  from  the  old  propellant  are 
presented  in  Figure  ik.  In  Figure  15  are  presented  the  data  on  the 
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Figure  14.  Growth  constants  for  old  Utah  F-propellant .  The 
upper  curves  are  for  double-ended  firings.  All  points  were 
obtained  at  a  mean  pressure  of  14.5  atm  (nominal). 


new  propellant;  the  dashed  lines  represent  the  older  data  and  the 
solid  lines  are  parallel  to  the  dashed  lines  but  through  the  new  data. 
Only  a  few  double-ended  firings  were  made  with  the  new  propellant 
but  these  are  adequate  for  the  comparisons  to  be  made. 

The  response  function  may  be  estimated  from  these  growth-constant 
data  in  a  manner  suggested  by  Coates  (Appendix  E).  He  suggests  that 
the  damping  be  sssumed  equivalent  for  both  double-  and  single-ended 
firings.  The  growth  constant  due  to  one  burning  surface  is  thereby 
approximated  by  the  difference  between  the  two  kinds  of  growth 
constants.  The  response  function  may  then  be  calculated  from  the 
result  given  by  Cantrell  (Equation  (2.3O)),  with  the  assumption  that 
the  acoustic  field  is  nearly  that  of  the  zero-order  field  being  made. 
This  result  is  consistent  with  that  of  Coates  (Appendix  E).  However, 
it  will  be  more  consistent  with  this  work  to  use  Equation  (2 .15). 

Thus 

Q  -  a 

aL  =  —  fr-P*  (4.6) 

or  tnrough  Coates'  suggestion 

ao  '  ai 

aL  =  —  -1-  (b.7) 

When  bL  =  *  then  Coates'  result  is  obtained,  i.e., 


The  response  function  estimates  obtained  from  application  of  this 
equation  are  shown  in  Figure  l6.  The  solid  line  and  dashed  line 
correspond  to  the  old  and  new  propellants,  respectively.  The  dashed 
line  has  been  drawn  parallel  to  the  solid  line. 
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Figure  Ip.  Growth  constants  for  new  Utah  F-propellant .  The 
upper  curves  are  for  double-ended  firings.  All  points  were 
obtained  at  a  mean  pressure  of  14.5  atm  (nominal). 
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ire  16.  Calculated  response  function  (real  part)  from  growth  constants 
Points  shown  were  calculated  from  the  growth “constant  curves. 
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It  should  be  noted  that  Coates  used  the  data  shown  in  Figure  14, 
with  a  few  of  his  own,  in  his  comparison  of  growth- constant  techniques 
(Appendix  E).  He  was  able  to  obtain  good  agreement  between  the  double- 
and  single-ended  growth- constant  estimate  of  the  response  function 
and  estimates  obtained  from  two  strictly  double-ended  techniques.  Thus 
the  response-function  estimates  shown  in  Figure  l6  should  be  reasonably 
good.  These  estimates  may  now  be  compared  with  those  obtained  from 
the  relative  amplitude  measurement  and  phase-angle  or  decay-rate 
measurements . 

The  comparison  is  facilitated  if  the  estimated  values  of  b  are 
smoothed.  The  calculated  values  of  b  obtained  by  the  three  methods 
of  the  last  section  are  plotted  on  Figure  17  as  a  function  of  reciprocal 
chamber  length.  Straight  lines  fit  the  data  adequately.  These  curves 
may  be  used  to  define  a  value  of  b  at  a  particular  chamber  length  for 
calculation  purposes. 

In  order  to  estimate  the  response  function  we  must  estimate  the 
value  of  a  in  addition  to  b.  The  value  of  a  may  be  obtained  from  the 
phase-angle  (Equation  (2.13))  or  the  initial  decay  rate  (Equation  (2.14)). 
The  phase-angle  data  may  only  be  used  with  the  relative  amplitude  data. 

The  response  function  is  related  to  the  experimental  information 
by  the  approximate  equations  (2.31)  and  (2.32) 


Re 


(~)  [bL  +  tan  bL  (l  +  bL  tan  bL)] 


S  0 


(2.31) 


^  (-)  tan  bL 


psr0 


Ira 


(2.32) 


(cm-1  ) 


Figure  17*  Inverse  length  dependence  of  b  estimated  by  Methods 
(l)  and  (2). 
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The  values  of  p^,  co,  and  L  to  be  used  are  readily  obtained  from  the 
firing  data.  The  values  of  r^  and  pg  are  obtained  from  propellant 
data. 


The  response  function  was  estimated  for  the  case  where  the  analyti¬ 
cal  description  is  assumed  to  hold,  i.e.,  the  value  of  b  being  obtained 
from  the  relative  amplitude  (Method  (l))  and  a  from  the  phase  angle 
and  decay  ate.  The  real,  part  of  response  function  as  estimated  in 
this  manner,  Equations  (2.12),  (2.3l)>  and  (2.13)  or  (2.14)  being  used, 
is  shown  in  Figure  18  with  the  estimate  obtained  from  the  double-  and 
single-ended  growth  constants.  Clearly  the  former  are  larger  by 
about  a  factor  of  two  than  the  latter. 

The  response  function  was  also  estimated  for  the  cases  where  b 
was  approximated  as  —  and  7  (Equations  (2.14)  and  (2.31)  being  used). 
Because  the  value  of  a  obtained  from  the  decay  rate  is  independent  of 
b,  a  very  simple  comparison  may  be  made.  By  means  of  the  curves  for 


& 

e 


are  in  the  ratio 


corresponding  to  b  —  ~  are  plotted.  Also  shown  is  the  curve  for 


b  shown  in  Figure  17,  the  estimates  of  Re 
2.88:  1.00:  1.11  corresponding  to  b  as  estimated  from  Methods  (l), 

(2),  and  (3),  respectively.  In  Figur-  ^9  the  calculated  estimates  of 
Re 

the  double-  and  single-ended  growth-constant  estimate.  The  agreement 
is  seen  to  be  reasonably  good  (note  that  the  vertical  scales  are  dif¬ 
ferent  on  Figures  18  and  19 ) • 

These  considerations  indicate  that  (l)  the  acoustic  behavior  in 


the  neighborhood  of  the  burning  surface  is  similar  to  a  zero-order 
field,  (2)  the  overall  deviations  from  a  zero-order  field  must  be 
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Figure  18.  Calculated  response  function  (real  part)  from 
relative-amplitude  and  phase-angle  or  decay  rate  data. 
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caused  by  effects  other  than  the  combustion.  Thus  they  appear  to  sup¬ 
port  the  concept  of  an  asymmetric  reflection. 


Imaginary  Fart  of  the  Response  Function 


The  importance  of  Im 

E 

e* 

has  not 

has  not  been  considered.  Unless  bL 


very  nearly  equals  n  radians,  the  numerical  value  of  Im 


larger  than  Re 


1± 


is  much 


In  order  to  assess  its  importance  ve  shall  esti- 


le  J 

mate  the  growth  constant  due  to  the  burning  surface  by  means  of  the 
Cantrell  equation,  i.e., 


(^•9) 


Equation  (^.9)  is  obtained  directly  from  Equation  (2.22)  by  neglecting 
a  small  term  in  each  integral.  If  bL  is  not  nearly  equal  to  n  radians, 
then  the  third  term  in  the  surface  integral  may  not  be  neglected.  This 
term  augments  the  growth  constant,  whereas,  the  contribution  due  to 
the  first  two  terms  in  the  surface  integral  is  reduced  by  the  relative 
amplitude  squared.  The  growth  constant  corresponding  to  the  various 
approximations  for  b  was  estimated  from  this  equation,  the  integrals 
being  evaluated  for  a  frictionless  acoustic  field  for  simplicity.  The 
growth  constant  corresponding  to  b  —  ®  was  found  to  be  about  five  per 
cent  greater  than  that  for  the  b  —  ^  case.  The  growth  constant  corres¬ 
ponding  to  the  b  obtained  from  Method  (l)  was  about  thirty-five  per 
cent  greater  than  that  for  the  b  —  -  case.  Thus,  as  calculated  by 
the  several  methods,  the  growth  constant  due  to  the  burning  surface  is 
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not  greatly  changed  even  though  the  response  function  changes  fairly 
drastically  from  one  method  to  another. 

Related  Observations 

In  an  effort  to  learn  more  about  the  behavior  of  the  relative 
amplitude;  this  quantity  was  plotted  through  the  firing  for  a  run  from 
each  frequency  range.  The  curves  are  shown  in  Figures  20  and  21. 

These  figures  clearly  show  that  the  relative  amplitude  was  less  than 
unity  throughout  each  firing.  However;  early  in  the  run  the  pressure 
at  the  burning  end  of  the  motor  was  measured  through  about  1.2  cm.  of 
propellant  and  this  may  be  the  reason  for  the  low  measurement  early  in 
the  run.  These  figures  also  clearly  show  that  in  the  latter  part  of 
the  run  the  measurement  is  not  substantially  dependent  upon  the  pro¬ 
pellant  length.  It  is  interesting  to  note  the  drop  in  relative  ampli¬ 
tude  at  burn-out.  This  indicates  that  at  the  propellant  end  amplitude 
decreases  in  relation  to  the  other  end  at  burn-out.  It  may  be  associated 
with  the  decay  in  the  flow  field  or  mean- temperature  field. 

It  was  also  found  that  the  mean  pressure  was  very  strongly  affected 
by  the  oscillatory  amplitude.  In  Figure  22  are  plotted  the  amplitude 
and  mean  pressure  throughout  a  firing  for  a  low-frequency  range  firing. 
The  mean  pressure  is  nearly  a  mirror  image  of  the  amplitude.  This  is 
attributed  to  the  heat  loss  from  the  gas  associated  with  the  oscil¬ 
latory  behavior.  The  rate  of  heat  transfer  at  the  walls  is  apparently 
greatly  enhanced  by  the  oscillations .  In  Appendix  K  it  is  shown  that 
the  mean  temperature  is  approximately  proportional  to  the  mean  pressure 


RELATIVE  AMPLITUDE 


RELATIVE  AMPLITUDE 


FRACTIONAL  BURNING  TIME 


Figure  21.  Measured  relative  amplitude  throughout  a  firing.  Burn 
out  occurs  when  the  fractional  time  is  1.0.  (Run  No.  0-227, 
nominal  frequency  1 . U  kc . ) . 


FRACTIONAL  BURNING  TIME 


Figure  22.  Comparison  of  the  variation  of  pressure  amplitude  and 
mean  pressure  with  burning  time.  Note  the  similarity  in  the  two 
curves.  (Run  No.  0-227,  nominal  frequency  1.4  kc.). 
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if  the  flow  through  the  nozzle  is  isentropic  and  at  steady  state. 

Therefore,  the  mean  pressure  drops  when  the  rate  of  heat  loss  is 
increased. 

The  Asymmetric  Reflection 

It  is  difficult  to  suggest  a  plausible  mechanism  for  the  postu¬ 
lated  reflection.  The  mean  flow  near  the  nozzle  must  be  complex  and, 
therefore,  even  an  approximate  analysis  of  the  wave  behavior  is  difficult. 
However,  the  flow  does  exhibit  the  asymmetry  that  apparently  is  needed. 

In  addition,  there  may  be  vortices,  or  a  vortex,  generated  by  the  mean 
flow  and  these  could  be  significant  in  terms  of  a  reflection. 

There  is  a  related,  but  small,  convection  of  acoustic  energy 
into  the  non -burning  end  of  the  motor.  This  is  observed  from  consideration 
of  the  energy  equation,  when  perturbed  to  first  order,  Integrated  over 
the  volume,  and  time  averaged  (see  Equation  (B.U)).  These  considerations 
show  that  there  are  flux  terms  comparable  to  the  energy  flux  at  the 
burning  surface.  However,  if  attention  is  focused  on  either  wave, 
rather  than  the  composite  behavior.  It  xs  found  that  there  is  an  energy 
flux  associated  with  each  wave  which  is  much  larger  than  the  convective 
effect,  the  latter  being  smaller  by  a  factor  of  the  Mach  Number. 

Therefore,  it  is  not  likely  to  be  important. 

It  should  be  possible  to  mathematically  construct  a  mechanical 
or  electrical  analog  of  the  acoustic  system  which  would  exhibit  similar 
behavior  to  that  observed.  The  construction  of  a  realistic  model  is 
hampered  by  our  poor  understanding;  however,  such  considerations  could 
be  fruitful.  The  mechanical  system  shown  in  Figure  23  was  examined  in 


Figure  23.  Mechanical  analog. 


this  regard.  The  left  side  of  the  system  shown  is  associated  with 
the  flame  and  the  right  side  with  the  cold  end.  The  relative  ampli¬ 
tude  ( analog )  obtained  from  an  analysis  of  this  system  was  found  to  be 
on  the  order  of  .5  for  a  suitable  choice  of  parameters;  however,  it 
was  necessary  to  choose  an  unduly  large  damping  term  in  order  to  achieve 
this.  The  analysis  indicates  that  other  analogs  could  be  analyzed 
which  would  exhibit  the  desired  features  and  still  retain  realistic 
parameters.  Such  considerations,  although  tenuous,  could  lead  to  a 
better  understanding  of  the  behavior  in  the  motor. 

Thus,  we  find  ourselves  unable  to  suggest  a  satisfactory  explan¬ 
ation  of  the  postulated  reflection  and  must,  therefore,  argue  only 
that  it  must  exist  in  order  tu  explain  our  experimental  observations. 

Conclusions 

In  conclusion,  we  find  that  an  asymmetric  reflection  must  exist 
near  the  nozzle  in  order  to  explain  a  number  of  observations.  This 
reflection  is  such  that  a  portion  of  a  wave  traveling  toward  the 
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propellant  surface  is  reflected  so  as  to  reinforce  the  wave  traveling 
in  the  opposite  direction.  The  Justification  for  this  conclusion  is 
found  primarily  in  the  following  observations: 

1.  The  measured  relative  amplitude  is  consistently  found  to 
be  much  less  than  unity.  A  high  level  of  confidence  is 
placed  in  this  measurement  as  representative  of  the  ratio 
of  amplitudes  existing  at  the  two  ends  of  the  gas  cavity. 

The  low  value  is  thought  to  be  due  to  the  asymmetric  reflection 
because  (l)  calculated  values  for  Rej^iij  agree  with  the  best 
available  estimate  when  the  reflection  is  assumed;  (2)  other¬ 
wise  the  calculated  velocity  of  wave  propagation  is  unduly 
large  and  the  calculated  velocity  amplitude  at  the  propellant 
surface  is  unduly  large. 

2.  Analysis  of  possible  mean  temperature  effects  shows  that 
these  effects  are  in  opposition  to  the  low  relative  amplitude 
observed. 

3.  The  low  relative  amplitude  persists  through  most  of  the  decay 
transient.  This  is  expected  to  be  possible  only  if  there  is 
an  asymmetric  reflection. 

4.  The  low  relative  amplitude  observed  cannot  be  due  to  high 
attenuation  effects  because  of  the  low  aQ  and  persistence  of 
the  low  relative  amplitude  in  the  decay  transient. 

5 .  Measurements  made  in  the  curved  wall  and  near  the  propellant 
indicate  that  the  asymmetric  reflection  occurs  but  not  near 


the  propellant  surface. 


\ 
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6.  The  flow  pattern  near  the  nozzle  is  the  only  apparent  source 
of  the  necessary  asymmetry.  Comparison  of  the  oscillatory 
envelopes  from  double-  and  single-ended  firings  supports 
this  concept. 


CHAPTER  V 


LCW -FREQUENCY  RESPONSE  FUNCTION 

In  this  chapter  we  shall  discuss  the  measurement  of  the  acoustic 
admittance  under  conditions  where  the  gas -cavity-burning-surface  system 
is  actually  stable.  The  manner  in  which  these  measurements  were  made 
is  a  rather  simple  extension  of  standard  T-burner  technique,  but  the 
approach  has  not  been  previously  reported.  The  work  was  supported  by 
Stanford  Research  Institute  (SRl).* 

SRI  supplied  two  propellant  formulations  which  were  being  used 
in  their  work  and  requested  that  the  acoustic  admittance  of  these  pro¬ 
pellants  be  measured  at  35  atm.  and  in  a  motor  102  cm.  in  length. 

The  initial  experimental  approach  was  to  measure  double-  and 
single-ended  growth  constants.  Toward  that  goal  the  T-burner  was 
lengthened  to  102  cm.  and  several  firings  made.  These  firings  shewed 
that  it  would  be  impossible  to  attain  the  required  mean  pressure  during 
the  growth  transient  with  a  nozzled  burner.  The  burner  was,  therefore, 
adapted  for  surge-tank  operation  so  that  the  required  mean  pressure 
could  be  obtained  by  pre-pressurization.  When  the  burner  was  operated 
in  this  manner  it  was  found  that  the  pressure  pulse  from  the  igniter 
excited  the  fundamental  mode  of  oscillation.  However,  these  osci? lations 

# 

Mr.  Lionel  A.  Dickinson  of  SRI  very  kindly  gave  permission  for 
the  author  to  use  the  information  reported  here. 
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decayed  out  during  the  burning,  the  energy  supplied  by  the  flame  being 
insufficient  to  maintain  the  oscillations.  Therefore,  the  motor  was 
actually  a  stable  system. 

The  admittance  calculation  from  growth  constants  does  not  require 
that  the  growth  constants  be  positive.  Thus  double-  and  single-ended 
growth  constants,  even  though  they  are  negative,  may  still  be  used. 

In  addition,  the  fact  that  the  igniter  excited  the  gas  cavity  made  it 
possible  to  measure  the  damping  in  the  motor  by  firing  an  igniter  in 
an  unloaded  motor,  i.e.,  with  no  propellant  sample. 

Experimental  Work 

The  surge  tank  used  for  this  work  was  too  small  to  prevent  a 
fairly  significant  rise  in  mean  pressure;  however,  meaningful  results 
were  still  obtained.  The  tank  was  pre-pressurized  to  32  atm.  with 
nitrogen  and  the  subsequent  ignition  raised  the  pressure  to  35  am*  in 
the  period  when  the  negative  growth  constant  was  measured.  The  igniter 
was  made  up  of  F-propellant  sawdust  and  Pyrofuse  wire,  the  latter 
being  electrically  ignited.  Several  firings  were  made  with  each  of 
the  propellants  supplied  and  with  the  igniter  alone.  In  Figure  Sh¬ 
are  shown  pressure  traces  from  a  typical  firing. 

Calculated  Results 

The  data  obtained  were  sufficient  to  permit  estimation  of  the 
response  function  for  the  propellants.  The  data  are  summarized  in 
Table  1  and  the  calculated  values  in  Table  2.  The  SRI -111  propellant 
burned  rather  fast  so  that  when  it  was  used  in  a  double-ended  firing 
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Figure  24.  Typical  pressure  traces  from  102-cm.  T-burner 
firing.  Top  trace,  oscillatory  pressure  followed  by 
reference,  .68  atm/cm.  Lower  trace,  total  pressure, 

6.8  atm/cm.,  from  3^*8  atm.  initial  level.  Sweep  rate 
was  .5  sec/cm.,  from  left. 
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TABLE  1 

SUMMARIZED  DATA 


Number  of 
Firings 

Frequency 

cps 

°Ll 

Sec.  x 

Igniter  only 

3 

180 

-19* **8  +  3 

109  Propellant 
Single-ended 

3 

220 

-11.8  +2.4 

Double-ended 

3 

240 

-  3-8  +  6 

111  Propellant 
Single-ended 

4 

205 

-  8.4+  .8 

Double-ended 

1 

— 

0 

F  Propellant 
Single- ended 

1 

240 

-  8.9 

TABLE  2 

CALCULATED  RESULTS 


Propellant 

109 

111 

F 

ai 

GDL 

8.0  +  2.7 

11.4  +  l.l 

10.9 

£ 
r— l|OJ 

-  V 

8.0  +  .5 

9-9* 

— 

a2 

- 

al 

8.3  +  3-0 

8.4* 

— 

Re 

,e 

** 

6.5 

l6 

7.5 

*  C*2  =  0  used 

**  (a-j^  -  CtDL)  used 
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the  mean  pressure  rose  rapidly.  This  prevented  obtaining  a  reliable 
value  of  the  double-ended  grovth  constant  for  this  propellant.  One 
double-ended  firing  was  made  and  the  grovth  constant  measured  in  the 
early  (and  probably  nonlinear)  portion  of  the  decay,  the  value  being 
-3 .6  sec.”\  A  better  approximation  is  thought  to  be  zero. 

The  response  function  may  be  calculated  from  Equations  (2.15 ), 
(2.3l),  and  the  zero-order  field  approximation.  Thus 

Re 


vhere 

“b  ■  C2  -  “l 

the  density,  pg,  was  unknown  for  the  SRI  propellant  and  was  approxi¬ 
mated  as  I.65  gm/cm^. 

Significance  of  Work 

This  work  is  particularly  significant  because  it  offers  a  new 
means  of  estimating  the  acoustic  properties  of  relatively  stable  pro¬ 
pellants.  One  of  the  greatest  disadvantages  of  the  T-burner  method 
has  been  the  inability  to  study  propellants  too  stable  to  drive  the 
burner.  However,  from  this  work,  it  appears  that  stable  propellants 
can  be  studied  by  suitably  pressure  pulsing  a  T-burner.  It  is  quite 
likely  that  suitable  pulses  could  be  generated  to  excite  motors  much 
shorter  than  100  cm.  so  that  the  method  could  be  used  at  higher 
frequencies.  Horton  observed  oscillations  being  generated  by  his 
igniter  at  frequencies  less  than  600  cps.  in  a  T-burner  (29] .  The 


Vo 
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DL 
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accuracy  of  this  method  is  not  great  because  of  the  necessity  to 
subtract  two  kinds  of  growth  constants.  In  this  regard  ^  (a^  + 
gives  the  highest  accuracy.  However,  the  simplicity  of  the  method 
greatly  compensates  for  the  lack  of  accuracy.  One  may  compare  the 
simplicity  of  this  appraoch  to  the  work  of  Poner  [15]*  In  order  to 
apply  the  T-burner  to  relatively  stable  propellants,  he  has  found  it 
necessary  to  construct  a  very  complex  electronic  and  acoustic  system. 
The  method  we  have  used  at  least  will  permit  the  ranking  of  the 
instability  characteristics  of  various  propellants  which  are  too 
stable  to  maintain  oscillations  in  a  T-burner. 

The  magnitude  of  the  calculated  response  function  is  also  of 
significance.  These  values  are  as  large  or  larger  than  any  previously 
reported . 


*P 


CHAPTER  VI 


SOLID-PHASE  THEORY 

As  noted  in  Chapter  I,  the  solid  phase  work  developed  as  an 
extension  of  the  work  by  Coates  [ll]  and  Ryan  and  Coates  [4-9]. 

Coates  constructed  a  T-burner  in  which  a  long  grain  of  propellant 
was  advanced  as  it  burned  so  that  the  burning  interface  remained  at 
a  nearly  fixed  position.  A  sketch  of  the  burner  is  shown  in  Figure 
25.  The  propellant  fit  closely  into  the  motor,  the  curved  boundary 
being  lubricated  with  grease,  and  was  bonded  to  the  face  of  the 
piston.  Two  pressure  transducers  were  used;  one  measured  the  pres¬ 
sure  at  the  face  of  the  piston  and  the  second  measured  the  pressure 
at  the  non-burning  end.  Both  transducers  were  located  along  the  axis 
of  the  cylinder.  Coates  found  tl.it  the  solid  participated  strongly 
in  the  oscillatory  behavior  At  times  he  found  the  solid  dissipated 
all  of  the  acoustic  energy  '  ipplie^  by  the  combustion  and  the  oscil¬ 
lations  ceased.  The  cessation  apparently  was  associated  with  resonance 
phenomena  in  the  solid  He  was,  however,  encumbered  in  the  analysis 
of  his  data  by  the  inability  to  accurately  measur*  phase  angles. 

More  recently  a  theoretical  description  of  the  solic  phase  was 
developed  by  Ryan  and  Coates  [49].  Coates'  data,  which  had  been 
recorded  on  magnetic  tape,  were  re-evaluated  by  the  author  usi.ng  the 
electronic  phase  meter.  A  reprint  of  the  publication  of  this  work  may 
be  found  in  Appendix  L.  These  data,  in  terms  of  the  analysis,  suggested 
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ier  configuration  for  solid-phase  work 


"78 


a  severe  degradation  in  visco-elastic  properties.  In  order  to  study 
this  apparent  degradation  the  author  made  a  number  of  firings  using 
Coates’  burner. 

The  author’s  experimental  data,  though  similar  to  Coates’  data, 
led  to  some  apprehension  about  the  validity  of  the  one-dimensional 
analysis  used  by  Ryan.  Therefore,  a  three-dimensional  (cylindrically 
symmetric)  analysis  was  developed.  This  analysis  showed  that  the 
apparent  change  in  properties  could  equally  well  be  attributed  to  a 
change  in  conditions  at  the  curved  boundaries  of  the  solid.  The 
latter  possibility  is  thought  to  be  the  more  likely  case.  The  low 
relative  amplitude  found  in  the  gas-phase  work  must  also  be  considered. 

We  shall  proceed  by  developing  a  three-dimensional  analysis  of 
the  wave  behavior  in  the  solid  and  then  attempt  to  els s ess  the  impli¬ 
cations  of  the  experimental  data. 

Wave  Behavior  in  the  Solid 

As  in  the  gas-phase  work,  it  will  be  necessary  to  infer  the  wave 
behavior  from  dynamic  pressure  measurements  made  in  the  motor.  Also 
the  dimensions  and  density  of  the  propellant  will  be  known.  Unfortunately 
the  visco-elastic  moduli  are  not  known  for  this  propellant  and  these 
will  have  to  be  estimated. 

In  order  to  develop  a  theoretical  description  it  will  be  neces¬ 
sary  to  make  the  following  assumptions. 

1.  The  propellant  may  be  treated  as  homogeneous  and  isotropic. 

The  assumption  of  homogeneity  is  Justified  because  the  scale 
of  the  heterogeneity  is  on  the  order  of  100  microns.  The 
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composition  of  the  propellant  is  given  in  Appendix  H. 

2.  Cylindrical  symmetry  will  be  assumed. 

3.  The  frequency  will  be  assumed  to  be  constant. 

4.  Linear  viscoelastic  behavior  will  be  assumed,  i.e.,  the 
viscoelastic  constants  are  only  functions  of  the  frequency. 

The  viscoelastic  constants  will  be  treated  as  complex  in 
accord  with  the  usual  technique  [.13>59]* 

5.  The  wave  behavior  is  at  steady  state  even  though  the  length 
is  changing  slowly. 

Under  these  restrictions  the  wave  equations  for  the  solid  are  [35] 


a2s 


Pa?£m  (x  +  2tl)  5?  +  2,1  ~sl 


where 


a2e. 


p  i?* =  (x  +  *u)  $£ '  * 


1  d(rSr)  3£x 
r  ST"  +  ~5x 
as  as 

0  _ r  _ x 

2  ®q  =  “3£  ‘  37“ 


(6.1) 


(6.2) 


(6.3) 

(6.4) 


The  solutions  to  these  equations  are  given  by  Love  and  are  known  as 
the  Poc hammer -Chree  solutions  [l4,  35)*  Assume  solutions  of  the  form 

s  .  0  ej(7x  +mt)  (6.5) 

r 

£  .  w  ej(7x  *  mt)  (6.6) 

where  U  and  W  are  functions  of  r  only.  The  7  used  here  is  again  a 
propagation  constant.  Through  the  use  of  Equations  (6.5)  and  (6.6), 
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the  wave  equations  may  be  separated  into  two  equations  of  the  Bessel 
type  for  A  and  Equations  (6.3)  and  (6.4)  may  then  be  used  to 

obtain 


where 


U  “  A  Jo^hr^  +  B7J1 

W  =  AJ7Jq  (hr)  +  ^  trJ2.  (kr)J 


pcs 

X+2p 


(6.7) 

(6.8) 

(6.9) 

(6.10) 


If  7  is  replaced  by  a  negative  7 ,  the  solutions  are  the  same 
except  that  the  expression  for  W  is  replaced  by  the  negative  of  that 
given  in  Equation  (6.8).  For  this  analysis  we  shall  choose  the  origin 
at  the  face  of  the  piston.  The  piston  should  be  sufficiently  rigid 
so  that  £  will  be  zero  at  that  surface.  Therefore 


£x  =  W  e^ay*'  sinh  j7x 

t  -  U  cosh  J7x 

r 


(6.11) 

(6.12) 


The  pressure  measured  at  the  origin  is  equal  to  the  negative  of 
the  normal  stress  at  this  point.  Similarly,  the  normal  stress  at  the 
burning  end  of  the  propellant  may  be  estimated  from  the  pressure 
measured  at  the  non-burning  end  of  the  motor.  The  normal  stress 


is  composed  of 
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from  which  it  is  found  that 


(6.13) 

(6.14) 


p'f  =  |ps|  eJcut  cosh  irl  (6.15) 

Equation  (6.15)  forms  the  basis  of  the  analysis  as  it  provides 
the  means  for  calculating  7  from  the  experimental  data.  Ryan  obtained 
an  equivalent  expression  in  his  analysis.  When  his  expression  was 
applied  to  the  experimental  data.,  7  was  found  to  vary  with  time.  It 
was  this  variation  that  led  Ryan  and  Coates  to  the  conclusion  that 
the  propellant  properties  were  changing. 

Theoretically  7  is  obtained  by  applying  the  boundary  conditions 
at  the  curved  wall.  The  actual  stress  condition  at  this  boundary  is 
not  known.  Since  the  clearance  is  small  (—  .01  cm.)  and  the  gap  is 
filled  with  grease,  one  might  expect  a  limited  radial  motion.  However, 
when  a  radial  displacement  is  estimated  from  the  data  and  a  Poisson's 
ratio  of  .5,  the  magnitude  is  found  to  be  small  compared  to  the  gap. 
Thus  wc  shall  examine  three  limiting  conditions:  a  free  boundary,  a 
freely-supported  (no  radial  motion  but  free  axial  motion)  boundary, 
and  a  fixed  boundary. 

A  free  boundary  requires  the  stress  to  be  zero  at  the  wall. 

The  stress  components  are  given  by  Equation  (6.14)  and 


T 

rr 


X* 


(6.16) 
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The  free -boundary  requirement  leads  to 


(£—  -  272)  J±  (kr) 


‘  *  J0<hr)  oA 

2^  — —  x+2(i  JQ(hr) 


-  27 


3j0(hr)  r  dj-^kr) 


(6.17) 


r=r 


w 


Truncated  series  expressions  may  be  used  to  approximate  the  Bessel 
2 

functions.  If  r  is  neglected,  then 


2  pen 


(6.18) 


where 


E  « 


n(3X+2p) 


(6.19) 


A  freely  supported  boundary  requires  that  i r  and  be  n^l  at 
the  wall.  In  this  case  it  is  found  that 


djQ(hr)| 
~5r 


=  0 


(6.20) 


r»r 


*  w 


The  lowest  root  gives 


r2  =  -£SL 
r  X+2p 


(6.21) 


All  larger  roots  correspond  to  highly  damped  motions  and  should  not 
occur.  Equation  (6.21)  is  precisely  the  result  found  if  the  wave 
equations  are  restricted  to  one-dimensional  behavior. 

Finally  a  fixed  boundary  requires  that  6^  and  fcr  be  nil  at  the 
wall.  Here  it  is  found  that 
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[rJ^(kr)  ]  -  J-^kr)  JQ(hr)  yc 


(6.22) 


Again,  series  expansions  may  be  used  to  approximate  the  Bessel  functions. 

2 

If  r  is  neglected,  the  frequency  becomes  zero.  Therefore,  no  simple 
expression  is  obtained  in  this  case. 

Clearly  the  propagation  constant,  y ,  can  vary  if  the  boundary 
condition  changes  during  a  firing.  Later  the  magnitude  of  such  changes 
will  be  estimated.  In  a  one-dimensional  analysis  the  propagation  con¬ 
stant  is  only  dependent  upon  the  viscoelastic  properties.  Consequently 
the  analysis  requires  the  properties  to  change  if  the  propagation  con¬ 
stant  changes. 


Summary  of  Equations  to  be  Used 

The  following  equations  will  be  important  in  later  discussions: 


pf  =  lps| 

e*^  cosh  iyl 

(6.15) 

cvj  i 

3I» 

n 

CVI 

(free  boundary) 

(6.18) 

y2  „  £2>L 
'  X+2  u 

(freely  supported  boundary) 

(6.21) 
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CHAPTER  VII 


SOLID-PHASE  EXPERIMENTAL  APPROACH 

The  experimental  apparatus  used  in  this  work  was  that  constructed 
by  Coates.  His  motor  was  briefly  described  in  Chapter  VI.  From  that 
discussion  the  reader  will  recall  (see  Figure  25)  that  pressure  measure¬ 
ments  were  made  at  both  ends  on  the  motor,  i.e.,  under  the  propellant 
and  at  the  cold  end.  The  cold-end  measurement  is  used  to  infer  the 
pressure  at  the  surface  of  the  propellant.  This  inferred  value  and 
the  second  pressure  measurement  may  then  be  used  to  calculate  the 
propagation  constant,  7,  from  Equation  (6.15).  The  calculation  was 
facilitated  by  a  set  of  parametric  curves  which  were  developed  from 
Equation  (6.15).  These  curves  allowed  a  direct  calculation  of  the 
real  and  imaginary  parts  of  7  from  the  phase  angle  and  amplitude 
ratio  between  the  two  pressures.  The  curves  may  be  found  in  Appendix  M. 
Experimentally,  therefore,  It  was  necessary  to  measure  the  phase  angle 
and  amplitude  ratio  from  the  pressure  data  for  suitable  firings. 

The  experimental  system  differed  from  that  used  in  the  gas -phase 
study  only  by  the  inclusion  of  a  system  to  position  the  propellant. 

The  latter  system  was  comprised  of  a  hydraulically-driven  piston 
and  a  driving-rate  control.  The  sensing  element  for  the  control  was 
a  thermocouple  mounted  in  the  wall  of  the  motor  adjacent  to  the  burning 
interface.  The  control  varied  the  driving  rate  so  as  to  maintain  a 
fixed  thermocouple  temperature.  The  propellant  was  bonded  to  the  face 
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of  the  piston  and  against  the  pressure  transducer.  When  assembled, 
the  curved  wall  of  the  propellant  was  coated  with  grease  which  pro¬ 
vided  lubrication,  prevented  burning  along  the  side,  and  filled  the 
thin  gap  between  the  propellant  and  the  motor.  A  complete  description 
of  the  apparatus  and  procedures  may  be  found  in  Appendix  G. 

A  number  of  firings  were  made  in  this  mode  of  operation.  Many 
of  these  did  not  provide  usable  data  due  to  poor  position  control,  poor 
mean  pressures,  or  electronic  difficulties.  The  control  system  was 
found  to  be  adequate  only  for  small  corrections  in  the  driving  rate. 

It  was,  therefore,  necessary  to  adjust  the  average  driving  rate  until 
It  was  very  close  to  the  burning  rate.  Poor  mean  pressure  occurred 
when  the  grease  did  not  prevent  burning  along  the  walls  of  the  motor. 
Electronic  difficulties  frequently  arose  because  of  the  complexity  of 
the  system.  However,  in  spite  of  these  problems,  usable  data  were 
acquired. 


CHAPTER  VIII 


SOLID-PHASE  RESULTS  AND  DISCUSSION 

In  this  chapter  the  general  nature  of  the  experimental  data 
will  be  discussed.  An  apparent  variation  in  the  propagation  con¬ 
stant  with  burning  tim$  reported  by  Ryan  and  Coates  [49]  and  con¬ 
firmed  in  this  work,  will  be  considered.  Although  the  asymmetric 
reflection  may  be  responsible,  its  effects  cannot  be  assessed.  In 
what  follows,  it  will  tentatively  be  assumed  that  the  apparent  vari¬ 
ation  occurs  in  the  solid  phase.  The  importance  of  the  boundary  con¬ 
ditions  in  determining  the  propagation  constant  will  be  examined. 
Finally  it  will  be  concluded  that  though  little  can  be  said  with 
confidence,  the  apparant  variation  is  more  likely  due  to  wall  effects 
(or  the  asymmetric  reflection)  than  to  degradation  of  the  propellant. 

Character  of  the  Data  Obtained 

In  Figure  2 6  are  shown  the  recorded  data  from  two  firings.  These 
are  typical  of  the  kinds  of  data  obtained.  From  such  photographs  the 
ratio  of  amplitudes  and  phase  angle  were  measured  throughout  the  firing. 
Note  the  irregularities  in  the  phase-angle  traces.  Ihese  irregularities 
can  frequently  be  associated  with  changes  in  amplitude  so  that  it  is 
unlikely  that  they  are  caused  by  distortion  of  the  waveform.  Observe 
also  that  neither  does  the  amplitude  ratio  approach  one  at  burn-out 
nor  does  the  phase  angle  reach  180  degrees.  These  last  two  observations 
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were  frequently  made  and  they  motivated  the  gas-phase  work  described 
earlier. 

The  failure  of  the  phase  angle  to  approach  l80  degrees  is  probably 
due  to  the  instrumentation.  Unfortunately  no  reference  signal  was 
recorded  during  this  work.  Instead  the  phase  shift  due  to  the  instru¬ 
mentation  was  measured  and  found  to  be  less  than  five  degrees  for 
frequencies  less  than  four  kilocycles.  An  error  this  great  at  that 
time  was  considered  acceptable.  An  examination  of  Figure  26  shows 
that  the  final  phase  angles  are  much  less  them  175  degrees.  Phase 
angles  on  the  order  of  179  degrees  were  measured  in  the  gas -phase  work 
so  that  this  phase  angle  should  at  least  approach  175  degrees.  We 
offer  no  explanation  for  the  difference  but  recognize  that  it  must  be 
considered  as  a  possible  error. 

It  may  also  be  observed  In  Figure  26  that  the  relative  amplitudes 
at  burn-out  (about  .65)  are  not  as  low  as  the  .45  observed  in  the  gas- 
phase  work.  The  instrumentation  was  adjusted  so  that  each  signal 
would  be  amplified  the  same  amount.  Therefore,  the  Instrumentation  is 
not  expected  to  have  caused  the  entire  discrepancy.  It  appears,  then, 
that  we  must  correct  the  ratio  of  amplitudes  to  be  used  in  the 
calculations,  but  the  amount  of  correction  is  uncertain. 

We  may  observe  that  perhaps  we  should  not  expect  the  relative 
amplitude  to  reach  .45  at  burn-out.  It  Is  not  unlikely  that  the  amount 
of  asymmetric  reflection,  which  undoubtedly  occurs,  may  vary  during  an 
extended  firing.  Conceivably  the  observed  difference  in  relative  ampli¬ 
tudes  is  an  indication  of  such  a  variation.  Thus,  not  only  is  the 
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amplitude  correction  uncertain,  but  also  the  application  of  a  constant 
correction  is  questionable. 

The  cessation  of  oscillations,  which  Coates  [ll]  observed,  never 
occurred  in  this  work. 

Calculated  Propagation  Constant 

The  propagation  constant  was  calculated  from  the  data  by  three 
methods,  namely, 

1.  by  neglecting  the  difference  in  amplitude  at  the  two  ends  of 
the  gas  phase  and  assuming  a  phase  shift  in  the  gas  of  180 
degrees, 

2.  by  assuming  the  relative  amplitude  in  the  gas  phase  to  be 
.1+5  and  the  phase  shift  to  be  180  degrees, 

3.  by  assuming  corrections  for  the  amplitude  and  phase  angle 
such  that  the  amplitude  ratio  approached  unity  and  the 
phane  angle  approached  180  degrees  at  burn-out. 

The  calculated  propagation  constant  for  each  of  these  methods  varied 
with  the  propellant  length,  the  last  method  showing  least  variation. 

On  Figure  2J  are  plotted  the  calculated  values  of  al  and  bl  as  calcu¬ 
lated  by  Method  (3)  above.  These  are  typical  of  the  kinds  of  results 
obtained.  The  data  and  calculated  values  for  a  number  of  firings  may 
be  found  in  Appendix  M.  If  the  propagation  constant  were  independent 
of  length,  the  pointB  would  fall  near  a  straight  line  through  the 
origin.  Clearly  these  calculated  points  are  not  well  described  by  a 
non-length  dependent  propagation  constant.  Here  we  must  recognize 
the  possibility  of  an  apparent  variation  arising  from  changes  in  the 
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ratio  of  amplitudes  between  the  two  ends  of  the  gas  column,  i.e., 
changes  in  the  amount  of  asymmetric  reflection.  In  the  gas-phase  work 
the  reflection  was  found  to  develop  within  .3  second  from  ignition; 
therefore,  in  these  firings  it  is  expected  to  have  occurred  essentially 
throughout.  However,  the  separation  of  these  gas -phase  effects  from  the 
solid-phase  effects  is  impossible  without  a  better  understanding  of 
the  asymmetric  reflection.  Therefore,  lacking  anything  better,  we 
shall  tentatively  assume  that  the  propagation  constant  for  the  pro¬ 
pellant  is  changing  and  attempt  to  learn  what  we  can  from  the  data 
obtained . 

Importance  of  Boundary  Conditions 

In  the  theoretical  discussion  of  the  solid  phase  it  was  shown 
that  the  propagation  constant  is  determined  (l)  by  the  viscoelastic 
properties  and  (2)  by  the  conditions  at  the  curved  boundaries.  In  the 
analysis  of  Ryan  and  Coates  [49]  the  second  possibility  was  eliminated 
by  the  use  of  a  one-dimensional  analysis.  They  were,  therefore,  led  to 
the  conclusion  that  the  properties  must  change.  However,  we  shall 
shortly  demonstrate  that  changes  in  the  boundary  condition  are  of  great 
importance. 

It  may  be  noted  that  if  Ryan  and  Coates  [49]  had  used  Method  (3) 
rather  them  Method  (l)  in  their  calculations,  they  would  have  found 
the  variation  in  the  propagation  constant  to  be  significantly  reduced. 

The  fact  that  the  solid  boundary  condition  may  strongly  affect  the 
acoustic  behavior  hem  been  noted  by  Angelus  (l),  Horton  (26),  and 
Coates  (ll)  among  others.  Horton  ran  viscoelastic  tests  by  reflecting 
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a  weak  shock  •  ave  from  the  end  of  a  cylindrical  grain  of  propellant. 

He  found  that  the  propellant  behavior  was  strongly  affected  by  the 
boundary  conditions.  In  a  similar  experiment,  Coates  found  that  his 
propellant  would  not  oscillate  and  attributed  this  to  wall  losses. 

From  Horton’s  data  we  may  estimate  the  effect  of  changing  the 
wall  condition  from  a  free  to  a  fixed  boundary.  We  estimate  b  for 
the  free  boundary  to  be  about  three  times  that  for  a  fixed  boundary. 
Horton  reports  an  attenuation  term  only  for  the  fixed-boundary  case 
so  that  the  difference  in  a  may  not  be  estimated.  For  the  fixed 
boundary  a/b  =  .2  is  estimated. 

If  the  viscoelastic  constants  were  known  for  Utah  F  propellant, 
Equations  (6.18)  and  (6.21 )  could  be  used  to  estimate  the  magnitude 
of  possible  wall  effects.  We  shall,  therefore,  estimate  these  values. 

The  value  of  n  will  be  estimated  from  shear  data  [ h-6  ]  for  a  propellant 
similar  to  Utah  F.  The  value  of  X  will  be  estimated  from  compressibility 
data  for  a  propellant  somewhat  less  similar  to  Utah  F.  The  compressi¬ 
bility  at  zero  frequency  will  be  used  since  we  have  found  no  other 
data  in  the  literature;  but  in  the  frequency  range  of  intereso  the 
compressibility  should  differ  by  no  more  than  a  factor  of  two  [lU, 

By  means  of  this  approach,  we  estimate  that 

p  =»  2  (l+jKlO9)  dyne/cm2  (8 . 1 ) 

X  a  10  (l0^)  dyne/cm2  (8.2) 

For  comparison,  we  estimate  from  Horton's  data  that  n  =  5  (lO'j  dyne/cm' 
for  his  propellant,  Equation  (6.18)  and  u  «  X  being  used.  By  means  of 
Equations  (6.l8)  and  (6.2l)  it  is  found  that 
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7  =  .8  -  j . 3  free  boundary  (8.3) 

7  =>  .09  -  j.0009  freely  supported  boundary  (8.4) 

From  these  considerations  it  is  clear  that  vail  effects  can  be  large 
enough  to  account  for  the  observed  changes  in  the  propagation  constant. 

It  seems  likely  that  the  thin  grease  film  at  the  wall  of  the 
propellant  could  give  rise  to  additional  wall  effects  in  the  form  of 
viscous  losses.  In  this  case  the  shear  stress  at  the  wall  should  be 
proportional  to  the  rate  of  strain  at  the  wall.  An  examination  of 
the  theory  (Equations  (6.1l),  (6.12),  and  (6.14))  shows  that  in  this 
case  the  wall  stress  varies  with  x.  Therefore,  although  the  solutions 
of  Chapter  VI  are  invalid  for  this  case,  it  appears  that  such  wall 
effects  could  give  rise  to  an  apparently  variable  propagation  constant. 

Indications  of  Wall  Effects 

It  appears  that  there  is  no  way  to  separate  the  effects  of  the 
boundaries  and  the  effects  of  the  viscoelastic  properties;  however, 
there  are  a  few  experimental  observations  which  seem  to  favor  the 
former . 

1.  The  data  show  some  rapid  fluctuations  in  phase  and  amplitude 
which  are  probably  due  to  fluctuations  in  the  boundary 
conditions . 

2.  The  calculated  values  of  a  are  larger  than  expected  from 
those  estimated  for  the  various  boundary  conditions.  This 
may  well  be  an  indication  of  strong  viscous  losses  at  the 
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wall.  On  the  other  hand,  one  successful  firing  with  an 
increased  clearance  at  the  wall  did  not  produce  lower  a 
values  (Run  O-189,  also  see  Figure  27. ) 

3.  The  calculated  propagation  constant  frequently  changes  in 
the  direction  to  indicate  increased  strength  rather  than 
weakening  of  the  propellant.  An  example  of  this  is  shown  in 
Figure  27  (the  lower  curve),  where  the  curve  drops  sub¬ 
stantially  and  indicates  reduced  values  of  a.  Such  changes 
t  too  large  to  be  due  to  burning  away  of  weaker  material. 

On  the  basis  of  this  discussion  it  appears  more  likely  that  the 
apparent  variation  in  the  propagation  constant,  if  it  occurs  in  the 
solid  phase,  is  due  to  wall  effects. 

Conclusions 

In  conclusion,  we  find  that  little  can  be  said  with  confidence 
about  the  solid-phase  behavior.  It  will  be  necessary  to  develop  an 
understanding  of  the  asymmetric  reflection  in  order  to  gain  much  infor¬ 
mation  from  experiments  similar  to  those  described.  Given  this  under¬ 
standing,  it  will  be  necessary  to  separate  wall  effects  from  changes 
in  properties;  perhaps  this  could  be  achieved  by  independently  evaluating 
the  likelihood  of  changes  in  properties. 


CHAPTER  IX 


SUMMARY  OF  CONCLUSIONS 

During  the  gas -phase  work,  data  were  obtained  which  rather  con¬ 
clusively  show  the  acoustic  behavior  in  a  single-ended  T-burner  to  be 
greatly  different  from  any  previous  description  of  the  behavior.  Most 
significantly,  the  acoustic -pressure  amplitude  was  found  to  be  lower 
at  the  burning  end  of  the  motor  by  a  factor  of  .1+5  than  the  amplitude 
at  the  opposite  end.  High  confidence  is  placed  in  this  observation. 
Consideration  of  various  observations  has  led  to  the  conclusion  that 
an  asymmetric  reflection  is  occurring  in  the  vicinity  of  the  nozzle, 
the  reflection  being  such  that  a  considerable  portion  of  a  wave  traveling 
toward  the  propellant  is  reflected  so  as  to  reinforce  the  wave  traveling 
in  the  opposite  direction.  There  appears  to  be  no  other  way  of 
explaining  the  data.  It  is  significant  that  the  observed  ehavior 
appears  to  be  associated  with  the  acoustics  of  the  cavity  ratK'  •*  than 
the  character  of  the  burning  surface. 

In  the  low-frequency  T-burner  work,  a  method  was  devised  to  esti¬ 
mate  the  response  function  for  a  burning  propellant  when  the  system  is 
actually  stable,  i.e.,  the  growth  constant  is  negative.  Itoe  method 
made  use  of  the  pressure  pulse  from  the  igniter  to  excite  oscillations 
in  the  cavity.  It  appears  likely  that  the  approach  could  be  extended 
to  much  higher  frequencies  in  order  to  study  propellants  under  stable 
conditions. 
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The  data  from  the  solid-phase  study  proved  be  intractable. 

It  was  found  that  changes  were  occurring  in  either  the  propellant 
properties,  the  boundary  conditions  at  tb*“  wal'  of  the  propellant, 
or  in  the  amount  of  asymmetric  reflection.  It  was,  however,  impossible 
to  determine  which  of  these  was  actually  changing. 
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APPENDIX  A 


ANALYTICAL  EXPRESSIONS  FOR  THE  PERTURBED 
PRESSURE  AND  VELOCITY  IN  A  T- BURNER 

Analytical  expressions  are  desired  for  the  oscillating  pressure 
and  velocity  which  may  he  used  to  infer  the  acoustic  behavior  at  the 
burning  surface  from  measurements  made  elsewhere  in  the  motor.  Rayleigh 
[45]  gives  an  analysis  of  wave  motion  in  a  cylinder,  for  which  he 
credits  Kirchhoff.  His  treatment  includes  viscous  and  thermal 
effects  for  a  cylindric ally- symmetric  cavity  with  no  mean  flow  and 
assumes  thermal  equilibrium  between  the  gas  and  the  walls.  His 
analysis  will  be  applied  to  a  T-burner  with  special  consideration 
given  to  approximations  necessary  to  obtain  the  final  results. 

It  will  be  assumed  that 

1.  the  disturbing  effect  of  the  nozzle  may  be  neglected, 

2.  the  main  body  of  the  gas  is  approximated  by  a  homogeneous 
and  non-reactive  ideal  gas, 

3.  the  effect  of  the  steep  temperature  gradient  at  the  walls 
may  be  accounted  for  in  the  parameters  of  the  final 
solution, 

4.  the  effect  of  turbulence  is  only  to  modify  the  effective 
viscosity  and  thermal  conductivity, 

5-  the  amplitude  is  sufficiently  small  to  permit  a  linear 

perturbation  analysis  (see  Appendix  C  for  the  limitations 
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of  a  linear  analysis ), 

6.  and,  in  accord  with  Rayleigh  [45],  the  dissipation  function  will 
he  neglected. 

The  effect  of  mean  flow  will  be  neglected  initially  and  then  approxi¬ 
mated  as  a  uniform  flow  (i.e.,  uniform  throughout  a  cross  section)  in 
order  to  examine  its  effect. 

Consequently,  solutions  are  required  for  the  following  equations: 
motion- -x  direction 


tfu 

p0Dt  + 


2  ~  1  d 

nV  u  *  *  ^  ^ 


(A.l) 


continuity 


P0V-  u 


(A.2) 
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Nov,  vith  the  aid  of  Equation  (A. 6)  and  the  ideal  gas  lav  in  perturbation 
form,  Equation  (A.l)  becomes 

hu-u’V2u»||  (A. 7) 


where 

P  =  (s  +  Q)  +  u"  hs  (A. 8) 

p0  v 

Since  the  solution  to  these  equations  is  lengthy  and  since  the 
reader,  with  some  persistence,  may  obtain  the  development  from  Rayleigh 
[45],  it  will  only  be  very  briefly  outlined  in  this  treatment. 

Equations  (A. 5),  (A. 6),  and  (A. 7)  are  combined  to  form  a  single 
equation  in  9.  The  method  of  combination  is  not  obvious.  Roughly, 
the  divergence  of  the  vector  equation,  of  which  Equation  (A. 7)  is 
the  x  component,  is  formed  and  Equations  (A. 5)  and  (A. 6)  are  used  to 
replace  the  divergence  of  velocity  and  s.  The  resulting  equation  is 
fourth  order  in  ©,  The  solution  to  this  equation  is  obtained  by 
assuming,  as  a  trial  solution,  the  function  Q'  which  satisfies  the 
equation 


V  Q'  =  XQ’ 


(A.9) 


A  characteristic  equation  is  obtained,  it  being 


o 

h 


^  +  h  (u’ V') 

L  0 


Vo 


+  h  (p'+u"+'D) 


X  +  h  =  0 


(A. 10) 


If  X^  and  Xg  are  taken  as  the  roots  of  Equation  (A.IO),  then 


® 


(A. ll) 
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where  Q1  and  Qg  are  the  solutions  to  Equation  (A. 9)  which  correspond 

to  X^  and  Xg,  respectively.  The  constants  and  Ag  are  to  be  determined 

by  the  boundary  conditions .  Equation  (A.io)  will  be  of  interest  later. 

The  reader  is  referred  to  Rayleigh  for  the  remainder  of  the 

1 7x 

solution.  Upon  assuming  an  exponential  dependence  on  x,  i.e.,  e  , 
he  obtains  the  following  solutions: 


u  -  AQ  + 

(A. 12) 

~  _  -.17 a 

+  a  rYh 

.  a.  rv' 

3 ~ 

+  *£  V  Xg  , 

)  3r 

(A. 13) 

p  =  pQ  (3  +  0)  (A. 14) 

v 

s  »  ©  -  “  2  0  (A. 15) 

n 


where 


J0(kQr)  eJ<Dt  sinh  J7x 
JQ(kir)  e'^a)t  cosh  j  7x 
J0(kgr)  e^a>t  cosh  J7x 


-  +  7^) 


(A.l6) 
(A. 17) 
(A. 18) 
(A. 19) 
(A.20) 
(A.21) 


in  order  to  obtain  Equations  (A.l6),  (A.17),  and  (A.l8),  an  axial 
dependence  of  e”^7x  was  incorporated,  TT  was  set  equal  to  zero  at  the 
origin,  and  a  harmonic  time  dependence  was  assumed. 
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Consider  now  the  roots  of  the  characteristic 
(A. If  .  If  a  harmonic  time  dependence  is  assumed, 
roots  are  given  explicitly  by 


equation.  Equation 
i.e.,  h«jo>,  the 


(A. 22) 


whe.  u 

A2  "  ° 

Bo  «  "  (»*'  +  H"  +  t>) 

3  H0 


C,  * 


Numerically  it  is  found  that  4  * ,  u",  o  «  o>,  e  so  that  B~  »  4  A_C_ 

3  3  3 


and 


X  « 


2A_C_ 

-B  +  (p  — 73 .. 3. . ) 

3  -  '  3  B_  1 

_ _ I _ 2 _ 

2A„ 


(A.23) 


to  an  excellent  approximation.  Thereby  it  is  found  that  the  two 
roots  are 


2 

(0 

“2 

c 


(A. 24) 


7 ^ 

»s 


(a. 25) 


These  are  valid  even  for  very  turbulent  behavior  with  effective  vis¬ 
cosities  on  -che  order  of  100  times  the  normal  viscosity.  It  is  highly 
unlikely  that  UA^C^/B.^  wo: >ld  ever  be  as  large  as  10  :  for  a  T- burner. 
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In  the  light  of  the  results  of  Equations  (A.24)  and  (A. 25)  and 
the  fact  that  7  is  on  the  order  of  unity,  it  is  found  that  k^  is  very 
small  and  and  kg  are  very  large.  Consequently  in  the  vicinity  of 
the  curved  boundary,  J^k^r)  and  J^(kgr)  may  be  adequately  approxi¬ 
mated  by  asymptotic  formulas  and  J^(k^r)  may  be  adequately  approxi¬ 
mated  by  a  truncated  series.  The  approximations  are  again  very  good. 

The  propagation  constant,  7,  may  be  obtained  by  requiring  that 
the  perturbed  velocities  and  temperature  become  zero  at  the  curved 
boundary.  Even  though  the  result  vill  not  be  expected  to  apply  exactly, 
it  is  of  interest  to  obtain  an  approximate  expression.  Through  the 
use  of  the  approximate  representations  for  the  Bessel  functions  the 
propagation  constant  is  found  to  be 


7 


+ 


(1  -  .1) 

r  'sW* 
w 


[u 


1/2  i7l/2 -r:  1/2)  +  (h 


'T 


T 


(2.7) 


This  result  may  be  compsLred  with  an  expression  obtained  by  Cantrell 
et  al.  [6J.  They  have  considered  a  somewhat  similar  problem  but  have 
included  the  effects  of  mean  flow  and  the  steep  temperature  gradient. 
By  means  of  a  variational  approach  they  have  obtained  an  approximate 
expression  for  7  as  follows : 


where  ^  here  Indicates  a  spatial  average  Except  for  the  temperature 
ratios  the  expressions  are  equivalent.  There  appears  to  he  good 
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experimental  agreement  with  Equation  (A, 27)  for  the  case  to  which  it 
really  applies  [12,45]. 

We  shall  examine  the  relative  importance  of  each  of  the  terms 
in  the  solution  for  u^  Equation  (A.  12),  because  the  complete  solution 
is  cumbersome.  If  it  is  required  that  the  perturbed  temperature  and 
radial  velocity  be  zero  at  the  curved  boundary,  two  equations  are  ob¬ 
tained  which  may  be  used  to  express  A  and  A^  in  terms  of  A^.  When  a 
particular  numerical  approximation  was  made,  it  was  found  that 


A  - 


It  is  thereby  concluded  that,  except  very  near  the  curved  wall, 

2 

u  =  -  -2—-  A1  J^k^r)  e^0*  sinh  J7x  (A. 27) 

where  ^  has  been  neglected  in  comparison  with  one.  The  perturbed 
c 

pressure,  Equation  (A.15)>  may  be  written  as 

V  =  ■rTP0  JQ(k}r)  e^at  cosh  j7x  (2.5) 

(00 

where  again  has  been  neglected  in  comparison  with  one.  The  corres- 
c 

ponding  acoustic  admittance  at  the  burning  boundary  may  then  be  written 
as 

y  =  -  tanh  J7L  (2.6) 

p0 

It  is  also  of  considerable  interest  to  investigate  the  effects 
of  aiean  flow  on  the  propagation  constant.  If  the  mean  flow  may  \e 


r 
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represented  as  u^  which  is  independent  of  position,  then,  in  operator 
notation,  the  substantial  derivative  becomes 

Dt  ■  St  +  uo  SJ  (A-28) 

±  J7x 

Now  assume  an  exponential  dependence,  e  ,  on  axial  position  and 
a  harmonic  time  dependence.  Equation  (A. 28)  thereby  becomes 

5t  -  J  (“  *  \  u0)  <A-29> 

Thus,  effectively,  the  constant  h  in  Equations  (A.  5),  (A. 6),  (A. 7), 
and  (A. 8)  has  a  slightly  changed  significance  but  the  solutions  follow 
in  the  same  manner.  The  propagation  constant  is,  however,  different 
for  each  wave.  In  the  expression  for  the  propagation  constant.  Equation 
(2.7),  co  is  replaced  by  (D+7+  u^.  Since  7  —  ^  then  os  is  approximately 
replaced  by  cn(l+M).  Thereby  Equation  (2.7)  becomes 

7+  =  +  (l  +  M)  7  (2.8) 

Thus  the  velocity  of  a  wave  moving  with  the  mean  flow  is  increased  by 
an  amount  approximately  equal  to  the  mean  velocity.  Since  the  Mach 
Number  is  small  (M  —  .01)  in  this  work.,  the  effect  on  the  propagation 
constant  should  be  correspondingly  small  and  will  be  neglected.  It 
is  worthy  of  comment  that  Cantrell  eb  al.  have  concluded  from  their 
work  that  the  only  effect  of  low  Mach  Number  mean  flow  is  to  produce 
the  steady-state  boundary  layer  [6]. 

Finally,  consider  the  transient  periods  during  the  growth  or 
decay  of  oscillations.  During  these  periods  the  angular  frequency 
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becomes  complex  and  the  amplitude  is  modulated  by  an  exponential  time 
dependence,  l.e.,  during  the  grcrwt.h  transient,  o>  is  replaced  by 
o>  -  JOf.  If  ve  let  P  denote  this  complex  frequency,  then,  for  the 
transient  periods,  the  expression  for  the  propagation  constant 
(Equation  (2.7))  becomes 


(A. 30) 


where 


w 

During  stable  oscillations. 


[»1/£  (rT1/2  -  r;  1/2)  ♦  r1/£] 

where  y  S  b-ja,  ve  find 


b 


CP1/2. 


(A.  31) 


(A. 32) 


however,  the  expression  for  a,  at  least,  is  not  quite  correct  because 
non-linear  effects  have  been  omitted.  Now  during  the  growth  transient 
we  find  from  Equation  (A.30) 


bG 


CD 
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1  + 


(A. 33) 


If  it  may  be  presumed  that  the  behavior  of  the  flame  is  linear,  the 
subscript  G  may  be  dropped  from  Equation  (A. 33)  since  b  is  expected 
to  be  essentially  the  same  during  the  growth  and  stable  periods. 

It  will  shortly  be  shown  that  a  L  =  rjar'  so  that,  from  Equation  (A. 33). 


a  +  a 
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The  decay  transient  Is  dealt  with  in  a  somewhat  different  fashion. 
After  burn-out  it  is  expected  that  the  wave  reflection  at  the  (formerly) 
burning  end  will  be  nearly  perfect  because  the  wave  then  impinges 
upon  a  rigid  steel  flange.  Consequently  it  is  required  that  the  acoustic 
velocity  (axial  direction)  be  nil  at  this  point.  This  requirement 
leads  to 

tanh  a^L  +  j  tan  bdL  =  0  (A. 3*0 


Necessarily  then 

7d  -  ^  n  **  1,2,3,  •  •  •  (A.  35) 


the  imaginary  part  being  zero.  With  this  result,  Equation  (A.30) 
may  be  used  to  obtain  the  real  and  imaginary  parts  of  i.e.. 


cu 
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where  t]  «  (■=—■)  has  been  used .  In  general,  the  fundamental 
frequencies  before  and  after  burn-out  will  not  be  the  same.  A  Fourier 
series  type  of  representation  is  required  to  meet  the  general  initial 
conditions.  However,  if,  during  the  decay,  all  frequencies  higher 
than  the  fundamental  are  filtered  out,  then  the  decay  rate  that  is 
measured  is  that  of  the  f undamental .  If  the  initial  decay  constant 
is  listed,  then  q  may  be  eliminated  between  Equations  (A. 32)  and  (A. 36) 


to  yield 
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(A.  37) 


Since  the  frequency  dependence  in  this  nonlinear  region  is  questionable 
and  since  o>  and  should  not  differ  greatly,  we  shall  approximate 


Equation  ( A . 36 )  as 


a 

*1.  =  2f  (2.1*0 

It  may  be  noted  that  if  the  frequency  dependence  of  a  were  some 
power  other  than  .5,  then  the  ratio  of  frequencies  in  Equation  (A. 37) 
would  appear  to  that  power.  Therefore,  the  non-linearities  should 
not  affect  the  approximation  in  Equation  (2.114-)  greatly. 

The  analysis  given  in  this  appendix  provides  most  of  the  relation¬ 
ships  used  In  the  gas -phase  work. 


APPENDIX  B 


ANALYTICAL  EXPRESSION  FOR  THE  RATE 
OF  OSCILLATORY  GROWTH 

There  are  several  treatments  available  from  which  one  may 
estimate  the  acoustic  admittance  of  the  burning  boundary  in  a  T- 
burner  from  growth  and  decay  constants  (e.g.,  [10]  or  Appendix  E). 
Probably  the  most  sound  and  most  elegant  treatment  is  that  of  Cantrell 
and  Hart  [7]»  Their  development  will  be  outlined  here. 

The  approach  is  to  perturb  the  energy  equation  to  second  order, 
integrate  over  the  volume,  time  average  over  harmonic  terms,  and 
then  observe  the  consequences  with  respect  to  the  growth  constant. 

The  perturbed  equations  of  motion  and  continuity  are  used  to  eliminate 
the  second-order  terms  from  the  fined  expression. 

The  energy  equation  for  isentropic  behavior  of  an  ideal  gas  is 

It  [p  «vT  +  \  u2)j  =  -  V-PU  (cv  T  +  p  +  |  u2)  (B.l)» 

The  equation  of  motion  for  irrotational  and  inviscid  flow  is 

^  2 

P  5t  +  p^^2")  +Vp  =  0  (B.2)* 

The  continuity  equation  is 

+  V*(pu)  =  0  (B.3)* 

•These  equations  may  be  found  in,e.g.,  Reference  3. 
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The  energy  equation  when  perturbed,  integrated  over  the  entire  volume, 
and  time  averaged  over  harmonic  terms  becomes 


It  /v  [po  <cvT  +  4>2  +  Oi  (°vT  +  r\ +  °2  KT  +  f  v] 

f8  “  ’ .  ”0  (h  +  +  \  (h  +  r  \  +  *2 (h  + f-  )J 


=  -  < 


where  the  definitions 


dV  > 

dS  > 
(B.4) 


C  T  +  ^  =  h  and  pu  =  m 

v  p 

have  been  used.  The  divergence  theorem  (see,  e.g.,  [6l]). 
used  in  obtaining  Equation  (B.4).  It  is  now  necessary  to  eliminate 
the  second-order  terms  from  Equation  (B.4)  by  means  of  the  equations 
of  continuity  and  motion. 

One  of  the  needed  expressions  is  obtained  by  taking  the  dot 
product  of  Sq  and  the  equation  of  motion,  Equation  (B.2), 
and  subsequently  integrating  over  the  volume  and  time  averaging.  The 
corresponding  result  is 


< 


h 

3t 


u  dV  >  = 


< 


(h 


u 


dS  > 


(b.5) 


The  divergence  theorem  has  again  been  used. 

A  second  needed  expression  is  obtained  by  multiplying  the 

2 

continuity  equation  by  (h  +  .  Upon  integration  and  time  averaging 

0 

the  result  becomes 

hT  2  j >  2 

<  J  p(h  +  £-)  d  V  >  =  ~  <  J  n  •  (h  +  )  dS  >  (B.6) 
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The  final  expression  is  obtained  by  re-expressing  the  integral 


L  p(0vT  +  ¥ 


dV 


in  terms  of  perturbed  densities,  pressures,  and  velocities,  i.e., 

2 

P  (CT  +  \ 

V 


Lb  (cvt  +  dV  >  - 


<  **  Jvl.po  -  s 


(% 


v  -2  ^2 

'T  P- 


u. 


s-  5  +  ^  '  “o +  r  )po 


+  Px  (“”p  *  uQ)  +  P2  (UQ  +  ) 


dV  > 


(B.?) 


The  substitutions 


U 


P,  ur 


and 


J_  (a  +  Irf  fill 

T--vo  _  -y  o  'pp  +  p  o  ' 

/T0  to  m0 


u„ 


/U  \  *♦  1 

=  u0  *  ^  +  ~ 

have  been  used  in  obtaining  Equation  (B.7)» 

Equations  (B.4),  (b.5)>  (B.6),  and  (B.7)  are  then  combined  to 

form 


^  a  r  !  i  pi  i  2  /-  -  \ 

<  3t  J„  I  2  —5  +  2  P0  “l  +  P1  (ul  ’  V 


V  L  -  P0C 


dV  > 


L  n 


"Vi  +  TT  uo  +  poui  (ui  •  uo 


+  Plu0  (  \  ■  V 


dS  > 


(B.8) 
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The  substitutions 

hl  -  —  “d  =.  ^ 

have  been  used  in  obtaining  Equation  (b.8). 

If  the  time  dependence  of  the  amplitude  may  be  expressed  ae  an 
exponential  in  time,  i.e.,  eGt,  and  a  «  a>,  then  the  differentiation 
with  respect  to  time  may  be  carried  out  explicitly.  For  example 


<1 i  Li 


poc 


dv  >  =  act  < 


1  pi 


\r  2  2 

v  poc 


dV  > 


(B.9) 


In  this  case  Equation  (B.8)  may  be  rewritten  as 


2a 


“o  +  Vl  '  "o)  +  Vo 


1  2  ~*  \ 

+  2  P0U1  +  P1  (ul  •  V 


dV  > 

(2.22) 


dS  > 


which  is  equivalent  to  the  result  given  by  Cantrell  and  Hart.  The  neutral 
stabi lity  criterion,  given  by  Cantrell  and  Hart,  may  be  obtained  directly 
by  setting  a  =  0,  i.e.. 


P?A  +  ~2  “o  +  Vl  '  V  +  Vo  (“x  ' 


v 


dS  >  =  0 


(B.10) 


Equation  (2.22)  is  too  cumbersome  in  its  present  form  to  be  very 
useful.  However,  if  good  approximations  to  the  acoustic  field  can  be 
made,  the  integrations  can  be  carried  out  to  obtain  a  more  useful 
expression.  For  low -velocity  mean  flow,  to  which  the  equation  applies, 
the  last  term  in  the  integrands  of  both  numerator  and  denominator  may 


'"’i 


\ 


-  n6  - 


be  neglected.  If  the  real  part  of  the  acoustic  velocity  at  the  burning 
boundary  is  not  very  small,  then  the  third  term  in  the  integrand  of  the 
surface  integral  must  be  retained.  Cantrell  and  Hart  have  assumed  that 
the  zero-order  field  is  an  adequate  approximation  and  carried  out  the 
integrations.  In  this  case  the  third  term  in  the  integrand  of  the 
surface  integral  is  negligible.  The  result  of  this  approximation  is 


a  -  a. 


DL 
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(B.12) 


where  co  =  —  and  the  definition  of  the  virtual  acoustic  admittance 
have  been  used. 

Equation  (B.12)  may  be  obtained  much  more  simply,  but  with  more 
restrictive  assumptions.  It  is,  however,  reassuring  that  the  same 
result  is  obtained,  since  this  approach  is  reasonably  general.  Equation 
(2.22)  applies  to  an  acoustic  field  which  does  not  approximate  the 
zero-order  field,  as  well  as  the  simpler  case,  and,  therefore,  is  of 
greater  interest  in  this  work. 


APPENDIX  C 


LINEAR  OR  NONLINEAR  BEHAVIOR 


Tiie  assumptions  have  been  made  in  the  analysis  of  the  gas  phase 
that  (l)  the  behavior  of  the  flame  is  linear  and  (2)  that  the  linearized 
conservation  equations  are  an  adequate  description.  Obviously  these 
assumptions  warrant  scrutiny.  Hart  and  McClure  [21]  have  indicated 
that  for  amplitudes  such  that 


Po 


<  .1 


the  gas  may  be  treated  as  linear  and  the  important  nonlinear  effects 
occur  at  the  walls.  A  criterion  such  as  this  is  necessarily  rather 
arbitrary  and  requires  a  jtidgment  as  to  how  much  error  will  be  permitted. 
An  estimate  of  this  error  will  be  made  in  this  section.  Also,  the 
reasons  for  nonlinear  effects  at  the  wall  will  be  considered. 

It  has  been  assumed  that  boundary  condition  at  the  flame  zone 
was  linear,  i.e.,  the  acoustic  admittance  was  independent  of  amplitude. 
The  justification  for  this  assumption  is  that  the  oscillations  grow 
exponentially,  i.e.,  with  the  growth  constant  being  a  true  constant, 
frequently  to  greater  amplitudes  than  those  occurring  when  our  measure¬ 
ments  were  made.  This  is  true  of  both  single-  and  double-ended 
firings.  An  exponent'  owth  Is  a  strong  indication  of  linear 
behavior  [27].  Thus,  his  is  the  case  in  the  amplitude  region 

to  which  the  analyst  ’l  1  applied,  i.e.,  where  the  measurements  were 
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made,  then  the  linear  assumption  is  thought  to  be  valid.  Because 
there  is  no  reason  to  expect  the  behavior  of  the  flame  to  change  with 
burning  time,  changes  in  amplitudes  during  burning  are  thought  to  be 
due  to  other  effects. 

.Now  it  is  advantageous  to  examine  the  isentropic  relationship. 
For  first-order  perturbations 


•"W 

£_  =  i_  £_ 

po  7t  p0 


(c.l) 


but  rigorously  the  relation  should  be  written  as 

l/>v 

f-  .  (1  +  B-)  -  1  (C.2) 

P0  Po 

For  the  case  where  7 T  =  1.25  and  p/pQ  =  .1,  the  ratio  of  the  calcu¬ 
lated  density  variations  from  these  two  equations  is  8. 00:7*93*  The 
error  Incurred  by  using  the  first-order  perturbation  approximation  is 
less  than  one  per  cent  at  this  amplitude  level. 

The  error  introduced  by  neglecting  second-order  terms  in  the 
equations  of  motion  arid  continuity  will  be  estimated  by  obtaining 
second-order  approximate  solutions  and  comparing  these  with  the  first- 
order  approximations.  The  approach  to  be  used  for  obtaining  solutions 
to  the  nonlinear  equations  is  given  by  Rayleigh  [ U5  j .  The  effects  of 
viscosity  and  mean  flow  will  be  neglected.  The  motion  will  be  assumed 
to  be  one-dimensional  and  isentropic .  For  this  case  the  equations  of 
motion  and  continuity  are 


du  ,  —  du 
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Op  ~  dp  du  ^ 

St  +  u  S  +  p  3^  a  0 


du 


(C.4) 


where  p  »  +  p.  The  first-order  approximations  are  obtained  by- 

neglecting  the  nonlinear  terms  as  small.  After  application  of  the 
boundary  conditions,  these  are  found  to  be 


=  A  e 


lent  tux 
cos  — 
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^  .  -  te-  eJmt  sir. 

1  p0°  c 


tux 


(C.5) 

(C.6) 


The  subscripts  on  £  and  u  here  indicate  the  order  of  approximation. 

These  solutions  are  now  substituted  for  the  nonlinear  terms  in  Equations 
( C . 3)  and  (C.4)  the  result  being 
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^  4  cos  — 
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(C.8) 


The  relationship  for  the  cosine  of  the  sum  of  two  angles  has  been  used 
to  obtain  Equation  (C.8).  Equations  (C.7)  and  (C.8)  may  be  combined 
to  eliminate  the  perturbed  velocity  to  yield 
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dx  c 
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P  =  2  ~2  2  cos  T 

c  P«c 


(C.9) 


where  harmonic  behavior  has  been  assumed.  The  solutions  to  these 
equations  are 
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If  A  =  |pg| ,  7^  =  1.25,  and  |  pj  /p^  «  1,  then  the  ratio  of  coefficients 
of  the  second  to  first  terms  is  .053  for  £  and  .107  for  u.  Obviously 
the  above  procedure  could  be  repeated  to  obtain  higher  order  approxi¬ 
mations;  however,  this  should  suffice. 

The  time-average  energy  density  may  be  calculated  from  the  second- 
order  approximations  as 


2p0c 


>  =  i 


P0C 


V  3  ^C2  3  QP2 


(C.12) 


Again  for  A  =■  J  p J  ,  7t  =  1.25,  and  | p  | /p0  =  .1,  the  energy  density 
dtie  to  the  second  tern  is  about  one  per  cent  of  the  total. 

On  the  basis  of  the  above  analysis  it  seems  reasonable  to  assume 
that  a  linear  description  is  adequate  for  amplitudes  such  that 

<  .1 

p0 


with  an  error  of  less  than  ten  per  cent.  The  greatest  error  appears 
to  arise  in  the  perturbed  velocity;  the  error  being  about  equivalent 
to  j Pg  /pq.  It  was  found  in  the  author’s  investigation  that  this 
ratio  was  about  <  .05  on  the  average  for  operation  in  the  single-ended 
mode. 


We  shall  now  consider  the  kinds  of  nonlinear  effects  that  might 
arise  at  the  curved  boundary.  From  the  energy  equation  it  is  observed 
that  energy  will  be  lost  at  the  wall  through  thermal  conduction  and 
work  against  viscous  forces  which  will  give  rise  to  the  terms  of  the 


type 
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(C.13) 


If  the  system  behaves  in  a  linear  fashion,  the  energy  loss  terms  depend 
upon  the  square  of  the  amplitude  [5]*  Since  both  the  viscosity  and 
thermal  conductivity  depend  upon  temperature  (absolute)  to  about  the 
.6  to  1.0  power,  they,  too,  depend  upon  the  y  coordinate.  It  is 
observed  experimentally  that  the  oscillations  enhance  the  heat  transfer 
to  the  walls  and  thus  the  mean  temperature  and  temperature  gradient 
must  depend  upon  amplitude.  Thus  it  is  easily  seen  that  the  energy 
loss  may  become  nonlinear.  It  also  appears  that  variations  in  the 
mean  temperature  and  velocity  field  may  strongly  affect  the  wall 


losses. 


APPENDIX  D 


VIRTUAL  ACOUSTIC  ADMITTANCE 


McClure  [39]  has  pointed  out  that  there  is  a  convection  of 
acoustic  energy  into  the  acoustic  system  at  the  burning  boundary. 
Effectively  the  acoustic  velocity  is  augmented  over  that  which 
would  exist  in  the  absence  of  transpiration.  In  order  to  show  how 
this  arises  we  shall  consider  a  one-dimension  description  for  isen- 
t.ropic  and  inviscid  behavior.  The  equations  of  motion  and  continuity 
are  for  this  case, 


(D.l) 


dp 


"  uo  5^ 


(D.2) 


Equations  (D.l)  and  (D.2)  may  be  combined  to  eliminate  g—;  and  terms 
2 

involving  u^  may  be  neglected  as  small.  The  corresponding  result  is 
found  to  be 


jo»0  (u  +  u0  “  0 


(D.2) 


where  a  harmonic  time  dependence  has  been  assumed.  In  the  absence  of 
mean  flow  the  quantity  in  parentheses  in  Equation  (D.3)  would  be  u. 

Thus  the  effective  acoustic  velocity  is  "u  +  Uq  ji*—  .  If  one  proceeds 
with  obtaining  a  solution  to  Equations  (D.l)  and  (D.2),  the  propogatlon 


constant  is  found  to  be 
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7+  =  +  ~  (l  +  M)  (D.4) 

By  obtaining  the  solutions  for  £  and  u  one  can  also  find  the  effective 
velocity. 

Since  the  corrections  for  mean  flow  have  been  neglected  in  this 
study,  the  acoustic  admittance  calculated  from  the  pressure  distribution 
is  expected  to  be  a  virtual  admittance.  By  assuming  that  the  propagation 
constant  was  of  the  form  of  Equation  (D.4)  and  estimating  the  Mach 
Number,  it  would  be  possible  to  obtain  the  actual  admittance.  However, 
in  the  latter  case  the  pressure  distribution  is  no  longer  of  the  cosine 
type  and  a  much  more  elaborate  calculation  would  be  necessary.  The 
small  gain  in  information  does  not  appear  to  be  worth  the  additional 
work. 
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APPENDIX  F 


THE  EFFECTS  OF  A  GRADIENT  IN  MEAN  TEMPERATURE 
ON  THE  WAVE  BEHAVIOR 


The  approach  that  will  he  used  is  to  assume  a  mean  temperature 
distribution  and  then  examine  the  wave  behavior.  We  shall  assume 
that 

1.  the  wave  behavior  is  one-dimensional, 

2.  the  gas  is  ideal  and  inviscid, 

3.  there  is  negligible  mean  flow, 

4.  and  the  wave  behavior  is  isentropic. 

Under  these  restrictions  the  equations  of  motion  and  continuity, 
Equations  (2.2)  and  (2.3),  may  be  differentiated  and  combined  to  give 


^2 — 
2-E  _ 


2  d 


2~ 


(F.l) 


dt*  dx2 
where  c  is  position  dependent. 

The  mean  temperature  distribution  should  be  obtained  from  an 
analysis  of  the  heat  transfer  problem  for  the  cavity;  however,  not 
enough  is  known  about  the  heat  loss  at  the  walls  to  carry  out  the 
analysis.  Therefore,  it  will  be  assumed  that  the  temperature  distri¬ 
bution  can  be  represented  as 

-r(L-x) 


T0  =  Tf  e 


(F.2) 
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where  r  is  an  arbitrary  constant  and  x  is  taken  as  zero  at  the  non¬ 
burning  end  of  the  motor.  A  second  distribution  will  be  assumed  as 

T0  -  (-3) 

A  solution  to  Equation  (F.l)  may  be  obtained,  with  Equation  (F.2) 
and  the  isentropic  assumption  being  used  for  c,  as 


P  -  eJa)t  [A1  JQ(0  +  Ag  Yq  (5)1  (FA) 

where 

f  2  r(L-x)/2 

^  “  c  r  e 

cf 

The  solution  can  be  completed  by  application  of  the  boundary 
conditions  to  this  result;  however,  the  behavior  of  the  solution 
is  more  easily  visualized  if  the  temperature  gradient  is  assumed  to  be 
small.  For  a  small  gradient,  T  is  small  and  £  is  large;  therefore, 
asymptotic  formulas  may  be  used  for  the  Bessel  functions.  By  means 
of  these  formulas.  Equation  (F.4)  may  be  rearranged  to  give 

p  =  1^2  eja)t  A3  cos  (£  +  A4)  (F.5) 

The  acoustic  velocity  must  be  nil  at  the  origin;  therefore, 

tan  (£c  +  A^)  =  -  gT-  (F.6) 

A4  -  -  C  -  gf-  (*-7) 

bc 

where  £c  =  CX_Q  and  advantage  has  been  taken  of  the  largeness  of  £. 

A  second  boundary  condition  is  that  pressure  at  the  origin  be  equal 
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to  the  measured  pressure  at  that  point,  i.e., 


C”  cos  (j^-)  =  jpj  eJa)t  (F.8) 

from  which  is  obtained.  Finally,  since  T  is  small,  the  exponential 
in  the  definition  of  £,  Equation  (F.4),  may  be  approximated  by  a  trun¬ 
cated  expansion.  The  acoustic  pressure  may  then  be  approximated  as 

7-  |Fg|  e1*/1*  cos  [f  (x  -  ^)]  (F.9) 


where 


c 

c 


e-rL/2 


From  this  result  it  is  seen  that  the  amplitude  increases  with  mean 
temperature .  Thus,  if  the  mean  temperature  is  higher  at  the  burning 
end  than  at  the  non-burning  end  of  the  motor,  the  amplitude  should 
be  greatest  at  the  propellant  end.  The  same  conclusion  was  reached 
when  Equation  (F.l)  was  solved  with  Equation  (F.2)  as  the  temperature 
distribution  on  an  analog  computer. 

The  second  temperature  distribution.  Equation  (f.3)>  also  leads 

In  this  case,  the  solution  is 

(r)  +  c2  j.u(c)]  (f.io) 

e-rx/2 


to  an  equation  of  the  Bessel  type. 


~  _  Ojo>t  r 
p  -  e  [Cx  J ^ 


2  cu_ 
J  ~  _ 

1  Cf 


V 


CD 

Cf  r 


I  K1^ 


where 
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If  r  is  large,  so  that  the  temperature  distribution  approximates  the 
thermal  boundary  layer  at  the  cold  end,  the  Bessel  functions  may  be 
approximated  by  truncated  series.  In  this  case  Equation  (F.IO)  is 
approximated  by 


= 


cox 

“J 


°3e 


1  - 


(0.2 
'frK 


1/2 


-rx 


4  (l  +  u) 


+  c. 


0>X 

c 


CD. 

C. 


1  - 


Mi  -  u) 


(F.ll) 


Here  it  is  observed  that  as  long  as  the  wave  length  is  large  compared 
with  the  thickness  of  the  boundary  layer,  the  terms  within  [ ]  in  this 
equation  are  nearly  unit.  Therefore,  for  a  thin  boundary  layer  the 
solution  is  the  same  as  that  obtained  if  the  layer  is  neglected  and  the 
layer  has  no  effect. 

If  T  is  small  and  asymptotic  formulas  are  used  to  approximate 
the  Bessel  functions  in  Equation  (F.IO),  an  equation  of  the  same  form 
as  (F.9)  is  obtained.  Therefore,  this  analysis  also  leads  to  the  con¬ 
clusion  that  a  higher  tuean  temperature  at  the  burning  end  will  cause 
a  higher  pressure  amplitude.  These  last  two  conclusions  were  also 
verified  by  solving  Equation  (F.l),  vith  (F.3)  as  the  temperature 
distribution,  on  an  analog  computer. 


APPENDIX  G 


APPARATUS  AND  PROCEDURES 

In  this  appendix  the  details  of  M  e  experimental  system  will  be 
discussed. 

T- Burner :  Gas  Phase  Study 

The  T-burner  used  in  the  gas -phase  work  is  'hown  in  Figure  28.  A 
sketch  of  the  configuration  is  shown  in  Figure  1.  The  motor  was  con¬ 
structed  entirely  from  steel,  except  for  a  replaceable  carbon  nozzle. 

The  central,  or  nozzle,  section  was  coupled,  and  sealed  by  means  of 
O-rings,  to  approximately  equal  lengths  of  tubing  (l.27-cm.  wall 
thickness,  nominally),  the  ends  of  whlc.  in  turn  were  sealed.  The  inside 
diameter  was  3*96  cm.  throughout.  End  closures  accommodated  the  pressure 
transducers  and  were  bonded  firmly  in  place  for  each  firing  with  an 
epoxy  adhesive.  Heavy  steel  flanges  were  used  to  rigidly  clamp  the 
motor  together  for  a  firing. 

Propellant  samples  were  cut  to  the  desired  length  (l.27  cm., 
nominally)  from  cylindrical  castings.  Before  cutting,  these  grains 
were  prepared  by  scraping  the  sides  to  remove  a  parting  agent  acquired 
during  casting  and  leaching  to  inhibit  the  sides  and  provide  a  clean 
surface  for  bonding  purposes.  The  ends  of  the  cut  samples  were  scraped 
with  a  razor  blade  to  remove  irregularities.  The  curved  wall  of  the 
sample  was  coated  with  epoxy,  as  was  the  top  of  the  end  closure  except 
for  the  surface  of  the  transducer.  A  small  amount  of  epoxy  was  spread 
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Figure  28.  Two  views  of  the  T-burner  used  in  the  gas -phase  work. 
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around  the  inside  of  the  end  of  the  tube  because  of  a  tendency  for  the 
adhesive  to  be  "viped-off"  the  sample  upon  insertion.  The  sample  was 
then  placed  on  the  end  closure  and  inserted  into  the  end  of  the  motor. 
Firm  transducer  contact  was  obtained  by  placing  a  rubber  stopper  and 
weight  on  the  free  surface  of  the  sample  while  the  epoxy  cured. 

Converging  nozzles  were  made  from  carbon  rod  stock,  the  throat 
diameter  being  determined  by  the  mean  pressure  desired  and  length  of 
motor  to  be  used. 

T-Burner:  Solid-Phase  Study 

The  motor  used  for  the  solid-phase  work  was  constructed  by  and 
is  described  by  Coates  [11 ].  A  sketch  of  the  configuration  is  given 
in  Figure  25.  Except  for  the  piston  arrangement  used,  to  advance  the 
propellant,  the  motor  was  the  same  as  the  motor  for  the  gas -phase  work. 
One  end  closure  was  replaced  by  the  piston  which  advanced  the  propellant 
so  that  the  burning  interface  remsiined  at  a  nearly  constant  position. 

The  piston  was  sealed  by  means  of  an  O-ring.  During  burning  the 
position  of  the  interface  was  sensed  by  a  thermocouple  in  the  wall  of 
the  motor.  The  output  of  the  thermocouple,  in  conjunction  with  a  con¬ 
trol  unit,  was  used  to  control  the  rate  of  propellant  movement.  The 
piston  was  hydraulically  driven,  the  rate  being  varied  by  an  electri¬ 
cally  operated  control  valve. 

Propellant  samples  were  cut  to  the  deBired  length  from  castings 
which  fit  into  the  motor  with  a  low  wall  clearance  (  .01  cm.).  These 

samples  were  securely  bonded  to  the  face  of  the  driving  piston.  Firm 
transducer  contact  was  obtained  by  placing  a  weight  on  top  of  the  sample 
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while  the  adhesive  cured.  The  curved  wall  of  the  sample  was  liberally 
coated  with  hydrocarbon  grease  before  assembling  the  motor. 

Igniter  System 

The  igniter  system  was  composed  of  a  short  piece  of  Pyrofuze  wire 
(3-5  cm..)  and  a  small  amount  of  F-propellant  sawdust.  Pyrofuze  wire 
(Pyrofuze  Corporation,  Mount  Vernon,  New  York)  is  a  bi-metallic  wire 
which  reacts  exotherm! cally  when  heated.  Short  copper  leads  were 
attached  to  the  Pyrofuze  wire  which  was  subsequently  sprayed  with  an 
electrically  insulating  resin.  The  Pyrofuze  wire  and  propellant  sawdust 
were  placed  adjacent  to  the  propellant  in  the  motor  and  the  copper 
leads  passed  through  a  rubber  stopper  inserted  in  the  nozzle.  Ignition 
was  achieved  by  electrically  heating  the  Pyrofuze.  The  function  of  the 
rubber  stopper  was  to  carry  the  igniter  leads  out  of  the  motor  as  the 
pressure  rose. 

Ins  trumentation 

The  instrumentation  is  summarized  by  Figure  29  in  which  a  block 
diagram  of  the  electronic  system  is  shown.  The  system  shown  is  that 
most  frequently  used  in  the  gas  phase  work  but  is  little  different 
from  that  used  at  other  times.  Switch  (h)  permitted  a  reference 
signal  to  be  recorded  following  the  firing  data  without  stopping  the 
tape  recorder.  The  RC  filters  and  voltage  dividers  (5)  provided  the 
necessary  separation  of  the  oscillating  component  from  the  total  sig¬ 
nal  and  permitted  the  adjustment  of  the  recorded  voltage  to  the  proper 
level.  During  data  acquisition  the  tape  recorder  (j)  was  operated  at 
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Figure  29-  Block  diagram  of  instrumentation. 

1.  Kistler  Model  401  pressure  transducers  (or  Model  6l6  water-cooled 
transducer) 

2.  Kistler  Model  5&5  charge  amplifier 

3-  Heath-kit  Model  AG-9A  audio  oscillator 
U.  Double-pole-double- throw  switch 
5*  RC  filters  and  voltage  dividers 
6.  Tektronix  Model  122  preamplifiers 

7>  Precision  Instrument  Model  207  tape  recorder  with  two  direct  and 
one  FM  record -reproduce  channels 

8.  Tektronix  Model  535  oscilloscope  with  Type  M  plug-in  unit 

9-  Hewlett-Packard  Model  150A  oscilloscope 


REPRODUCE  SYSTEM 


Figure  29  (continued) 

5.  RC  filters  and  voltage  dividers 

6.  Tektronix  Model  122  preamplifiers 

7*  Precision  Instrument  Model  207  tape  recorder 
8.  Tektronix  Model  535  oscilloscope 

10.  Spencer-Kennedy  Laboratories  Model  302  variable  electronic  filter, 

dual  channel 

11.  Ad-Yu  Model  405  electronic  phase  meter 
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60  inches  per  second  with  two  AC  and  one  DC  signal,  being  recorded. 

The  original  data  were  also  recorded  photographically  on  the  two 
oscilloscopes  (8,9)* 

The  data  were  reproduced  from  the  tape  recorder  (7)  at  7  l/2  inches 
per  second  (usually).  The  filter  (lO)  reduced  the  distortion  and  noise 
levels.  Oscillatory  amplitudes  and  frequencies  were  measured  by  dis¬ 
playing  the  signals  on  the  oscilloscope  (8)  and  photographing  the  traces. 
The  phase  meter  (ll)  was  used  to  measure  phase  angles  between  the  data 
signals,  the  meter  output  being  displayed  on  the  oscilloscope  simul¬ 
taneously  with  the  oscillatory  envelopes.  The  output  of  the  tape 
recorder  was  amplifed  (5*6)  to  the  necessary  input  level  for  the  phase 
meter. 

The  experimental  system  used  was,  for  the  most  part,  satisfactory. 
Most  of  the  difficulties  that  arose  were  due  to  failures  of  the  instru¬ 
ments  or  wiring,  rather  than  any  weakness  of  the  overall  system. 


APPENDIX  H 


PROPELLANT 

The  propellant  used  in  substantially  all  of  the  work  reported 
here  is  known  as  Utah  F.  It  was  produced  in  small  batches  in  this 
laboratory.  The  tion  was: 

1.  40  weight  per  cent  coarse  ammonium  perchlorate  (American 
Potash  size:  -48  +100), 

2.  40  weight  per  cent  fine  ammonium  perchlorate  (American 
Potash  size:  50  per  cent  less  them  15  microns), 

3.  2  weight  per  cent  copper  chromite  (Harshav  Chemical,  Cu-0202  p), 

4.  15-3  weight  per  cent  polybutadiene-acrylic  acid  copolymer 
(Thiokol), 

5.  2.7  weight  per  cent  epoxy  resin  (Shell,  Epon  828). 

The  propellant  was  mixed  for  60  minutes  at  a  pressure  of  7  to  8  mm.  Hg 
absolute  and  a  temperature  of  65 °C,  nominally.  It  was  cast,  at  the 
same  pressure.  Into  a  cylindrical  mold.  It  was  cured  for  7  days  at 
80 °C,  nominally. 

F  propellant  has  the  following  pertinent  properties: 

hoc 

1.  Linear  burning  rate,  rQ  ■  .394  (pQ)'  7  cm/sec  where  pQ  is 

in  atmospheres  (absolute).* 

# 

The  rate  expression  was  determined  from  strand  bomb  data.  These 
data  may  be  found  in:  Beckstead,  M.  W.,  unpublished  Ph.D.  thesis, 
University  of  Utah,  Dept,  of  Chemical  Engineering  (Salt  Lake  City, 

June,  1985)- 
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2.  Density,  P  ■  1.63  gm/cra^ 

3»  Adabiatic  flame  temperature  is  2600°K  at  68  atm. 

4.  Calculated  ratio  of  heat  capacities,  7^=1 -24  at  68  atm.  and 
2600°K. 

5-  Calculated  molecular  weight  is  23-5  gm/mole  at  68  atm.  and 
2600°K. 

6.  Calculated  isentropic  velocity  of  sound  is  1.0?  (lO^)  cm/sec 
at  68  atm.  and  2600°K. 

7.  Heat  capacity,  C  *  . 316  cal/gm  K. 

P 

p 

8.  Thermal  diffusivity  is  .001697  cm  /sec. 

9-  Heat  of  combustion  (old  P  propellant)  is  2.563  Kcal/gm. 
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VALIDITY  OF  MEASUREMENTS 

The  experimental  observation  that  the  pressure  amplitude  at  the 
burning  surface  vas  much  less  than  that  at  the  non-burning  end  was  a 
source  of  concern.  Consequently  a  rather  thorough  investigation  was 
carried  out  in  an  effort  to  explain  it. 

It  is  felt  that  the  low  measurement  was  not  due  to  the  instru¬ 
mentation  for  the  following  reasons : 

1.  The  transducers  and  information  channels  were  interchanged 
frequently  and  the  observation  was  unaffected. 

2.  The  transducer  calibrations  were  frequently  checked. 

3.  A  reference  signal  was  recorded  immediately  after  each 
firing.  The  reference  signal  was  switched  into  both  infor¬ 
mation  channels  and  passed  through  each  complete  system 
with  the  exception  of  the  Kistler  charge  amplifiers  and 
transducers.  The  signal  came  from  a  single  audio  oscillator 
which  was  set  at  a  frequency  close  to  that  of  the  data. 

The  reference  was  subsequently  used  to  compensate  for  any 
differences  in  the  gain  of  the  two  electronic  systems. 

4.  The  thermal  effects  on  the  transducers  should  be  small.  The 
output  of  the  Kistler  401  transducer  increases  with  tempera¬ 
ture  at  the  rate  of  .02  per  cent  per  degree  K.  Thus,  unless 
the  temperature  of  the  transducer  increases  enormously,  the 
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effect  vill  be  insignificant .  The  DC  level  of  the  output 
drops  with  heating  due  to  reduced  stress  on  the  crystal, 
the  crystal  being  mounted  tinder  stress  in  the  transducer. 
From  the  kinds  of  DC  shifts  encountered  during  operation, 
it  is  estimated  that  the  temperature  rose  about  20°K. 

5-  The  401  transducer  has  a  recess  around  the  sensing  diaphragm 
on  the  face  of  the  transducer.  Normally  this  recess  was 
filled  with  silicone  grease.  It  was  thought  that  under 
pressure  the  propellant  might  be  inadequately  supported  in 
this  region  and  cause  a  slightly  dome-like  surface  over  the 
sensing  diaphragm.  In  order  to  test  for  this  problem  the 
recess  was  filled  with  epoxy  and  a  firing  was  made.  The 
measured  relative  amplitude  was  .455* 

6.  As  a  check  on  the  ability  to  accurately  measure  a  pressure 
through  the  propellant,  a  thin  sample  of  propellant  (-35  cm.) 
was  bonded  into  a  short  test  section  in  the  configuration 
used  during  a  firing.  The  test  section  was  then  adapted 
to  a  dead-weight  gauge  tester  and  the  transducer  output  was 
measured  under  these  conditions.  The  mean  pressure  measured 
was  about  four  per  cent  higher  than  that  measured  without 
the  propellant.  The  small  increase  in  output  is  attributed 
to  an  increased  effective  area  for  the  sensing  diaphragm, 
i.e.,  the  recess  causes  more  of  the  load  to  be  supported  by 
the  sensing  diaphragm  than  would  occur  if  the  measurement 
were  made  in  a  fluid. 
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7*  Since  no  technique  was  available  for  dynamically  calibrating 
the  transducers,  a  crude  dynamic  comparison  was  devised. 

Two  transducers  were  mounted  in  an  assembled  motor  with  no 
propellant.  An  igniter  and  a  small  amount  of  propellant 
sawdust  were  ignited  near  the  midpoint  of  the  motor.  The 
outputs  of  the  two  transducers  were  then  compared.  A  rate 
of  pressure  rise  equivalent  to  200  cps  was  obtained.  No 
difference  in  output  could  be  detected. 

8.  A  great  dead  of  care  was  exercised  in  assembly  of  the  pro¬ 
pellant  sample  and  transducer  to  ensure  proper  contact  with 
the  sensing  diaphragm.  The  transducer  was  mounted  so  that 
the  sensing  diaphragm  protruded  about  .08  cm.  above  the  end 
plate.  This  allowed  for  a  reasonable  glue  line  of  epoxy. 
After  assembly  the  propellant  vas  held  in  firm  contact  by  a 
small  weight  resting  on  a  rubber  stopper.  It  was  felt  that 
proper  contact  was  thus  obtained.  If  poor  contact  had  been 
a  problem,  it  should  have  appeared  through  distortion  of  the 
output  signal.  A  comparison  of  the  waveforms  at  the  two 
ends  of  the  motor  was  favorable. 

9.  Since  Kistler  transducers  are  acceleration  sensitive,  it 
was  felt  that  this  might  be  a  problem.  An  output  due  to  a 
vibrational  acceleration  could  be  in  phase  with  the  pressure 
and  thus  augment  the  amplitude  of  the  total  output.  A  test 
for  such  effects  was  devised.  The  cold  end  of  the  motor 
was  closed  by  bonding  a  .  3-cra.  steel  plate  over  the  end. 
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The  motor  was  assembled  as  usual  except  the  transducer  in 
the  cold  end  was  not  mounted  in  contact  with  the  .3-cm. 
plate.  The  transducer  was  thereby  prevented  from  exposure 
to  the  gases  but  still  exposed  to  any  accelerations  present. 
When  a  firing  was  carried  out  in  this  manner  the  output  of 
this  transducer  showed  no  signal  above  the  noise  level 
even  though  strong  oscillations  were  present  in  the  motor. 

10.  Because  the  output  of  the  Tektronix  122  preamplifiers  used 
in  each  information  channel  is  not  flat  in  the  frequency 
range  of  interest,  a  comparison  of  the  relative  output  was 
made.  A  signal  was  passed  through  each  amplifier-filter 
combination  over  the  frequency  range  of  1.0  to  6.0  Kc. 

The  variation  in  relative  output  was  completely  negligible. 

On  the  basis  of  the  above  instrumentation  evaluation  it  was  con¬ 
cluded  that  the  relative  amplitude  measured  was  that  experienced  by  the 
transducers. 

One  area  of  obvious  concern  is  the  neglect  of  attenuation  in 
the  propellant,  i.e  ,  the  assumption  that  the  dynamic  pressure  measured 
under  a  thin  slab  of  propellant  is  equivalent  to  that  at  the  surface 
of  the  propellant.  The  validity  of  this  assumption  was  checked  by  two 
routes.  The  first  was  to  measure  the  pressure  on  the  non-burning  end 
of  the  motor  through  a  .2-cm.  slab  of  propellant  whose  surface  was 
inhibited  to  prevent  burning.  The  average  relative  amplitude  for  two 
runs  made  in  this  manner  was  .460.  Ti.-1  second  route  was  to  examine 

the  decay  of  oscillations  after  burn-out.  It  was  found  that  the  decay 
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occurred  smoothly  and  that  a  non-unity  relative  amplitude  persisted 
until  the  signals  reached  the  noise  level.  It  was  thereby  concluded 
that  the  pressure  measured  was  equi\  lent  to  that  at  the  propellant 
surface  during  burning. 

The  effects  of  the  boundary  layers  at  each  end  of  the  motor  on 
the  pressure  measurement  must  be  considered  in  order  to  complete  the 
examination.  At  the  non-burning  end  it  is  expected  that  there  is  a 
steep  temperature  gradient  between  the  transducer  and  the  main  body 
of  the  gas .  Because  the  mean  flow  is  expected  to  be  small,  the 
effects  of  this  thermal  boundary  layer  can  be  approximated  analytically 
without  great  difficulty.  In  Appendix  F  an  analytical  solution  is  ob¬ 
tained  for  a  mean  temperature  distribution  that  is  thought  to  approxi¬ 
mate  that  due  to  the  boundary  layer.  This  solution  (Equation  (F.ll)) 
indicates  that  the  pressure  amplitude  is  unaffected  as  long  as  the 
boundary  layer  Is  thin  compared  to  the  wave  length.  It  is  expected  that 
the  boundary  layer  may  be  neglected  in  the  problem  at  hand. 

The  effect  of  the  boundary  layer  associated  with  the  flame  is 
considerably  more  complicated  to  treat.  It  is  usually  argued  that  the 
amplitude  must  be  unaffected  in  traversing  the  flame  zone  because  the 
flame  zone  Is  much  thinner  than  the  wave  length.  However,  because  the 
flame  zone  is  very  complex.  It  is  difficult  to  assure  oneself  that  this 
argument  is  adequate.  Consequently  a  rather  extensive  approximate 
analytical  treatment  of  the  wave  behavior  in  the  flame  zone  was  developed. 
Because  of  its  complexity  it  will  be  given  in  a  separate  appendix 
(Appendix  j).  The  approach  was  to  perturb  the  one-dimensional  equations 


of  energy,  motion,  and  continuity;  with  mean  flow  and  mean  temperature 
terms  being  retained.  The  resulting  equations  were  combined  and  solved 
numerically.  The  solutions  showed  no  indications  of  any  serious  re¬ 
ductions  in  amplitude  in  passing  through  the  flame  zone  to  the  solid. 

It  was  concluded  that  pressure  amplitude  on  the  gas  side  of  the  flame 
zone  was  essentially  the  same  as  at  the  solid  surface. 

Consideration  was  also  given  to  the  effects  of  gradient  in  mean 
temperature  along  the  motor.  Analytical  solutions  were  obtained  for 
the  wave  behavior  with  assumed  mean  temperature  distributions  (see 
Appendix  F).  This  analysis  showed  the  pressure  amplitude  to  be  higher 
at  the  hottest  end  of  the  motor.  Since  the  non-burning  end  cannot  be 
hotter  than  the  burning  end,  it  was  concluded  that  a  gradient  in  mean 
temperature  would  oppose  the  low  amplitude  at  the  propellant  end. 

Several  firings  were  made  with  a  newly  acquired  water-cooled 
transducer  mounted  in  the  curved  wall  of  the  motor.  The  intent  was 
to  compare  the  amplitude  measured  by  this  transducer  with  the  ampli¬ 
tude  inferred  from  end  measurements.  The  data  obtained  were  rather 
erratic  so  that  a  good  comparison  could  not  be  made;  however,  the  ampli¬ 
tudes  measured  in  the  wall  and  near  the  propellant  were  compatible  with 
the  amplitudes  measured  under  the  propellant. 

Two  firings  were  made  with  the  nozzle  shifted  toward  the  propellant 
from  the  midpoint.  The  damping  effects  of  the  nozzle  should  be  mini¬ 
mized  when  the  locations  of  the  nozzle  and  the  pressure  node  coincide. 

The  data  from  these  firings  indicated  that  the  amplitude  was  reduced 
when  the  nozzle  was  closer  to  the  propellant  than  x/L  »  .65.  When  the 
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nozzle  was  located  In  the  range  x/L  =  -55  to  .65,  the  amplitude  appeared 
to  be  normal.  In  both  cases  the  relative  amplitude  was  about  the  value 
observed  before.  The  average  relative  amplitude  was  .475-  Distortion 
seemed  to  be  higher  than  usual.  These  measurements  indicate  that  the 
pressure  node  is  located  near  the  midpoint  of  the  motor. 

The  conclusions  of  importance  may  be  summarized  as  follows: 

1.  The  measured  amplitude  ratio  exists  within  the  motor.  In 
particular,  it  is  not  introduced  by  the  instrumentation; 

it  is  not  introduced  by  either  the  thermal  boundary  layer  at 
the  cold  end  or  the  boundary  layer  associated  with  the  flame; 
and  it  Is  not  due  to  the  propellant  through  which  the  low 
amplitude  is  measured. 

2.  A]. though  the  gas  may  be  cooler  at  the  non-burning  end  of  the 
motor,  the  effect  is  in  opposition  to  that  measured. 

3-  Pressure  measurements  made  in  the  curved  wall  of  the  motor 
and  near  the  propellant  were  compatible  with  those  made 
under  the  propellant. 

4.  Measurements  made  with  the  nozzle  shifted  toward  the  pro¬ 
pellant  from  the  midpoint  of  the  motor  indicate  that  the 
pressure  node  is  in  the  vicinity  of  the  midpoint. 

5.  The  low  relative  amplitude  persists  through  most  of  the  decay 


transient. 


APPENDIX  J 


ANALYSIS  OF  THE  WAVE  BEHAVIOR  IN  THE  FLAME  ZONE 

At  one  point  in  the  gas -phase  work  the  data  appeared  to  indicate 
that  the  pressure  amplitude  was  significantly  lower  on  the  propellant 
side  of  the  flame  zone  than  on  the  gas  side.  Consequently,  an  elabo¬ 
rate  analysis  of  the  wave  behavior  in  the  flame  zone  was  developed. 

Later  consideration  of  the  data  showed  that  the  observed  amplitude 
behavior  could  not  be  due  to  the  flame  zone  and,  therefore,  the 
analysis  was  really  unnecessary.  However,  the  analysis  will  be 
outlined . 

The  approach  was  to  perturb  the  equations  of  motion,  <.  ontinuity, 
and  energy  for  one-dimensional  and  inviscid  behavior.  The  steady-state 
equations  were  subtracted  from  the  perturbed  equations  and  the  resulting 
equations  combined  to  for  t  ' Ifferential  equation  in  the  perturbed 
pressure.  Because  ol  the  -lability  to  calculate  the  mean  temperature 
distribution,  various  disorib  tions  were  assumed.  The  differential 
equation  was  solved  numerically.  It  was  assumed  that  the  analytical 
expression  for  the  acoustic  pressure,  Equation  (2.5),  was  sufficient 
outside  the  flame  zone.  The  analytical  solution  was  used  -.0  calculate 
initial  values  for  the  numerical  integration.  The  numerical  solution 
then  gave  the  pressure  amplitude  at  the  propellant  surface,  i . e . ,  under 


the  flame  zone. 
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The  perturbed  equations  of  motion  and  continuity,  subject  to  the 
above  conditions  and  after  subtraction  of  the  steady- state  equations, 
are 

motion  -  x  direction 


du 

Pn  JT  +  JET  + 


Su 


0  ~ 


du„ 


(J.l) 


continuity 


(J.2) 


dp 

The  quantity  has  been  neglected  since  it  can  be  shown  to  be  small 
2 

when  PqUq  «  pQ  and  the  variation  in  molecular  weight  is  neglected. 

A  harmonic  time  dependence  was  assumed.  Equations  (j.l)  and  (j.2) 

may  then  be  combined  to  solve  for  \x.  This  expression  is  differentiated 
Su 

to  obtain  ~  .  Both  of  these  expressions  may  then  be  inserted  into 
Equation  (j.2)  to  obtain  a  differential  equation  in  p'  and  p*.  It  was 
then  necessary  to  appeal  to  the  energy  equation  to  relate  p"  and  p*. 

In  perturbation  from  the  energy  equation  is 


pocp  I  +  cCp 


DT0  D|- 
Dt  "  Dt 


(J.3) 


where  the  steady-state  equation  has  been  subtracted  and  the  reaction 
terms  and  perturbed  conduction  terms  neglected.  A  harmonic  time 
dependence  and  a  constant  molecular  weight  sire  assumed.  Equation  ( J . 3) 
then  yields  an  expression  in  terms  of  the  zeroth  and  first  derivatives 
of  p*  and  p”.  This  expression  is  then  combined  with  the  expression 
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obtained  from  the  equations  of  motion  and  continuity  to  eliminate  p. 

A  third-order  equation  is  obtained  with  variable  and  complex  coefficients. 
The  coefficients  are  very  complicated. 

The  third-order  differential  equation  was  solved  with  a  Runga- 
Kutta-Gill  integration  technique.  The  mean  temperature  wan  assumed  to 
be  described  by  functions  of  the  exponential,  sine,  parabolic,  and 
cubic  types.  In  addition,  the  effects  of  the  chemical  reactions  were 
incorporated  by  making  some  very  crude  approximations. 

A  number  of  variations  were  tried  and  extensive  computer  time 
used  in  carrying  out  the  analysis;  however,  in  no  case  were  we  able 
to  predict  a  pressure  amplitude  under  the  flame  zone  that  was  signifi¬ 
cantly  less  than  that  outside  the  flame. 


APPENDIX  K 


ESTIMATION  OF  THE  ISENTROPIC  VELOCITY  OF 
SOUND  FROM  THE  NOZZLE  CONDITIONS 


The  isentropic  velocity  of  sound  in  the  main  body  of  the  gas 
may  be  quite  simply  estimated  from  the  nozzle  conditions  if  it  is 
assumed  that 

1.  the  flow  through  the  nozzle  is  steady, 

2.  the  conditions  in  the  main  body  of  the  gas  approximate  the 
stagnation  conditions  for  the  gas, 

3-  the  flow  is  isentropic, 

4.  and  the  contraction  coefficient  for  the  nozzle  is  unity. 

All  four  assumptions  should  be  quite  good.  Since  the  Mach  Number 
at  the  throat  of  the  nozzle  must  be  unity,  the  gas  temperature  at  the 
throat  is  related  to  the  gas  temperature  in  the  main  body  of  the  gas  as 


L0t 


(K.l) 


V1 


Now,  if  the  mass  generation  rate  at  the  propellant  is  equated  to  the 
mass  flow  rate  through  the  nozzle,  then 


(5*1 

2 

po 

1  g  1 

K  / 

r_o  'T 

C  a 

IV1/ 

7t+1 

2(7j-i; 


(K.2) 


Everything  on  the  right  side  of  this  equation  is  known,  can  be  measured, 
or  can  be  estimated  reasonably  well.  The  linear  burning  rate  expression 
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(see  Appendix  H)  may  be  substituted  for  so  that  c  may  be  expressed 
as  a  function  p^  only,  the  function  being  a  proportionality  to 


Equation  (K.2)  was  used  to  estimate  the  velocity  of  sound.  For 
firings  in  the  2.5  to  4.5  kilocycle  range,  the  throat  diameter  was 
.521  cm.,  whereas  in  the  1  to  2.5  kilocycle  range  the  throat  diameter 
was  .510  cm.  At  68  atm.  and  2600°K,  7^  =  1.24.  Since  the  pressure 
and  temperature  were  actually  somewhat  lower  than  these  values, 

7t  =  1.26  was  used.  The  calculations  gave 

c  =  .967  (105)  cm/sec  D  =  .521  cm. 

c  =  .911  (105)  cm/sec  =  .  510  cm. 

If  7t  -  1.24  is  used,  c  is  lower  by  a  factor  of  .985* 

The  above  calculated  values  were  used  in  Chapter  IV. 
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APPENDIX  M 


PARAMETRIC  CURVES 

This  appendix  contains  the  parametric  curves  developed  from 
Equation  (6.15)  which  were  used  in  the  solid-phase  work.  Use  of  the 
curves  requires  a  knowledge  of  the  ratio  of  amplitudes  and  the  phase 
angle  between  the  oscillating  pressures  at  the  ends  of  the  propellant. 
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APPENDIX  N 


TABLES  OP  DATA 
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TABLE  3 

GROWTH  CONSTANT-FREQUENCY  DATA  FOR  OLD  UTAH-F  PROPELLANT 


Nominal  mean  pressure,  14.5  atm. 
Nominal  initial  sample  length,  1.3  cm. 


Run  Number 

a 

(sec-1) 

f 

(kc.) 

Notes 

Double-ended 

0-40 

63 

2.3 

(a) 

firings 

0-44 

72 

2.0 

0-45 

169 

6.7 

0-53 

167 

6.9 

(b,  c) 

0-57 

43 

1.0 

0-58 

64 

1.7 

0-59 

90 

1.7 

0-04 

73 

1.7 

(c) 

0-66 

78 

2.2 

(c) 

0-67 

91 

2.9 

(c) 

0-68 

84 

2.9 

(a) 

0-69 

81 

2.2 

(c) 

0-70 

79 

2.3 

0-71 

102 

2.9 

(b) 

0-72 

43 

1.1 

0-73 

31 

•  70 

(b) 

0-74 

37 

.73 

0-75 

31 

.98 

0-76 

183 

7.0 

(b,  c) 

0-77 

117 

4.7 

O-78 

110 

4.8 

(b,c) 

0-79 

167 

7-0 

0-82 

42 

.78 

(c) 

0-83 

113 

4.8 

0-101 

64 

1.7 

(b) 

0-102 

47 

1.1 

Single -ended 

0-103 

14 

.96 

firings 

0-105 

27 

1.9 

0-107 

48 

3.8 

0-108 

52 

5-1 

0-109 

56 

4.7 

0-110 

52 

3.5 

0-111 

23 

2.4 

0-112 

25 

1.9 

0-113 

21 

1.6 

0-114 

15 

•  99 

0-115 

9 

.70 
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TABLE  3  (continued) 


Run  Number 

a 

(sec"l) 

f 

(kc.) 

Notes 

0-124 

25 

1.7 

0-128 

9 

.63 

(b) 

0-131 

20 

1.7 

0-149 

64 

9-8 

0-151 

56 

10.2 

Notes : 

(a)  High  mean  pressure,  the  mean  pressure  exceeded  the  nominal 
value  of  14.5  atm.  by  more  than  3*5  atm. 

(b)  Low  mean  pressure,  the  mean  pressure  was  less  than  the 
nominal  value  of  14.5  atm.  by  more  than  3*5  atm. 

{  Marginal  mean  pressure,  the  mean  pressure  trace  Indicated  a 
non -uniform  burning  area. 
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TABLE  4 

GROWTH  CONSTANT-FREQUENCY  DATA  FOR  NEW  UTAH-F  PROPELLANT 

Nominal  mean  pressure,  14.5  atm. 

Nominal  initial  sample  length,  1.3  cm. 


Run  Number 

a 

(sec”^) 

f 

(kc.) 

Notes 

Double -ended 

firings 

0-264 

137 

3.2 

0-265 

105 

3.3 

Single-ended 

0-208 

49 

4.9 

firings 

0-209 

48 

5.0 

0-210 

73 

4.3 

(a) 

0-218 

49 

3-5 

0-219 

44 

5.0 

0-220 

37 

3.5 

0-221 

55 

4.3 

0-223 

26 

4.1 

0-224 

54 

5.0 

0-225 

50 

4.7 

0-226 

60 

4.9 

0-227 

43 

1.8 

0-228 

38 

2.3 

0-230 

42 

1.7 

0-231 

37 

1.7 

0-232 

38 

2.3 

0-233 

37 

1.7 

0-234 

38 

2.2 

0-235 

35 

2.2 

0-236 

34 

1.8 

0-237 

34 

1.7 

0-238 

4l 

2.3 

0-240 

39 

2.2 

0-241 

4l 

1.7 

0-242 

35 

2.3 

0-243 

75 

5.1 

0-244 

75 

4.5 

(a) 

0-246 

70 

4.5 

0-248 

6l 

4.6 

(a) 

0-249 

53 

5.2 

0-250 

58 

4.6 

0-251 

60 

5-4 

0-252 

64 

4.7 

0-253 

80 

5-4 

0-254 

69 

4.4 
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TABLE  4  (continued) 


Run  Number 

0 

(sec-^-) 

f 

(kc.) 

Notes 

0-255 

66 

5-2 

0-256 

72 

4.1 

0-257 

66 

5-3 

0-253 

57 

4.4 

0-260 

61 

5.3 

0-266 

53 

2.4 

0-267 

56 

2.6 

Notes : 

(a)  Marginal  mean  pressure,  the  mean  pressure  trace  indicated  a 
non-uniform  burning  area. 
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TABLE  5 

RELATIVE- AMPLITUDE  AMD  DAMPING  DATA  FOR  GAS -PHASE  STUDY 

Nominal  mean  pressure,  14.5  atm. 

Nominal  initial  sample  length,  1.3  cm. 


f 

fl 

$ 

L 

"aD 

Run  Number  (kc.) 

Pg 

(degrees ) 

(cm) 

(sec“l)  Notes 

0-195 

- 

- 

- 

- 

- 

(a,b) 

0-200 

3.78 

.505 

- 

10.2 

- 

(c) 

0-208 

3.76 

.481 

- 

11.8 

- 

0-209 

3.82 

.540 

- 

11.0 

- 

(c) 

0-210 

3.86 

.555 

- 

11.8 

- 

0-211 

4.01 

.526 

11.0 

- 

0-212 

3-70 

.610 

11.8 

- 

0-214 

3.70 

.425 

— 

11.0 

- 

0-217 

3.69 

.445 

- 

11.0 

78 

0-218 

2.66 

.445 

- 

15-4 

37 

0-219 

3.69 

.410 

178.5 

11.0 

75 

(c) 

0-220 

2.66 

.312 

178.5 

15.4 

54 

0-221 

3-42 

.424 

177.5 

11.6 

63 

0-223 

3.78 

.518 

177.0 

11.8 

47 

0-224 

3-90 

.453 

178.8 

11.1 

42 

0-225 

3.60 

.403 

177.7 

11.8 

59 

0-226 

3.90 

.465 

178.3 

11.0 

42 

(e) 

0-227 

1.42 

•  398 

178.6 

26.2 

- 

0-228 

1.80 

.449 

179-5 

21.1 

- 

0-230 

1.44 

.505 

178.5 

26.2 

- 

0-231 

1.39 

.440 

178.9 

26.2 

- 

0-232 

1.8l 

.425 

179.2 

21.1 

- 

0-233 

1.48 

.463 

179.2 

26.2 

- 

0-234 

1.85 

.446 

179-2 

21.1 

- 

0-235 

1.96 

.380 

179.1 

21.1 

- 

0-236 

1.45 

•  374 

178.8 

26.2 

- 

0-237 

1.46 

.427 

179.2 

26.2 

- 

0-238 

1.94 

.445 

178.9 

21.1 

- 

0-241 

1.43 

.448 

178.9 

26.2 

- 

0-242 

1.90 

.418 

179.2 

21.1 

- 

0-243 

4.36 

.525 

- 

10,5 

38 

M 

0-244 

3.51 

.474 

- 

11.4 

- 

(c) 

0-248 

3.51 

.455 

178.8 

11.2 

87 

(c) 

0-249 

4.25 

.613 

178.8 

10.3 

46 

0-250 

3.55 

.427 

179.0 

11.2 

52 
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TABLE  5  (continued) 


f 

m 

0 

L 

°D 

Run  Number 

(kc. ) 

ipd 

(degrees ) 

(cm) 

(sec”'*')  Notes 

0-251 

4.24 

.479 

179*0 

10.3 

69 

0-252 

3*50 

*455 

177*4 

11.2 

79 

0-253 

4.24 

.425 

177.9 

10.3 

43 

0-254 

3.60 

.467 

178.9 

11.2 

51 

0-255 

4.09 

.518 

178.2 

10.3 

63 

0-256 

3*52 

.465 

178.6 

11.2 

54 

0-257 

4.01 

.425 

178.4 

10.3 

46 

0-258 

3*40 

.585 

179.1 

11.2 

68 

0-260 

4.15 

*455 

- 

10.5 

48 

Notes : 

(a)  Omitted  run  numbers  correspond,  in  general,  to  runs  with  very- 
poor  mean  pressure  behavior,  including  0-195* 

(b)  No  reference  signal  was  used  before  0-219  and,  therefore, 
the  phase  angle  could  not  be  accurately  measured  for  these 
early  runs. 

(c)  Marginal  mean  pressure,  the  mean  pressure  trace  indicated  a 
non-uniform  burning  area. 

(d)  Unusually  low  amplitude. 

(e)  No  Ojj  was  measured  in  lower  frequency  range. 
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TABLE  6 

DATA  FROM  MEASUREMENTS  IN  THE  CURVED  WALL  OF  THE  MOTOR 

Nominal  mean  pressure,  14.5  atm. 

Nominal  initial  sample  length,  1-3  cm. 


Run  Number 

Position 
of  Wall 
Transducer  (a) 

Ratio  of 
Amplitudes  (b) 

f 

(kc.  ) 

L 

(cm) 

Notes 

264 

•  75 

1.0 

3.2 

16.5 

(c) 

265 

•  75 

1.0 

3-3 

16.5 

(c,f) 

266 

.94 

1.1 

2.4 

(d) 

.92 

1.2 

.88 

1.2 

267 

.80 

.60 

2.6 

16.5 

(c) 

.78 

•45 

.77 

.44 

.76 

•  34 

270 

.86 

.19 

1-7 

21.6 

(e) 

.85 

.18 

.84 

.18 

271 

.86 

.22 

1.8 

21.6 

(e,  f ) 

.85 

.21 

.82 

.17 

272 

.81 

.18 

2.50 

16.5 

(e) 

.80 

.17 

.78 

•1? 

.76 

.15 

296 

•  98 

.97 

2.7 

16.0 

(a) 

.96 

.91 

•  95 

.81 

.9* 

•  99 

.93 

•  91 

299 

.81 

.15 

2-5 

16.5 

(e,f) 

.79 

.07 

.76 

.08 

300 

.81 

.07 

2.5 

16.5 

(e,f) 

.79 

.05 

.76 

.06 
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Notes : 

(a)  Position  of  the  wall  transducer  expressed  as  a  fraction 
of  the  length  of  the  gas  cavity. 

(b)  Ratio  of  the  amplitude  measured  at  the  wall  to  that  measured 
at  the  end. 

(c)  Double-ended  firing. 

(d)  Single-ended  firing  with  end  measurement  under  the  propellant. 

(e)  Single-ended  firing  with  end  measurement  at  the  non-burning 
end. 

(f )  Mean  pressure  trace  indicated  that  the  burning  area  was 
irregular. 
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TABLE  7 

DATA  FROM  LOW-FREQUENCY  FIRINGS 

Nominal  mean  pressure,  35  &tm. 

Nominal  initial  sample  length,  1.3  cm. 


Run  Number 

-a  . 

( s  ec  ~  ^ ) 

f 

(cps) 

Initial 

Pressure 

(atm) 

Igniter  only 

0-285 

20.0 

180 

33.4 

0-286 

20.1 

175 

35.2 

0-287 

19.4 

184 

40.3 

Utah-F 

Single-ended 

0-283 

8.9 

244 

32.2 

SRI -.109 

0-281 

11.2 

220 

31.8 

Single-ended 

0-282 

8.9 

230 

32.5 

0-290 

15.4 

210 

32.2 

double-ended 

0-284 

4.4 

260 

32.8 

0-288 

4.2 

222 

31-5 

0-289 

2.9 

- 

32.2 

SRI-111 

Single-ended 

0-291 

8.3 

202 

32.2 

0-292 

9.4 

200 

32.2 

0-294 

8.7 

214 

32.5 

0-295 

7.4 

20 6 

32.2 

SRI -111 

Double-ended 

0-293 

3-6 

358 

32.5 
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TABLE  8 


DATA  FROM  SOLID  PHASE  WORK 

Nominal  mean  pressure,  14.5  atm. 

Run  Number 

l 

(cm) 

Pal 

SI  . 

♦ 

(degrees) 

f 

(kc) 

Notes 

159 

4.0 

•53 

165 

3.9 

(a) 

3.2 

•53 

166 

- 

2.4 

.54 

168 

- 

1.6 

•57 

168 

- 

.8 

•57 

168 

- 

162 

5.6 

1.0 

160 

3-9 

(a) 

4.8 

.94 

l6l 

- 

4.0 

.91 

159 

- 

3.2 

1.1 

129 

- 

2.3 

.84 

129 

- 

1.7 

.62 

150 

- 

.9 

.87 

139 

- 

168 

9-8 

1.0 

142 

3-9 

(a) 

8.9 

.87 

149 

8.1 

.81 

152 

- 

7-2 

.84 

153 

- 

6.4 

.84 

150 

- 

5-5 

•  70 

156 

- 

4.7 

.60 

162 

•  - 

3-8 

.50 

162 

- 

3.0 

•  57 

166 

- 

2.1 

.58 

166 

- 

1-3 

•  71 

160 

- 

.4 

•  71 

167 

- 

169 

15.3 

1.2 

38 

3.81 

(b) 

14.4 

.83 

60 

3.66 

13.5 

.70 

95 

3.47 

12.6 

•  63 

126 

3.42 

11.7 

•  77 

139 

-  r,2 

10.8 

•  97 

156 

1 

9-9 

1.3 

169 

0 

9.0 

1.4 

170 

Z^3 

8.2 

1.5 

176 

3.44 

7-3 

1.5 

180 

3.45 

6.4 

1.7 

178 

3.46 

5-5 

1.7 

176 

3.56 

4.6 

1.8 

174 

3.62 

3-7 

1.8 

175 

3.61 

2.8 

1.8 

174 

3-55 

-  171a  - 


TABLE  8  (continued) 


Run  Number 

i 

(cm) 

lpsl 

1*3 

♦ 

(degrees) 

f 

(kc) 

Notes 

0-177 

13-3 

1.4 

60 

3.80 

(b) 

12.5 

1.4 

62 

3-71 

11.7 

1.7 

94 

3-77 

10.9 

1.5 

112 

3.66 

10.1 

1.2 

128 

3.63 

9-3 

1.2 

132 

3.60 

8.6 

1.1 

136 

3.72 

7.8 

1.1 

133 

3.87 

7.0 

•  93 

143 

3.90 

6.2 

.87 

151 

3.82 

5-4 

.98 

150 

3.89 

4.6 

.92 

151 

4.10 

3-8 

.83 

157 

4.07 

3.1 

.76 

160 

4.02 

2.3 

.69 

163 

4.03 

1.5 

.68 

164 

4.01 

.7 

.64 

167 

4.00 

0-179 

14.6 

1.3 

43 

3.47 

(b) 

14.2 

1.7 

59 

3.42 

13-8 

1.8 

62 

3.44 

13-4 

1.9 

54 

3.59 

12.9 

1.4 

29 

3.68 

12.5 

1.9 

49 

3.76 

12.3 

1.6 

36 

3.81 

11.7 

1.7 

36 

4.20 

10.8 

1.8 

37 

4.30 

10.0 

2.2 

58 

4.31 

9.1 

2.8 

67 

4.29 

8.3 

3.5 

105 

4.26 

7-4 

2.5 

126 

4.27 

6.7 

2.7 

136 

4.29 

5-9 

2.1 

150 

4.31 

5-1 

1.7 

160 

4.33 

4.2 

.1 .4 

166 

4.33 

3-3 

1.3 

165 

4.33 

0-182 

10.7 

•  71 

46 

4.58 

(b) 

9.8 

1.1 

31 

4.52 

9.0 

1.6 

4l 

4.46 

8.2 

1.6 

53 

4.46 

7.4 

2.7 

90 

4.45 

6.6 

2.4 

126 

4.38 

5-7 

1.9 

143 

4.37 

T 
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TABLE  8  (continued) 


Run  Number 


0-184 


0-187 


O-189 


i 

(cm) 

We  1 

<J> 

(degrees) 

f 

(kc) 

4.9 

1.6 

153 

4.36 

4.0 

1.4 

155 

4.35 

3.2 

l.l 

162 

4.34 

2.4 

1.0 

164 

4.35 

2.0 

•  97 

j  64 

4.36 

1.5 

•  95 

lo2 

4-33 

1.1 

1.0 

161 

4.32 

.7 

1.0 

166 

4.30 

15. c 

.96 

44 

4.10 

14.5 

.77 

30 

4.05 

13-5 

•  77 

38 

3.90 

13.0 

.88 

43 

3.95 

12.1 

.88 

46 

4.06 

11.2 

.76 

58 

4.11 

10.2 

•89 

72 

4.13 

9.3 

.92 

90 

4.11 

8.4 

1.0 

104 

4.05 

7-5 

1.2 

138 

3*97 

7-1 

1.1 

127 

3-98 

6.5 

1.1 

136 

4.04 

6.2 

1.1 

128 

4.09 

5-2 

.96 

136 

4.19 

4-3 

.92 

146 

4.22 

3.4 

•  94 

161 

4.21 

2.5 

•  77 

168 

4.20 

1-5 

.65 

164 

4.24 

12.7 

1.4 

44 

4.0 

12.1 

1.7 

51 

11.4 

\  .6 

47 

10.7 

1.5 

54 

10.0 

1.5 

66 

9.4 

1.1 

60 

14. i 

1.2 

4l 

3.9 

13.3 

1.4 

73 

12.6 

1.6 

96 

11.8 

1.4 

122 

11.1 

1.2 

134 

10.3 

1.0 

140 

9.6 

.92 

148 

8.8 

.86 

151 

8.1 

1.0 

154 

Notes 


(b) 


(a) 


(a) 


Run 


0- 


0- 
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TABLE  8  (continued) 


Run  Number 

l 

(cm) 

15=1 

lr6l 

4> 

(degrees) 

f 

(kc) 

Notes 

7-3 

•95 

157 

6.6 

•98 

160 

5.8 

•92 

160 

5-1 

.89 

159 

4.3 

.82 

161 

3.6 

.86 

162 

2.8 

.81 

162 

2.1 

.82 

163 

1.4 

.83 

160 

.6 

•  73 

l6l 

0-191 

14.0 

.82 

21 

4.0 

(a) 

13-3 

l.l 

27 

12.6 

1.4 

36 

11.9 

2.1 

48 

11.2 

2-3 

79 

10.5 

2.5 

95 

9-8 

2.2 

96 

9.1 

2.0 

98 

8.4 

2.0 

92 

7-7 

1.8 

90 

7-0 

1.7 

ill 

6-3 

1.6 

132 

5.6 

1.6 

144 

4.9 

1.1 

l6l 

4.2 

1.0 

165 

3-5 

•  93 

166 

2.8 

•  93 

168 

2.1 

•93 

169 

1.4 

.90 

173 

•  7 

.88 

172 

0-194 

3-7 

•  58 

166 

3-4 

(a) 

2.9 

•  59 

164 

2.1 

•  58 

166 

1.3 

•  57 

166 

•  5 

.63 

167 
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Notes : 

(a)  Nominal  frequency  for  entire  run. 

(b)  Frequencies  tabulated  were  obtained  from  a  plot  of  frequency 
measurements.  These  measurements  were  made  by  sampling  the 
recorded  data  at  equal  time  increments.  The  tabulated  values 
are  accurate  to  two  places  only;  however,  a  third  place 

was  retained  to  indicate  the  frequency  control. 
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Roman  (cont’d) 


£ 

energy  density 

o 

ergs /cm3 

f 

frequency 

cycles/sec 

f(x,y,z) 

unspecified  function  of  position 

-- 

h 

solid-phase  work,  defined  by 
Equation  (6.9/ 

cm-1 

h 

Appendix  A,  constant 

-1 

sec 

h 

Appendix  B,  C^T  +  p/p 

ergs/gm 

lm[  ] 

imaginary  part 

-- 

i 

indexing  variable 

-- 

Jn 

Bessel  function  of  first  kind 
and  order  x> 

dimensionless 

J 

(-d1/s 

dimensionless 

K 

arbitrary  constant  in 

Equation  (F.3) 

dimensionless 

k 

gas  phase,  thermal  conductivity 

cal/cm  sec 

k 

solid  phase,  defined  by 

Equation  (6.9) 

-1 

cm 

defined  by  Equations  (A. 19/, 

(A. 20),  and  (A.2l) 

-1 

cm 

L 

length  of  gas  cavity 

cm 

l 

length  of  propellant 

cm 

M 

Mach  Number 

dimensionless 

— * 

m 

mass -flux  vector 

gm/crn^  sec 

m 

mass  flux 

gm/cm^  sec 

n 

integer 

dimensionless 

-* 

n 

unit-normal  vector 

dimensionless 

P 

defined  by  Equation  (A. 6) 

2  /  2 
cm  /sec 

-  l80  - 


Prandtl  Number,  pC^/lc 

total  pressure 

time-average  pressure 

real  part  of  first-order 
perturbation  in  pressure 


p  complex,  first-order 

perturbation  in  pressure 


oscillating  pressure  at  the  cold 
end  and  at  center  line 


p  oscillating  pressure  under  the 

propellant  and  at  center  line 

p^.  oscillating  pressure  at  the 

burning  Interface  and  center  line 


-[] 


solutions  to  equations  of  the 
type  of  Equation  (A. 9) 

real  part 


r 


radial  position 


unperturbed  burning  rate 

r  at  the  wall 

surface  area 

condensation  P*/Pq 

total  temperature 

perturbed  temperature,  first 
order  and  complex 

time-average  temperature 

time-average  temperature  in 
main  body  of  gas 

vail  temperature 


dimensionless 

atm 

atm 

atm 

atm 

atm 

atm 

atm 

dimensionless 

dimensionless 

cm 

cm/sec 

cm 

2 

cm 

dimensionless 

K 

K 

K 

K 

K 
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Roman  (cont'd) 


T 

Ot 

mean  temperature  at  throat  of  nozzle 

K 

t 

time 

sec 

u 

solid  phase,  defined  by 

Equation  (6.7) 

cm 

u 

Appendix  B,  C^T 

erg/gm 

VU2 

Appendix  B,  first-  and  second- 
order  perturbations  of  U 

erg/gm 

u 

axial  velocity,  total 

cm/sec 

uo 

axial  velocity,  time  average 

cm/sec 

u 

axial  velocity,  perturbed  and 
complex 

cm/s  ec 

u 

velocity  vector,  total 

cm/sec 

“o 

velocity  vector,  time  average 

cm/sec 

— # 

U1 

velocity  vector,  real  part  of 
first-order  perturbation 

cm/sec 

"2 

velocity  vector,  real  part  of 
second-order  perturbation 

cm/sec 

V 

volume 

3 

cmJ 

V 

radial  velocity,  perturbed 
and  complex 

cm/sec 

W 

defined  by  Equation  (6.7) 

cm 

X 

position  coordinate 

cm 

y 

position  coordinate 

cm 

y 

acoustic  admittance,  complex 

cm'Vdyne  sec 

y' 

real  part  of  admittance,  y 

cm'Vdune  sec 

y" 

imaginary  part  of  admittance,  y 

cnr/dyne  sec 

z 

position  coordinate 
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Greek 


a 

growth  constant 

sec  ^ 

ai 

single-ended  a 

-1 

sec 

a2 

double-ended  a 

-1 

sec 

ab 

oc  for  burning  surface 

sec  ^ 

aD 

initial  decay  constant 

-1 

sec 

aDL 

linear  decay  constant 

sec  ^ 

Ct(A) 

amplitude  dependent  decay  constant 

-1 

sec 

“d 

decay  constant 

-1 

sec 

complex  angular  frequency,  cd-JO 

-1 

sec 

r 

arbitrary  constant,  Equation  (F.2) 

-1 

cm 

7 

propagation  constant,  complex, 
b-ja 

-1 

cm 

7t 

ratio  of  heat  capacities,  C^/C^ 

dimensionless 

A 

defined  by  Equation  (6.3) 

dimensionless 

6 

arbitrary  constant,  Equation  (4.5) 

dimensionless 

€ 

reduced  pressure  perturbation, 
p/p0 

dimensionless 

defined  by  Equation  (F.4) 

dimensionless 

V 

defined  by  Equation  (F.IO) 

dimensionless 

1 

T) 

0 

defined  by  Equation  (A. 31 ) 
lP0Cv\  - 

I — —  T,  defined  by 

\  po 

Equation  (A. 4b) 

X 

sec  2 

dimensionless 

9(x,y,z) 

position-dependent  phase  angle 

dimensionless 

X 

solid  phase,  Lam^  constant 

dyne/ cm^ 

X 

Appendix  A,  eigen  value 

-2 

cm 

V*2 

roots  of  Equation  (A.IO) 

-2 

cm 

*"  T* . 


-<•  - 


V 
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p 

solid  phase,  Lame  constant 

dyne/ cm^ 

p 

gas  phase,  viscosity 

gm/cm-sec 

p' 

**/P0 

2  , 

cm  /sec 

p" 

Hi 

3 

2  , 

cm  / sec 

“p 

reduced  mass-flux  perturbation 
m/m0 

dimensionless 

X) 

,  defined  by  Equation  (A.  4c) 

p0  v 

2  / 

cm  / sec 

axled  displacement 

cm 

5r 

redial  displacement 

cm 

Jf 

3.14159 

dimens _onl ess 

P 

total  density 

o 

gm/ cm3 

p0 

time  average  density 

gm/ cm^ 

P 

perturbed  density,  first  order 
and  complex 

gm/cm^ 

P1 

perturbed  density,  first  order 
and  real  part 

o 

gm/cm3 

P2 

perturbed  density,  second  order 
and  reed  part 

gm/cm^ 

Pq  at  wedl  conditions 

gm/cm^ 

P8 

density  of  solid 

gm/ cm^ 

xx  rx'  rr 

stress  components 

dyne/ cm^ 

♦ 

phase  angle  between  jTg  and  p”g 

dimensionless 

-1 

03 

angular  frequency,  real  part 
of  0 

sec 

“e 

defined  by  Equation  (6.4) 

dimensionless 
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